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Absolute intensity measurements of the transitions in ammonia vapor with maxima at 
52,500, 66,000, and 77,000 cm™ are reported. Calculations are made of the intensities of the 
first and third transitions, using non-bonding nitrogen orbitals, of the form suggested by Slater. 
The calculated values are considerably lower than the measured values, but may be improved . 
by variation of the effective nuclear charge. Similar calculations for the second transition are 
less conclusive, but seem to indicate that it also may originate from a nearly non-bonding 
nitrogen orbital. Some other possible origins of this transition are considered and rejected on 
the basis of either their energy or their intensity. 





INTRODUCTION 


REVIOUS investigation of the absorption 

spectrum of ammonia vapor"? has revealed 
the presence of four electronic levels in the 
frequency range 45,000—120,000 cm~. We report 
here absolute intensity measurements of transi- 
tions from the normal level to the two lowest of 
these electronic levels, designated as transitions 
I and II, with intensity maxima at about 52,500 
and 66,000 cm. In addition, the intensity of 
the transition III with maximum at about 77,000 
cm~! may be estimated by comparison with I and 
II using previously obtained appearance-pressure 
data. There appears to be little hope of obtaining 
a direct measurement of the intensity of III or 
IV with present experimental techniques. 

As is well known, an experimental value of 
the absolute intensity of absorption associated 
with an electronic transition in a pure vapor is 

* Part of a dissertation presented to the Faculty of the 
Graduate School of the University of Rochester in partial 
fulfillment of the requirements for the degree of Doctor 
of “gual 


1A. B. F. Duncan, Phys. Rev. 47, 822 (1935). 
? A. B. F. Duncan, Phys. Rev. 50, 700 (1936). 


found by evaluation of /{a,dv over the entire 
absorption region. At a single frequency, a, 
=(1/lo) log. (Io/I), where Io is the intensity of 
the incident light and J is the intensity after 
passage through a column of the absorbing gas 
of length J) measured at 0°C and 760-mm pres- 
sure. The light source used should emit sharp, 
well-resolved lines of known frequency, well dis- 
tributed over the wave-length region in question. 

Values of a, in cm~ are plotted against wave 
number, and the value of /a,dv determined 
with a precision planimeter. The value of the 
integral multiplied by mc*/nme?(=4.2 X 10-* cm?) 
gives the f number or oscillator strength of the 
transition. In order to compare experimental 
intensities with calculated values, we may com- 
pute the square of the dipole moment integral 
(in square angstroms) from the approximate 
relation 


O?(exptl) =(3hf/82?vmc) X 10'*(9.12 XK 104f/¥m), 


where v», is the wave number of the intensity 
maximum of the transition. 


573 





574 ie & 


Two methods of photographic photometry 
may be used to evaluate J)/I at a single known 
frequency from measurements of photographic 
density. In both methods the intensity of the 
source must be kept as constant as possible. In 
the first method (intensity scale method) the 
intensity through the absorption cell is varied 
by a system of calibrated diaphragms. On each 
plate a series of exposures for equal time inter- 
vals is made through each of the diaphragms 
with the cell evacuated, then one or more expo- 
sures are made through one of:the diaphragms 
with the cell filled with vapor at a measured 
pressure. 

For every line within the absorption band, a 
series of density values corresponding to the 
series of exposures is recorded with a Moll type 
microphotometer. From a plot of density against 
log 10 (Jo/Z) for each line, [o/J for the exposure 
through the vapor, and hence a,, may be ob- 
tained directly. 

In the second method (time scale method) the 
intensity through the evacuated cell is kept 
constant and a similar series of exposures is made 
for different time intervals. In this case, density 
—log time curves give logio (¢o/t), which is equiv- 
alent to logyo (Jo/I) by the usual photographic 
reciprocity law. The exposure through the vapor 
is then made for time fo. In this work we have 
used both methods. 

The general method outlined in the preceding 
five paragraphs should give correct intensity 
values in the case of strictly continuous spectra. 
While the absorption spectrum of ammonia is 
not truly continuous, it approximates continuity 
sufficiently under the experimental conditions 
used here, since a, does not change anywhere 
very abruptly with frequency. At considerably 
lower pressures than were used in the present 
work and with slit widths 0.008-0.015 mm, 
vibrational structure is easily resolved, and in 
transition II, rotational structure is partly re- 
solved.' In the present work, with slit 0.2 mm 
wide, and with pressures 5-10 times that used 
to resolve fine structure, vibrational structure 
can be recognized as broad low maxima in the 
absorption curves for transition II. The fre- 
quencies, at which a, were measured, are suffi- 
ciently close together to show the detailed shapes 
of these curves. Not even vibrational maxima 
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are clearly recognized in the curves for transition 
I, and no rotational structure is apparent in any 
case. 

The integrated absorption coefficients for tran- 
sitions I and II are independent of the pressure 
in the range used, or rather, any pressure de- 
pendence is concealed and included in the total 
experimental error in the f numbers (5-10 
percent). 


EXPERIMENTAL 


The ammonia gas used in this study was 
purified by a large number of distillations in 
vacuum through KOH pellets, and subsequent 
fractionation. It was dissolved in dried NH,NO; 
and NH,SCN in the molar proportions 1:3 for 
storage, the resulting solution giving a convenient 
NH; pressure of about 180 mm. The desired 
pressures were obtained by expansion through 
calibrated volumes, and the pressure in the cell 
was read with a McLeod gauge. 

The apparatus used to study transitions | and 
Il was in all essentials the same as that of 
Ladenburg and Van Voorhis* and need not be 
described in detail. The background for observa- 
tion of absorption was a discharge in pure CO, 
for the region 45,000 to 61,000 cm (transition 
1) and in pure He for the region 60,000 to 71,000 
(transition I]). In the former case the well-known 
A'II —X'Z bands of CO appeared with great in- 
tensity and the heads were sufficiently sharp to 
serve as lines for photometric measurement. The 
wave-lengths given by D. N. Read‘ were used. 
In the second case the rotational lines in the He 
bands of the Lyman ('Z,+—'!Z,*) system® were 
used. 

The discharge tube was operated at a constant 
direct current of 0.065 ampere through a voltage 
regulator and a transformer-rectifier unit. Less 
than 0.001-ampere variation in current was ob- 
served with the longest exposures. The light 
from this discharge was made parallel by a LiF 
lens, passed through an absorption cell 17.7 cm 
long, and was focused on the slit of the vacuum 
spectrograph by a second LiF lens. The spectra 
were photographed from a 2-meter radius, 

3R. Ladenburg and C. C. Van Voorhis, Phys. Rev. 43, 
315 (1933). 


4D. N. Read, Phys. Rev. 46, 571 (1934). 
5 T. Lyman, Astrophys. J. 23, 181 (1906). 
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ELECTRONIC TRANSITIONS IN AMMONIA 


TABLE I. 








Frequency 
maximum Assumed 
Transition (cm!) type 


Q/ao 


obs 





I 52500 2p.—3s 
I] 66000 2p.—3d. 
il 77000 2p.—4s 
IV 85000 2p.—5s 


0.90 0.0134 
0.22 0.0017 
0.52 0.0001 








120,000-line grating, used in the first order.* The 
dispersion was about 4.15A per mm. The expo- 
sure times were 40 to 60 minutes. Eastman III-O 
and I-O ultraviolet sensitized plates were used, 
and were processed with the special precautions 
required in photographic photometry. 

The decrease in effective transmission of the 
LiF lenses became quite appreciable below 
1400A (70,000 cm) when used in the optical 
system described. To observe transition II] 
(and IV, with maximum at 85,000 cm) would 
have required prohibitively long exposures. We 
considered the possibility of using the spectro- 
graph as an absorbing column and moving the 
discharge close to the slit, with a very thin LiF 
plate separating the two. However, we did not 
believe that the pressure of ammonia could be 
maintained sufficiently constant in the spectro- 
graph, especially at the low pressures required. 


RESULTS 


Seven plates, four on the intensity and three 
on the time scale, were used to obtain fa,dv for 
transition I. The pressure of ammonia ranged 
from 0.150 to 0.875 mm, the most useful range 
being 0.2 to 0.5 mm. For transition II eight 
plates on the intensity scale and one on the time 
scale were used. The pressure range was 0.3 to 
3.0 mm, the most useful range being between 1.0 
and 2.0 mm. The results from the intensity and 
time scales checked within the experimental 
accuracy. Reciprocity law failure was negligible 
in these experiments. 

Average experimental values of the dipole 
strength Q and the f number of transitions I, I], 
and III are given in Table I, together with two 
sets of calculated values for transitions I, II, II], 


and IV. These calculations will be discussed in. 


the next section. do is 0.528A, the Bohr radius 
of the normal hydrogen atom. 


——. 


°A. B. F. Duncan, J. Chem. Phys. 8, 444 (1940). 


The experimental results for transitions I and 
II are estimated to be accurate to within 5 to 10 
percent. The indirect values for transition III 
are subject to somewhat larger error, although 
this is not necessarily much larger. The probable 
accuracy of these values is difficult to estimate. 


DISCUSSION OF RESULTS 


The electronic configuration of the ground 
state of ammonia, omitting the 1s electrons of N, 
is given by Mulliken’ ® as (sa,)*(me)*(za1)". The 
orbital sa; is primarily a N atom °2s orbital 
which, however, contributes a little to the three 
N—H bonds. The two ze orbitals are strongly 
bonding. The orbital za; is non-bonding and is 
approximately just a 2p, orbital of the N atom. 
Mulliken suggests that the lower excited states 
correspond to the excitation of one 2a; electron 
from this configuration, and that the lowest 
electronic transition should be 2za;—3sqa,. 


Dipole Moment Integrals for Transitions 
I, III, and IV°® 


We have assumed that the transitions (2a) 
—(nsa;) in ammonia may be approximated, for 
the purpose of intensity calculations, by the 
transitions 2p.—ns in the N atom. The dipole 
moment integrals will be of the form 


Q a Q. = (2) f vencbnadr 


and may be evaluated in spherical polar coordi- 
nates, with z=r cos 6. The factor (2)! allows for 
the fact that there are two equivalent 2), elec- 
trons, either of which may make the transition. 
The functions y2p, and y,, may be approximated 

7™R. S. Mulliken, J. Chem. Phys. 1, 492 (1933). 

8R.S. Mulliken, J. Chem. Phys. 3, 506 (1935). 

® The general method used here is discussed by R. S. 


Mulliken, J. Chem. Phys. 7, 14 (1939); Reports on Progress 
in Physics (1941), Vol. VIII. 
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by the radial functions given by Slater,!® multi- 
plied by the appropriate angular factors. The 
normalized radial parts of the functions are: 


(2Z*/n*)"+8 


Rulr) = 
ao\(I'(2n*+1))! 





(r/ao)"— 


Xexp [—(2*r/n*ap) ]. 


The complete 2p, function is obtained by multi- 
plying R2,(r) by (6) cos 6/2(27)!; the complete 
¥ns function is obtained by multiplying R,,. by 
(2)!/2(22)}. ao is the radius of the Bohr orbit on 
normal hydrogen. Substitution of these functions 
and integration yields 
Q (221*/my*)""+4(2Z2* /m2*) 0244 
= = (2)8(31/3)| ered 
do [1 (2n,*+1)0(2n.*+2) }} 

| I'(ni* +n2* +2) | (1) 

[(Za*/mr*) + (Za*/nat) Jorrtner st } 








The factor (3)!/3 is the result of integration 
over the angles @ and g. Subscripts 1 and 2 refer 
to the lower (2f,) and upper (ms) states. 

The values of the parameters Z* and n* in the 
states considered are given by Slater as: 


State 2p 3s 4s 5s 
zZ* 3.9 1.6 1.0 1.0 
n* 2.0 3.0 3.7 4.0 


The values of Q/ao calculated with these values 
of the parameters are listed in Table I. The 
values of f were calculated for the transitions 
2p.—mns on the assumption that these correspond 
to the observed transitions I, III, and IV, which 
were found earlier? to fit into a Rydberg series. 
Table I shows that the experimental values of 
the dipole moment integrals associated with 
transitions I and III and the values calculated 
by means of Slater functions differ by factors of 
3 and 26, respectively. The experimental value 
of the mean oscillator strength f (which is pro- 
portional to the intensity of the transition) is 10 
times the calculated value for transition I and 
550 times the calculated value for transition III. 
This disagreement between theory and experi- 
ment may be considered from two viewpoints. 
First, the electron involved in the transitions I 
and III may not be completely non-bonding. 


10 J. C. Slater, Phys. Rev. 36, 57 (1930). 
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Second, the functions used in the numerical 
evaluation of the dipole moment integrals are 
probably crude approximations to the true func- 
tions. Even if the form of the wave functions is 
approximately correct, the appropriate values 
of the parameters Z* and n* for our use may 
not be those given by Slater. 

The first viewpoint suggests that improvement 
in the calculated values may result from admix- 
ture of the orbitals 2s(N), 2p.(N), and 1s(H) for 
the normal state wave function (as suggested by 
Mulliken’). Perhaps the upper states in the 
transitions should not be regarded as pure N 
atom states and a linear combination of orbitals 
should be used there also. Calculations employing 
such mixed functions appear to offer considerable 
difficulty and may be inconclusive in view of the 
undetermined coefficients appearing in these 
linear combinations. We have not made such 
calculations. 

The second viewpoint suggests that we might 
vary Z* in the Slater function to obtain better 
agreement between calculation and experiment. 
The variation of Slater’s Z* has been discussed 
by Mulliken."' We have varied the Z* in Eq. (1) 
and find that if Z:,* is to be near the Slater 
value 3.9, then Z;,* must be increased to about 
3.5 to get agreement with experiment in transi- 
tion I. This is indeed a large increase, in com- 
parison with the increases suggested by Mulliken 
in other cases. 

This large increase in Z3,* appears improbable 
on the basis of simple energy considerations, 
such as have been used by Slater. Using the 
method of reference 10, the lowest ionization 
potential of ammonia, taken as the energy differ- 
ence of the configurations 1s?2s?2p?— 15s?2s?2p’ is 
close to the experimental value (11.56 ev obs., 
12.88 calc.). Similarly, with Z3,* = 1.6, the Slater 
value, the difference between 15?2s?2p? and 
15?2s?2p?3s is 3.85 ev (obs. in ammonia, 5.05). 
An increase of Z;,* to 1.8 gives exact agreement 
with the observed term value, while an increase 
to 3.5 gives a term value greater than the normal 
state. It is uncertain how seriously these energy 
calculations should be considered. 

It is well known that the upper level of 

4R.S. Mulliken, J. Chem. Phys. 8, 234 (1940); Reports 


on Progress in Physics (1941), Vol. VIII, pages 244, 250, 
269. 
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transitions I in ammonia is predissociated. The 
origin of this predissociation has not been exam- 
ined theoretically in detail. If the perturbation 
causing the predissociation is merely of a vibra- 
tional nature, it should have no effect in a 
comparison of experimental and theoretical in- 
tensities. If the perturbation is caused by another 
electronic level, the experimental intensity may 
well be higher than the intensity calculated for a 
single, atom-like transition. It appears at least 
that no simple explanation can be offered at the 
present time for the intensity of transition I. 


Discussion of Transition II 


It was suggested previously? that this transi- 
tion might be to a triplet level associated with 
the upper level of transition II]. From the 
present quantitative measurements, the experi- 
mental ratio of intensities of III to II is about 
7:1, which appears to be too low for a singlet- 
triplet ratio. In some previous unreported work, 
we looked for transitions at frequencies lower 
than transition I, in an absorption tube 10 
meters long with ammonia at one atmosphere 
pressure. None was found, indicating a very low 
intensity for any possible triplet transition asso- 
ciated with I. We would expect the triplet 
associated with III to have an even lower in- 
tensity. Accordingly, we withdraw this previous 
suggestion. 

We will discuss briefly other possible explana- 
tions for this transition. Discussion will be 
limited to parallel transitions, since analysis of 
the rotational structure of the bands of transition 
I] showed it to be of this type.” 

Of the other atom-like transitions of a za, 
electron we considered 22a; — 3da;, which we have 
assumed may be approximated by the N atom 
2p.—3d. transition. The dipole moment integral 


© A. B. F. Duncan and G. R. Harrison, Phys. Rev. 49, 
211 (1936). 


was calculated in the same way as before. The 
value of the integral is given by Eq. (1), except 
that integration over the angles now gives 
2(15)3/15. The calculation was made using 
Slater’s value of Z2,*=3.9, Z3a*=1.0, m2p* =2.0, 
and m3q*=3.0. The experimental and calculated 
results are compared in Table I. If Z3a* is 
increased from 1.0 to 1.5, complete agreement 
with experiment is obtained, without adjustment 
of any other parameters. By analogy with the 
energy levels of the nitrogen atom, this transition 
should occur at a higher frequency (about 80,000 
cm) in ammonia. However, as Mulliken points 
out,® the lower members of the Rydberg series 
in a polyatomic molecule may show abnormalities 
because of their relationship to orbitals in the 
core. For example, the frequency of transition | 
in ammonia is about 13,000 cm lower than 
would be expected by analogy with the corre- 
sponding transition in nitrogen. The transition 
2p.—3d, should therefore be considered as a 
possible explanation for transition II. 

Transition Il may, however, be caused by 
excitation of another electron in ammonia. The 
estimated ionization potential of an sa; electron 
is given as 27 volts,’ which indicates that even 
the first member of the (sa;—za,) series should 
lie above the ionization potential of the za, 
electron, and should not appear in the spectral 
region investigated. The ionization potential of 
a me bonding electron is estimated to be about 
16 volts, so it is quite possible that the first 
transition due to this electron may be in the 
region of transition II. The possible transitions 
consistent with parallel rotational structure are 
the N—V, bonding to anti-bonding transition 
and the (xe)—(ze)’ bonding to bonding transi- 
tions in which the upper state orbitals may 
possess atom-like characteristics. It should be 
interesting to attempt calculations of the in- 
tensities of such transitions in ammonia. 





THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 14, 


NUMBER 10 


OCTOBER, 1946 


The Experimental Determination of the Intensities of Infra-Red Absorption Bands 


I. Theory of the Method* 


E. BriGHT WILSON, JR. AND A. J. WELLS** 
Mallinckrodt Chemical Laboratory, Harvard University, Cambridge, Massachusetts 


(June 28, 1946) 


With spectrographs of available resolving power, the apparent integrated absorption coef- 
ficients of infra-red bands usually differ greatly from the true coefficients because the spectro- 
graph yields at each setting a weighted average of the fractional light transmitted in a band 
of frequencies, whereas what is desired is the unweighted average of the logarithm of the frac- 
tional light transmitted. It is shown that true absorption coefficients can be obtained by (a) 
eliminating the violent fluctuations in intensity with frequency by broadening the rotational 
lines with a non-absorbing foreign gas and (b) eliminating the error due to the intensity varia- 
tion of the envelope by extrapolating the apparent integrated absorption coefficient divided 
by the partial pressure to zero partial pressure of the absorbing gas. These two steps permit 
vibrational intensities to be measured to a reasonable accuracy even with a spectrograph of low 


resolving power. 





EW accurate measurements on the intensities 

of infra-red absorption bands have been 
published in spite of the utility of such data. 
Many of the results which have been obtained 
have been in error by as much as several hundred 
percent. This error is mainly caused by the lack 
of adequate resolving power in the spectrographs 
used. In this paper a method will be described 
which permits fairly accurate results to be ob- 
tained on most substances, even with a spectro- 
graph of low resolving power. This method is 
primarily useful in handling the vibration-rota- 
tion bands of gases whose rotational spacing is 
not too wide. 


I, DEFINITIONS 


It will be assumed that the absorption of 
strictly monochromatic light is governed by the 


familiar law 
I=I,exp (—apL), (1) 


where J is the light transmitted by a cell of 
length L containing a gas at pressure p, and Io 
is the incident intensity. @ is the absorption 
coefficient divided by the pressure and is as- 
sumed to be independent of p and L but de- 
pendent on the frequency.' The quantity which 


* Much of the material in this paper is taken from a 
thesis by A. J. Wells, submitted to the Department of 
Chemistry, Harvard University, in May, 1941. A short 
note was published in J. Chem. Phys. 9, 659 (1941). 

** Now at Experimental Station, E. I. du Pont de 
Nemours and Company, Wilmington, Delaware. 

1 Actually a@ is not strictly independent of p due to self- 
broadening, but this does not affect the later arguments, 
when the total pressure, including that of a foreign gas, 
is held constant. 
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is desired for comparison with theory is the 
integrated absorption coefficient for unit pressure 


A= fadv=(1/pL) fin (Io/I)dv, (2) 


where the integration is carried out over the 
frequency range covered by a given band. If a 
spectrograph of infinite resolving power were 
available, Io/I could be measured directly and 
there would be no problem. 

No actual spectrograph measures J or J for a 
single frequency but instead measures the inte- 
gral of the intensity in a range of frequencies. 
Let T be the apparent intensity registered by the 
instrument at a setting v’ (representing the 
central frequency admitted by the finite slit). 
Then 


T= f 1a v’)dv, (3) 


where g(v, v’) is the fraction of the light of actual 
frequency v which is registered by the instrument 
when its setting is v’. The integration is carried 
out over all values of v for which g(y, v’) #0. 
Since g will diminish rapidly outside of a narrow 
band centered at v’, the limits of integration can 
be extended to + without harm. A similar 
expression holds for JT» in which J» replaces /. 
It is convenient to define an apparent absorption 
coefficient 6 for frequency v’ by the equation 


pB=(1/L) In (To/T), (4) 
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and, for each band, an integrated apparent 
absorption coefficient B for unit pressure 


B= [ aav'=(1/pt) f In (To/T)dv’. (5) 
T and B are measurable quantities from which 


A is to be obtained. 


II. EXTRAPOLATION THEOREM 
It is desired to prove that 


Lim B=A. 


pl—0 


To do so examine the difference 


A—-B=(1 ‘ptf In (TIo/Tol)dv’ 


, Sf(v)glv, v’)dv 
=(1 pL) { In —— —dv', (6) 
F(X S gly, v'dv 





where 
f(v) =exp (—apL) =I/Ib, (7) 


and J) has been assumed to be constant over a 
slit width.? 

lf f were constant over a slit width, it could 
be removed from the inner integral and the 
expression would vanish, so that, for resolving 
power large compared with the variations of 
exp (—apL), A and B are equal. But in practice 
this exponential will go rapidly from nearly zero 
to nearly one through each rotation line, if there 
is enough gas to make the average absorption 
appreciable. Under these conditions A and B 
will diverge widely. 

Now consider the limit as pL—0. Differentiate 
the integral and the quantity pL in (6). The 
result is, after a few steps, 


sad 
S fedv 


‘fae 


since f—1 in the limit. 


Lim (A —B)=Lim f{e- 








.’ This assumption may break down in regions of absorp- 
tion by water vapor or carbon dioxide in the atmosphere 
unless these materials are excluded from the light path. 


If this expression vanishes, the desired theorem 
is proved. But it will vanish if the absorption 
coefficient a is constant over a slit width, a much 
less stringent condition than that its exponential 
be considered constant. Moreover, under certain 
conditions this expression will vanish regardless 
of the variations of a. Thus if g(v, v’) =g(v—»’) 
only, i.e., if the resolving power is constant over 
the width of the band, and if g(v—»v’) =g(v’—»), 
then 


for de= feo—r)do-¥1)=6, (9) 


which is independent of v’. Also 


J fesdear’= [af gdv'dr=G f adr, (10) 


so the expression in Eq. (8) will reduce to* 


feods’— fadr=0 (11) 


Consequently, if the incident intensity J» does 
not vary rapidly over a slit width and if either 
the resolving power is high compared with the 
variations in @ or the resolving power does not 
change much over the band, 

Lim B=A. (12) 


pL-0 


Therefore by measuring B at a series of values 
of pL (either by changing p or L) and extrapo- 
lating to zero pL, the true integrated absorption 
A can be found. 


III. PRESSURE BROADENING 


In practice the fluctuations of f=exp (—apL) 
will be very extreme on passing through each 
line of the rotational fine structure. Conse- 
quently, where such fine structure exists, it might 
be necessary to carry the measurements to im- 
practically low values of pL before a good limit 
could be estimated. To remedy this, the rota- 
tional fine structure can be eliminated by pres- 
sure broadening with the use of a sufficiently 
high pressure of a non-absorbing foreign gas, as 
was done by Bartholomé‘ for HCI. 

3 A similar argument has been used by D. M. Dennison, 
Phys. Rev. 31, 503 (1928), for a slightly different purpose. 


4E. Bartholomé, Zeits. f. physik. Chemie B23, 131 
(1933). 
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With spectrographs of high resolving power, 
the use of a foreign gas to eliminate the rotational 
fine structure may be sufficient to make B 
essentially equal to A without extrapolation 
but for poorer spectrographs the variation of 
f=exp (—apL) over the band envelope may 
still be sufficiently great to cause very large 
differences between A and B, as was pointed out 
by Kemble.’ The extrapolation procedure then 
takes care of this. 


IV. ABSORPTION AREAS 


Bourgin,® in studying the intensity of indi- 
vidual HCI rotation-vibration lines, extrapolated 
the integral (“‘absorption area’’) 


pLc= { ((T.—1)/ToMe’ (13) 


instead of B. By expansion of the logarithm in B, 
one gets 


B=(1 ipl) f(s Bx?+-+-)dy’, (14) 


where 

x=(T)—T)/To, (15) 
so that in the limit B and C are also equal. 
However, B is always larger than C and various 
arguments indicate that B is normally less than 
A so that B extrapolates to A more conveniently 
and rapidly than C. As shown by Dennison,* C 
is independent of the resolving power whereas B 
is not (at finite values of pL). Therefore the use 
of a better spectrograph assists the extrapolation 
when B is used but does not with C. 


V. EFFECT OF STRAY LIGHT 


In most spectrographs light of a frequency 
quite foreign to that being studied enters the 
detecting instrument, which if not taken into 
consideration will falsify the results. If the 
amount of stray is known, it can be corrected for, 
either by subtracting the stray intensity S from 
each observation 7” or by correcting B after 
integration. If the stray is greater than about 
10 percent of the total intensity Jo, it is more 
accurate to correct each observation. 


5 FE. C. Kemble, J. Chem. Phys. 3, 316 (1935). 
§ PD. G. Bourgin, Phys. Rev. 29, 794 (1927). 


AND A. J. 


WELLS 


Let T)=T) —S and T=T’—S so that the 
uncorrected integral B’ is 


B’ =(1/pL) f In [To'/T’)dv’ 


=(1 ie) f In [(To+S)/(T+S) dv’. 
If Sis independent of pL, then by differentiation 
one gets 


Lim B’ = —Lim fosr+sy cer ‘d( pL) |dy' 


dT /d(pL) 
= — Lim f dv 
L1+(S/To)]To 
(1/[1+(S/T>) ]) Lim B, 





I 


if 1+(S/7 >) is constant over the band. If the 
stray is not too large, it should therefore not 
interfere with the accuracy of the extrapolation 
method, since it can be corrected for. It should 
be noted, however, that the stray may also be 
absorbed by the material under study, thus 
introducing a further source of error. 


VI. SUMMARY OF METHOD 


The above considerations show that it should 
be possible under proper conditions to obtain 
values of the intensity of an infra-red vibrational 
absorption band (integrated over the various 
rotational lines) by measuring the integrated 
apparent absorption coefficient for unit pressure 
B (see Eq. (5)) at a series of partial pressures of 
the absorbing gas and extrapolating this quantity 
to zero partial pressure. The rotational fine 
structure should be eliminated by means of a 
sufficiently high pressure of some suitable non- 
absorbing gas. The extrapolated value of 8 
should then equal the true integrated absorption 
coefficient at unit pressure, A. 

Care needs to be taken that the presence of 
atmospheric absorption bands does not cause the 
intensity of the incident beam to vary rapidly 
in the frequency band admitted by the slit and 
also that the resolving power of the instrument 
does not vary excessively from one end of a band 
to the other. Correction needs also to be made 
for stray light. 
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An investigation has been made of the reactions of 
|-butene and 2-butene with mercury (*P;) atoms at 30°C 
in a static system. The main products of the 1-butene 
reaction are 2-butene and liquid polymer. A large number 
of other products are also formed especially at low pres- 
sures. The rate of isomerization increases rapidly with 
pressure in the region of complete quenching, while the 
rate of polymerization falls off rapidly. The over-all rate 
of consumption of 1-butene increases very slightly with 
increasing pressure. The over-all quantum yield of 1-butene 
consumption is 0.32 at 20-mm pressure. The principal 
steps in the mechanism proposed for the reaction are: 

1 —C,Hs+Hg(?P1)—1 —C,Hs*+Hg('So), 
1 —C,H,*+ 1 —C,Hs—2 —C,Hs+1 —C,Hsg, 
1 —C,Hs*—C,H;+H, 


followed by polymerization reactions of the type 
1 —C,Hs+C,H7—>CsH i, etc. 


It is therefore possible for the first time, to obtain a direct 
measurement of the rate of collisional deactivation merely 
by determining the rate of formation of 2-butene. 2-butene 
is much less reactive than 1-butene. The over-all quantum 
yield of 2-butene consumption has a maximum value of 
0.11 at an initial pressure of 1-2 mm. Above this pressure 
the rate of 2-butene consumption falls off rapidly. The 
main product of the reaction is a liquid polymer. Smaller 
amounts of hydrogen and methane are also formed. An 
activated molecule mechanism is proposed which accounts 
for the essential aspects of the reaction. 





INTRODUCTION 


N earlier papers from this laboratory it has been 
shown that ethylene,' butadiene,” propylene,* 
and isoprene’ all react with mercury (*P;) atoms, 
at room temperature, by an activated molecule 
mechanism. Isobutene,‘ on the other hand seems 
to decompose initially into C,H; radicals and H 
atoms. In other words, with isobutene, the acti- 
vated molecule formed by collision of the hydro- 
carbon molecule with the mercury (*P;) atom 
seems to have a shorter life than the interval 
between successive collisions. 

The present investigation on 1-butene and 
2-butene was undertaken in order to investigate 
further the relation between structure and 
mechanism of reaction in the unsaturated hydro- 
carbons. 


EXPERIMENTAL 


A general description of the apparatus em- 
ployed, and the experimental technique has 
already been given in a previous paper.” 


*Contribution No. 1414 from the National Research 
Laboratories, Ottawa, Canada. 

'D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 9, 
829 (1941). 

*H. E. Gunning and E. W. R. Steacie, J. Chem. Phys. 
12, 484 (1944). 

*H. E. Gunning and E. W. R. Steacie, J. Chem. Phys. 
14, 57 (1946). 

*H. E, Gunning and E. W. R. Steacie, J. Chem. Phys. 
14, 534 (1946). 


The cylindrical quartz cell had a volume of 
240 cc, and was immersed in a water thermostat, 
maintained at 30.00+0.01°C. The total volume 
of the system was 1050 cc. 

The analyses of the condensible products of 
the reaction reported in this paper were made 
on the mass spectrometer through arrangement 
with the National Bureau of Standards. The 
sampling technique employed, as well as the 
method of analysis of the products, non-con- 
densible in liquid air, have been previously 
described. 

The 1-butene and 2-butene were the C.P. 
grade manufactured by the Matheson Company, 
East Rutherford, New Jersey. After several 
trap-to-trap distillations, the middle fractions 
were collected and stored in one-liter flasks 
provided with Warrick-Fugassi valves.® 


TABLE I. Composition of original reactants. 


—————”—”~”s«dbutene”———iti‘éOCCé ten 


Component mole percent 


C3H, —_ 0.3 
Propylene — 0.4 
1-Butene 98.91 — 
2-Butene 0.87 99.0 
Butanes 0.13 0.2 
2-methyl-2-butene = 0.03 
Methyl pentenes 0.02 
Dimethyl pentanes 0.05 
Carbon dioxide -- 


5 E. Warrick and P. Fugassi, Ind. Eng. Chem. Anal. Ed. 
15, 13 (1943). 
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TABLE II. 1-Butene. 








B 
d(He) 
—— X10 

dt 


A 
d(CHa) 


—AP 

— X10? 
Run Po At 
No. (mm) mm/min. 


> 
~~ 
& 


moles/min. 


0.0186 
0.0290 
0.0427 
0.0504 
0.0599 
0.0585 
0.0536 
0.0422 
0.0426 
0.0375 
0.0416 
0.0377 
0.0429 
0.0273 
0.0227 
0.0208 


moles/min. 


0.0280 
0.0540 
0.0719 
0.0686 
0.0762 
0.0704 
0.0712 
0.0619 
0.0550 
0.0520 
0.0579 
0.0522 
0.0542 
0.0340 
0.0286 
0.0257 





2.12 
4,92 
6.03 
7.91 
8.66 
9.90 
15.35 
16.36 
17.49 
20.92 
22.52 
23.67 
24.74 
44.81 
62.38 
67.23 


1.45 
1.82 
1.73 
1.88 
1.83 
2.02 
1.83 
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76G 
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The compositions of the gases used, as deter- 
mined on the mass spectrometer, are given in 
Table I. 

RESULTS 


1-Butene 


1-butene shows a similar behavior to propy- 
lene,’ isoprene,’ and isobutene* in that the 
pressure decreases linearly with time even in the 
initial stages of the reaction. The polymer forms 
in colorless droplets on the window of the cell 
through which the light enters. The gas, non- 
condensible in liquid air, present at the end of 
each run’ was shown to be a mixture of hydrogen 
and methane. These data are summarized in 
Table II. 

From Table I1 we see that the rate of pressure 
decrease as well as the rates of formation of 
hydrogen and methane all reach a maximum at 
an initial pressure of 10 mm. Measurements of 
the actual amounts of methane and hydrogen 
formed, as functions of the time, were made at 
initial pressures of 15 mm and 30 mm. These 
quantity-time curves were all straight lines 
passing through the origin. 

The last column of Table II shows that the 
ratio [d(CHy,) |/dt/[d(He) |/dt is constant, within 
experimental error, at 1.39+0.12. 


THE PRODUCTS OF THE REACTION 


A detailed analysis of the products of the reac- 
tion has been made on the mass spectrometer at 
the National Bureau of Standards by Dr. A. 


AND E. 
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TABLE III. 1-butene. Mass spectrometric analysis— 
run 80G. Po=2.12 mm. 








Moles products formed per mole 
1-butene decomposed 


0.0166 
0.0250 
0.1372 
0.0230 
0.00204 
0.00353 
0.00676 
0.0366 
0.00101 
0.00605 
0.00524 
0.00232 
0.00030 
0.00434 
0.00040 


Product 





Hydrogen 

Methane 

Ethane 

Propane 

2-Butene 

n-Butane 

Pentenes 

Pentanes 

Benzene 

Methyl cyclopentane 
Hexanes 

Methyl! cyclohexanes 
Ethyl cyclohexenes 
Heptanes 

Octenes 





Keith Brewer® and these results are reported in 
Tables III, IV, V, VI, and VII. There are a 
large number of products at the lower pressures 
owing undoubtedly to radical recombination 
reactions. At the higher pressures, the main 
products are 2-butene and polymer. The polymer 
is a clear, viscous, colorless liquid, completely 
soluble in benzene. 

The data in Table VII on the composition of 
the polymer, and the rate of polymerization of 
1-butene, were obtained by a carbon-balance 
calculation. 

In Fig. 1 the over-all rate of consumption of 
1-butene, the rate of formation of 2-butene, 
and the rate of polymerization are plotted as 


TABLE IV. 1-butene. Mass spectrometric analysis 
run 81G. Po=4.92 mm. 








Moles product formed per mole 
1-butene decomposed 


0.0141 
0.0263 
0.1126 
0.0240 
0.2124 
0.0207 
0.00296 
0.0416 
0.0124 
0.00104 
0.0108 
0.0137 
0.00226 
0.00191 
0.00330 
0.00834 
0.00070 


Product 





Hydrogen 

Methane 

Ethane 

Propane 

2-Butene 

n-Butane 
Pentadienes 
Pentenes 

Pentanes 

Benzene 

Hexenes 

Hexanes 

Methyl cyclohexenes 
Ethyl cyclohexenes 
Heptenes 

Heptanes 

Octanes 








®For a description of the method of analysis, see A. 
Keith Brewer and Vernon H. Dibeler, J. Research Nat. 
Bur. Stand. 35, 125 (1945). 
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TABLE V. 1-butene. Mass spectrometric analysis— 
run 87G. Po=67.23. 








Moles product formed per mole 








Product of 1-butene decomposed 
Hydrogen 0.0091 
Methane 0.0113 
2-Butene 0.742 








functions of the initial pressure, using the data 


of Tables VI and VII. 
QUANTUM YIELD 


Four determinations were made of the quan- 
tum input to the cell from the mercury resonance 
lamp by filling the cell with ethylene at an initial 
pressure of 13 mm. LeRoy and Steacie! have 
shown that this reaction has a quantum yield of 
0.37 at this pressure. The mean quantum input 
to the cell, on this basis, was 7.07+0.1110-° 
einstein per minute. The values of the over-all 
quantum yield of 1-butene consumption are 
given for each pressure in the last column of 
Table VI. 

RESULTS 
2-Butene 


The results for the runs on 2-butene are sum- 
marized in Table VIII. 

It is apparent from Table VIII that 2-butene 
is considerably less reactive than 1-butene. The 
main product of the reaction is a liquid polymer. 
Smaller amounts of hydrogen and methane are 


TABLE VI. 1-butene. 














—d(1—CsHs) d(2—CaHa) 
Quantum 
Run Po dt yield 
No. mm, moles/min. 106 ¢(1 —C4Hs) 
80G 2.12 1.12 0.00229 0.159 
81G 4.92 2.06 0.437 0.291 
83G 8.66 1.71 0.526 0.242 
84G 15.35 2.23 0.967 0.316 
85G 22.52 2.25 1.133 0.318 
86G 44.81 2.13 1.285 0.301 
87G 67.23 2.28 1.693 0.323 








TABLE VII. 1-Butene. 








Rate of polymerization of 








Run Polymer 1-butene 

No. composition moles per min. X10 
80G n (C4H7.4) 0.92 

81G n(CyH7.2) 1.10 

83G n(C4Hzs) 1.16 

84G n(C4H7.8) 1.24 

85G n(C4H;.s) 1.12 

86G n(C4H7.9) 0.83 

87G n(C4H7.s) 0.58 
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also formed. The maximum rate of consumption 
of 2-butene seems to occur at an initial pressure 
of 1-2 mm corresponding to the maximum over- 
all quantum yield of approximately 0.11. 


DISCUSSION 
1-Butene 


Before attempting to devise a mechanism for 
this reaction, it might be well to summarize 
briefly the salient results of our investigation. 
These are: 


(a) The pressure vs. time curve is linear even 
in the initial stages of the reaction. 

(b) The rate of pressure decrease is maximum at 
10 mm. 

(c) The main products of the reaction are 
2-butene and a liquid polymer. Smaller 
amounts of a large number of other sub- 
stances are also formed, including methane 
and hydrogen. | 

(d) The rate of isomerization increases rapidly 
with pressure in the region of complete 
quenching, while the rate of polymeriza- 
tion falls off rapidly. 

(e) The over-all rate of consumption of 1-butene 
increases only very slightly with increasing 
pressure in the complete quenching region. 
The quantum yield is approximately con- 
stant in this region at 0.32. 

(f) The rates of formation of methane and 
hydrogen reach a maximum at 10-mm 
pressure. Furthermore the ratio 

[d(CHg) /dt]/(d(He) /dt } 
is constant throughout the entire pressure 
range at 1.39+0.12. These products are 
only of minor importance compared to 
isomerization and polymerization. 

(g) The polymer is a clear, colorless liquid, com- 
pletely soluble in benzene. 


The fact that the rate of isomerization in- 
creases with pressure, in the complete quenching 
region, while the rate of polymerization de- 
creases, suggests the sequence 


1—C,Hs+Hg(@P1) 1 —CuHs*+Heg('So), (1) 
1—C,Hs*->C,H;+H, (2) 
1—C,4H,*+1—C,Hs—>2—CyHs+1—CaHs, (3) 
1—C,H,*+1—CyHs—2(1 —CaHs). (4) 
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Polymerization will then be initiated by the 
radicals formed in (2) 


C,Hs+H—-CyHsg, (5) 
C,H;+1 —CyHs—CsH is, etc., (6) 
C,H,+1 — C,H s—-CsH 17, etc. (7) 


On this basis, isomerization would be vir- 
tually the only reaction occurring at high pres- 
sures. An examination of Table VI will show 
that this conclusion is consistent with our 
findings. At Po=15.35 mm, 0.967/2.23 or 0.43 
mole of 2-butene are formed per mole of 1-butene 
decomposed. While at Po=67.23, 1.693/2.28 or 
0.74 mole of 2-butene are formed per mole of 
1-butene decomposed. Consequently one would 
expect that at pressures above, say, 150 mm, 
free radical formation, and hence polymerization 
would be of negligible importance, the only reac- 
tion occurring being the isomerization of 1-butene 
to 2-butene. It is interesting to note at this point 
that LeRoy and Steacie’ found that 2-butene was 
a significant product in the cadmium (*P;) sen- 
sitized reaction of 1-butene. 

With increasing pressure, then, reaction (2) 
will eventually become completely suppressed 
and under these conditions the rate equations 
assume the asymptotic form 


d(2—C,4Hs) —d(1—C,Hs) 
stroll PO / tb 


dt dt 
Rg 
=h( —) = 0.32 
ks+k, 


7D. J. LeRoy and E. W. R. Steacie, J. Chem. Phys. 10, 
683 (1942). 








from Table VI or 


kg fki—0.32 
ks ( 0.32 ) 

Now ifall collisions between activated 1-butene 
molecules and normal 1-butene molecules lead to 
isomerization k,=0, and hence ki =0.32. This is 
not an improbable assumption in view of the 
fact that it has been shown that the initial step 
in the mercury photo-sensitized reaction of 
isobutene‘ is only about 16 percent efficient. In 
other words, if we assume that reaction (1) is 
only 32 percent efficient, the rate of formation 
of 2-butene becomes a direct measure of the rate 
of deactivation of the excited molecules. 

This conclusion is of considerable interest. 
Collisional deactivation has been postulated in 
a large number of reactions in order to explain, 
for example, a decrease in rate with increasing 
pressure. Now, for the first time, it becomes 
possible to follow deactivation directly simply 
by measuring the rate of formation of 2-butene. 

The small amounts of methane and hydrogen 
formed probably arise from the unimolecular 
decomposition of the activated butyl radical 
formed in (5). This radical will contain excess 
energy of hydrogenation‘ and may decompose to 
some extent by the reactions 


C,H s—H2+ CyH;, (8) 
C,Hy—>CHi+CsH;. (9) 


These reactions would explain the fact that 
the ratio of the rates of formation to methane to 
hydrogen is constant, since under these condi- 
tions 
se ks 


=—=1.39 from Table il. 
dt dt Rg 


At low pressures atomic cracking reactions of 
the type 





H+C,H»—-CH;+C3H;, (10) 
H +C,H s—2C2H;, (11) 


would become significant. The important prod- 
ucts listed in Tables III and [V could easily be 
accounted for by radical recombination reactions 
involving CHs, C2Hs, C3H;, C,H, and C,Ho 
radicals. 
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REACTIONS OF 1-BUTENE AND 2-BUTENE 


2-Butene 


The salient features of the 2-butene reaction 
are: 


(a) The rate of pressure decrease, as well as the 
rate of consumption of 2-butene both fall 
off with increasing pressure above 2 mm. 

(b) The maximum over-all quantum yield of 
2-butene consumption is 0.11. 

(c) The main product of the reaction is a liquid 
polymer. Smaller amounts of hydrogen 
and methane are also formed. 


The falling off in the rate of consumption of 
2-butene with increasing pressure again suggests 
an activated molecule mechanism. 


2-— C,Hs+Heg(*P;)— 
2—C.Hs*+Hg('So), (12) 


2—CyHe*+2—(C.Hs)2(2—CiHs), (13) 
2—C.H.*—>C.H7+H, (14) 
2—C.Hs+H—CiHy, (15) 
CyHo+2—CiHs CoH, etc., (16) 

C,H; +2-C,Hs CH ss, ete. (17) 


A comparison of Tables VII and VIII shows 
that the maximum in the polymerization rate 
curves occurs at an initial pressure of 1-2 mm 
for 2-butene, and at about 10 mm for 1-butene. 
Now since the quenching cross sections for the 
two molecules should not be significantly dif- 
ferent, we are lead to conclude that the activated 
state has a longer life in the case of 2-butene than 
in 1-butene, i.e., (1/ki4) >(1/k2). This assump- 
tion would also offer a plausible explanation for 
the fact that the rate of polymerization of 
2-butene is only about one-half that of 1-butene. 

The small amounts of hydrogen formed sug- 
gest the step 


H+2—C,Hs—-H24+ CyH;, (18) 


which would involve an activated C4Hg radical 
as an intermediate. 

The methane is present in significant amounts 
only at low pressures. The reaction 


CyHy—>CH4+C3Hs (19) 


would require that the rate of methane formation 
fall off even more rapidly with pressure than the 


TABLE VIII. 2-butene. 








—AP —d(2—CsHs) d(H2) d(CHa) 
ance GE deena 

Po At dt 
mm mm/min. 


2.02 0.93 
4.63 0.82 
8.32 0.65 

21.52 , 

32.36 

42.14 

49.02 

61.14 

70.64 


¢(2 —CaHs) 


0.041 0.047 0.11 
0.047 0.024 

0.038 0.009 

0.027 0.003 

0.025 0.001 

0.022 0.001 

0.022 0.0008 

0.020 0.0007 

0.020 0.0002 


dt 
moles/min. 106 
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rate of hydrogen formation, in agreement with 
our results. 
It is interesting to note that the reaction 


2 —C,Hs*+2 —C,H;—1 —C,Hs+2 —(C,Hs 


does not seem to occur, since no 1-butene was 
found in the products. Thermodynamically, of 
course, 2-butene is more stable than 1-butene, 
although their heats of formation differ only by 
1-2 kcal.’ 

CONCLUSIONS 


It is apparent that, with the striking exception 
of isobutene, the activated molecule mechanism 
seems to be common to all the lower olefines and 
alkadienes. 

The collisional isomerization reaction of 1-bu- 
tene to 2-butene is of considerable interest since 
it gives us a direct_means of determining the 
extent of collisional deactivation. Virtually every 
activated molecule formed seems to undergo a 
deactivating collision above an initial pressure of 
150 mm. We are therefore lead to the important 
conclusion that, with the exception of isobutene, 
free radicals play an important part in the 
mercury photo-sensitized reactions of the lower 
olefines and alkadienes only at relatively low 
pressures. 
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Force constants have been calculated for the out-of-plane bending motions of the vinyl 
halides and for cis- and trans- 1,2-dichloroethylene. Published vibrational assignments have 
been followed for the most part, but several changes were found necessary. A generally con- 
sistent set of force constants was obtained in which the values are also similar to those for the 


methyl halides. 





HILE bending force constants in general 
have received less attention than stretch- 
ing force constants, they are important in estab- 
lishing the assignments of the lower fundamental 
vibration frequencies. In thermodynamic calcu- 
lations the lower frequencies take on primary 
significance. In this investigation, an attempt 
has been made to evaluate the force constants 
for the non-planar vibrations in certain halo- 
genated ethylenes, and to find some correlation 
among them for this series of compounds. 
Complete assignments of fundamental vibra- 
tion frequencies for the vinyl halides have been 
made by Torkington and Thompson!? on the 
basis of their own infra-red spectra. These 
assignments appear to be straightforward and 
unambiguous except in the case of vinyl fluoride. 
Here some anomalies occur that are not yet 
satisfactorily explained. The fluoride, further- 


TABLE I. Non-planar frequencies and force constants of 
vinyl halides. 





CeHsI 
2.03 


CeoH;Cl CeHsBr 
1.69 1.85 





CeH3F 
1.36 





C—X distance, A 
Frequencies (cm~!) 
1 924 (860) 
we 860 (732) 
w3 715 (715) 
Assumed force 
constants 
Fi, (from ethylene) 
Fi. (= Fo) 
(ethylene) 
Fy3 (= Fs) 
Calculated force 
constants 


946 
909 
435 


940 
902 
497 


940 
895 
622 


0.39 


0.054 
0.00 


0.37 0.39 


0.045 
0.00 


0.39 


0.054 
0.00 


0.054 
0.00 


0.32 
0.09 
0.60 


0.48 
—0.21 
0.90 


0.43 
—0.10 
0.82 








1P. Torkington and H. W. Thompson, J. Chem. Soc. 
303 (1944). 

? P. Torkington and H. W. Thompson, Trans. Faraday 
Soc. 41, 236 (1945). 


more, is the only one of these substances for 
which Raman data are lacking. Wu’ and Herz- 
berg’ have summarized the data on the di- 
chloroethylenes. 

The calculations were carried out by the 
method of Wilson® for which the following co- 
ordinates were chosen : 
6,;=angle of wag of CHe group (angle between extension 

of C=C bond and HCH plane). 


62=angle of wag of CHX group. 
6; =angle of twist of groups at opposite ends of C=C bond. 


The elements of the reciprocal kinetic energy 
matrix then take the following form for the 
vinyl halides: 


Gy =2e@unt+uc(4e+4ep+2p?) ; 
Gi2=Go1= —ucp(3e+n+ 2p) ; 
Gi3=G31=0; . 
Goo = Cua tnuxtuc(e+n?+2p?+ 2en+ 2ep+ 2np); 
Ge3=Gs2= (eux i nex) /V3 ’ 
G33=@unt+3n'ex, 
where the u’s refer to reciprocal masses, and e, 
p, and 7 are the reciprocals of the C—H, C=C, 
and C—X bond distances, respectively. The 
C—C and C—H bond lengths (1.353 and 1.0714) 
are taken to be those in ethylene, as determined 
by Gallaway and Barker.* These authors report 
a bond angle of 119° 55’ for H—C—H in ethyl- 
ene, whereas electron diffraction measurements 
on vinyl chloride by Brockway, Pauling, and 
Beach’ indicate a C—C—Cl angle of 122+2°. 
3T. Y. Wu, Phys. Rev. 46, 465 (1934). ; 
4G. Herzberg, Infra-Red and Raman Spectra (D. Van 
Nostrand Company, New York, 1945), p. 330. 

6 E. B. Wilson, Jr., J. Chem. Phys. 7, 1047 (1939) and 
9, 76 (1941). See also A. G. Meister and F. F. Cleveland, 
Am. J. Phys. 14, 13 (1946). 

6 W. S. Gallaway and E. F. Barker, J. Chem. Phys. 10, 
88 (1942). 


7L. O. Brockway, L. Pauling, and J. Y. Beach, J. Am. 
Chem. Soc. 57, 2693 (1935). 
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For simplicity, 120° has been used for all bond 
angles in these calculations. The experimental 
value of 1.69A given by Brockway, Pauling, and 
Beach,’ for the C—Cl distance is accepted. This 
represents a shortening of the bond compared 
with 1.76A in CCly, and a corresponding relation- 
ship is assumed to exist in the other halides, 
i.e., the =C—X length is decreased with respect 
to that of —C—X. This is corroborated in the 
case of vinyl iodide, in which the C—I length is 
reported by Huggill, Coop, and Sutton’ as 
2.03A, while in Cl, it is 2.12A’ (Finbak and 
Hassel®). The C—X distances are included in 
Table I. 

In these calculations the force constant for 
CHe wagging, F1,, and the interaction constant, 
Fy, were assigned the values obtained from 
ethylene.” This proved satisfactory in all cases 
except vinyl fluoride, where the noticeably lower 
value of w; (compared with ethylene and the 
other halides) made it necessary to lower these 
constants slightly. All constants are expressed 
in “atomic weight-angstrom” units, which for 
angle bending constants is (ergs/radian*) X 1.698 
x 10+". The frequencies and constants are listed 
in Table I. In parentheses are given Torkington 
and Thompson’s? assignment for vinyl fluoride. 

In an attempt to confirm the F22 value for 
C.H;Cl, some calculations were made of the 
out-of-plane force constants in the cis- and trans- 
forms of 1,2-dichloroethylene. The appropriate 
equations are readily derived by the same meth- 
ods. In spite of some uncertainties in the assign- 
ments of these frequencies,** it is believed that 
the results obtained provide some measure of 
verification for the vinyl chloride wagging con- 
stant. The choice and classification of the fre- 
quencies is given by Table II, which also shows 
the force constants obtained from them. 

These molecules have higher symmetry than 
the vinyl compounds, so that the non-planar 
motions may be further subdivided. The Bz 
“unsymmetrical” wagging frequency in cis-1,2- 
dichloroethylene, and the B, “symmetrical” 
wagginig frequency in trans-1,2-dichloroethylene 
are both (by symmetry) free from interaction 





_*J. A. C. Huggill, I. E. Coop, and L. E. Sutton, Trans. 
Faraday Soc. 34, 1518 (1938). 

* C. Finbak and O. Hassel, Zeits. f. physik. Chemie 36B, 
301 (1937). 
‘© J. E. Kilpatrick and K. S. Pitzer, to be published. 


BENDING FORCE CONSTANTS 


TABLE II. 1,2-dichloroethylenes. 

















Frequency Force 
Vibration (cm!) constant 
Ccis- 
Bz (unsym. wagging) 694 0.45 
As (sym. wagging) 876 0.36* 
As (torsion) 406 0.78 
(A2 interaction constant) — 0.08 
trans- 
B, (sym. wagging) 765 0.36 
A, (unsym. wagging) 820 0.45t 
A, (torsion) 270 0.87 
(A, interaction constant) —0.17 
deutero-trans- 
B, (sym. wagging) 657 0.36 








* Assumed from By (sym.) in trans-. 
+ Assumed from Be (unsym.) in cis-. 


with the other non-planar vibrations. This leads 
to a rather simple and direct evaluation of these 
force constants, which in the case of the trans- 
molecule can also be checked with the deuterated 
compound. Then if the “‘symmetrical’’ wagging 
constant in trans- is accepted for the correspond- 
ing motion in cis-, and the “unsymmetrical’”’ 
wagging constant in cis- is accepted similarly for 
trans-, the torsional and interaction constants 
for the Az and.A, vibrations may be calculated 
and a set of values consistent with the most 
probable frequencies obtained. The greatest 
uncertainties in the assigned frequencies are in 
the A, vibrations of trans-1,2-dichloroethylene. 
While Wu’ and Herzberg‘ both suggest 620 cm7! 
for the ‘unsymmetrical’ wagging vibration of 
this group, it gives too low a force constant. 
Wood and Stevenson" have used 820 cm™ for 
this frequency in connection with studies of the 
cis-trans isomerization, and the present calcula- 
tions seem to confirm the latter value. The A, 
torsional frequency is too low to be observed. 
The 620 cm band, left unassigned by the choice 
of 820 cm= for the A, wagging fundamental, 
may now be assumed to be a combination of 
two low frequencies. One of these must be the 
350 cm-! C—CI planar bending, leaving a differ- 
ence of 270 cm~. This is taken as the twisting 
fundamental. 

The symmetrical and unsymmetrical wagging 
constants for the dichloroethylenes should be 
approximately Fo22—Fy. and F22+Fi2, respec- 
tively, from vinyl chloride. While the agreement 
with these combinations from vinyl chloride is 


11R. E. Wood and D. P. Stevenson, J. Am. Chem. Soc. 
63, 1650 (1941). 
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not very close, it is about what one might expect. 
Likewise, the torsion constants of the cis- and 
trans-molecules differ somewhat. However, the 
pattern of constants shows such general con- 
sistency as to lend considerable support to their 
validity. 

To examine the trend of the force constants in 
passing through the series of vinyl halides, it is 
informative to see whether or not a similar 
trend exists in some comparable series. The 
methyl halides afford such a comparison, and the 
force constants for bending of the H—C—X 
angles in these molecules are given by Crawford 
and Brinkley.” Their values, expressed in terms 


2B. L. Crawford, Jr. and S. R. Brinkley, Jr., J. Chem. 
Phys. 9, 69 (1941). 


ROBINSON 


of our units are as follows: 


CH;3F CH;Cl CH;Br CHI 
k 0.41 0.31 0.28 0.24 


The relative values of k may be properly com- 
pared with those of F22 in the vinyl compounds 
(see Table I). In both series, it is observed that 
progression from the fluoride to successively 
heavier halides leads to successively lower force 
constants. 

We are aware that our assignment for viny] 
fluoride will make it difficult to explain the 
published spectrum completely; however, any 
other assignment gives unreasonable force con- 
stants. Further investigation is needed here. 
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The relative partial molal entropy of water, S:—S,°, in aqueous solutions of seven amino 
acids is calculated from existing data. For all solutions, 8:—S,°=km? approximately; k is posi- 
tive for glycine solutions and negative for solutions of the two alanines and four aminobutyric 
acids. The S,—8S,° values are determined largely by the size and structure of the hydrocarbon 
part of the amino acid and only to a minor degree by its dipole moment. Positive S,—S,° values 
are interpreted as indicating a structure breaking effect of glycine on water; negative S,—S,° 
values indicate a structure strengthening effect which increases with size of the hydrocarbon 
residue. The data are consistent with the recent discussion of Frank. Heat capacity and 
viscosity values for some of these solutions qualitatively show expected behavior. 





TUDIES of the heats of dilution of aqueous 
solutions of amino acids'~* have shown that 
the dipole moment is not the main factor deter- 
mining the thermodynamic properties of these 
solutions. The quantity ®z2 (apparent relative 
molal heat content) is negative for glycine, has 
positive values for a-alanine and shows a 
further increment for a-amino-n-butyric acid, 
although the dielectric increments for these 
three substances are approximately equal ;‘ in 
1W. E. Wallace, W. F. Offutt, and A. L. Robinson, 
J. Am. Chem. Soc. 65, 347 (1943). 
2 Data for a- and f-alanine, L. S. Mason, H. Benesi, 
and A. L. Robinson, J. Am. Chem. Soc., to be published. 
3 Data for four aminobutyric acids, L. S. Mason and A. 
L. Robinson, J. Am. Chem. Soc., submitted for publication. 
4E. J. Cohn and J. T. Edsall, Proteins, Amino Acids, 


and Peptides (Reinhold Publishing Corporation, New York, 
1943), p. 146. 


contrast, a-amino-n-butyric and y-aminobutyric 
acids show nearly identical @z, values up to 
1.2 m although the dielectric increment of the 
y-acid is more than twice as large as that of the 
a-acid. The hydrocarbon residue is evidently an 
important factor. Its influence may be exerted 
through hydrocarbon-hydrocarbon, hydrocarbon- 
dipole, and hydrocarbon-solvent interactions. 
An examination of the differential entropies of 
dilution (S,—S,°) of aqueous electrolyte solu- 
tions’ permitted an attempt to correlate the 
sizes, shapes, and hydrogen bonding powers of 
ions with changes in the structure of water in 
5H. S. Frank and A. L. Robinson, J. Chem. Phys. 8, 
933 (1940). A recent paper by G. W. Stewart, J. Chem. 


Phys. 11, 72 (1943) is a particular interest in this connec- 
tion, 
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ENTROPY OF AQUEOUS AMINO ACID SOLUTIONS 589 


TABLE I. Differential entropy of dilution ($,—8S,°), in calorie-degree!-mole™, in aqueous solutions at 25°. 











\ m 0.2 0.4 0.6 0.8 1.0 1,2 

amino acid \_ 

glycine 0.0312 0.0338 0.0382 0.0015 0.0023 

a-alanine — 0.0.88 — 0.0335 —0.0;79 —0.0014 —0.0022 —0.0032 
6-alanine —0.0311 — 0.0336 —0.0;66 —0.0013 — 0.0020 — 0.0026 
a-amino-n-butyric —0.03;23 —0.0;91 —0.0021 — 0.0037 —0.0058 — 0.0084 
a-amino-i-butyric —0.0;30 — 0.0012 — 0.0027 —0.0048 —0.0076 —0.0107 
8-amino-n-butyric — 0.0325 —0.0;95 —0.0021 — 0.0037 — 0.0056 — 0.0080 
y-aminobutyric — 0.0013 — 0.0029 — 0.0050 — 0.0107 


— 0.0335 





— 0.0076 











the solutions. Recently Frank® has discussed 
entropies of vaporization of solutes from solu- 
tions and has interpreted the data to indicate 
that frozen patches or microscopic icebergs form 
around non-polar solutes and the non-polar parts 
of polar solute molecules in aqueous solutions. 
Data are now available for calculating 8,—8,° 
values for aqueous solutions of seven amino acids. 
These values lead to some possible conclusions 
regarding the effect of the hydrocarbon residues 
on the structure of the solvent, which are consist- 
ent with the interpretations of Frank. 

The §,—S,° values of Table I were calculated 
from 


S8:—8)°=(L,—RT In (a,/N))/T. (1) 


L, is the relative partial molal heat content of 
the water (H,—H),°) and a, is its activity, 
referred to pure water as the standard state. 
N,, the mole fraction of water in the solution, 
takes account of the entropy of mixing; the 
values of Table I give the non-ideal part of the 
differential entropy of dilution. For consistent 
comparison all activity values were derived from 
the vapor pressure measurements of Smith and 
Smith.7-* For the same reason, all LZ, values were 
derived from heat of dilution measurements 
made in this laboratory. Figure 1 shows the 
plots of the values of Table I against the molality. 

There is no satisfactory theory to explain the 
thermodynamic behavior of these dipolar ions in 
aqueous solution.'°- Purely electrostatic inter- 


°H.S. Frank, J. Chem. Phys. 13, 507 (1945). 

7E. R. B. Smith and P. K. Smith, J. Biol. Chem. 117, 
209 (1937). 

®E. R. B. Smith and P. K. Smith, J. Biol. Chem. 121, 
607 (1937). 

*E. R. B. Smith and P. K. Smith, J. Biol. Chem. 132, 
47 (1940). 

0G, Scatchard and J. G. Kirkwood, Physik. Zeits. 33, 
297 (1932). 

4 J. G. Kirkwood, Chem. Rev. 19, 275 (1936). 

= R. M. Fuoss, J. Am. Chem. Soc. 56, 1027 (1934). 

SR. M. Fuoss, J. Am. Chem. Soc. 58, 982 (1936). 





actions, calculated for various simplified models, 
do not suffice. From a consideration of activity 
data Cohn‘ has concluded that the nature and 
number of polar and non-polar groups, and their 
positions in the molecule, must be considered in 
order to obtain a better explanation of the facts. 

S,—S,° is the entropy change for the transfer 
of a mole of water from an infinitely dilute 
solution to a solution of finite concentration. 
Since entropy is related to structural order, the 
S,—S,° values should reflect the effect of the 
dipolar ions on the quasi-solid structure of water, 
positive values indicating structure weakening 
and negative values pointing to structure 
strengthening effects. For the seven amino acids 
studied the relation S,—S,°=km? holds fairly 
well, but only for glycine does the constant have 
a positive value.* The net effect of glycine is 
thus similar to that of sodium and potassium 
halides and nitrates at 1m while the behavior 
of the other six amino acids parallels the effect 
of ions such as H;O0+, OH-, NH,*, Lit, and F-, 
whose structure strengthening effect has been 
attributed to their small size or their hydrogen 
bridging capacity.° 

The structure strengthening effect of the 
amino acids obviously cannot be attributed to 
these same factors. The order shown by glycine, 
a-alanine, and a-amino-n-butyric acid, and also 
by 6-alanine and 8-amino-n-butyric acid, is the 
reverse of the order expected from constant 
dipole moment and increasing size. Another 
factor, the non-polar hydrocarbon groups, seems 
to play the dominant role, giving rise to a 
structure strengthening interaction with the 
water. The values of Table I are derived from 
accurate data and the differences shown by the 

* The dependence on m? is a reflection of the relative 
importance of the term in LZ; in Eq. (1). There is some 


uncertain evidence that I; for glycine! may also become 
negative at very small concentrations. 
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Fic. 1. Relative partial molal entropy of water, 8:—S,°, 
in aqueous amino acid solutions: 1, glycine; 2, 8-alanine; 
3, a-alanine; 4, 8-N-butyric acid; 5, a-N-butyric acid; 
6, a-iso-butyric acid; 7, y-butyric acid. 


members of an isomeric group are undoubtedly 
real, but it is interesting that AS, at 1m is 
— 0.0044 in going from glycine to the alanines 
(average) and almost the same, —0.0046, in 
passing from the alanines to the aminobutyric 
acids (average). Small differences are noted 
between a- and 6-n-acids; hydrocarbon branch- 
ing and the y position for the NH;* group give 
the most negative values of §,—8,°, but further 
study will be required to clarify these effects. 
The structure strengthening effect on water 
shown by the non-polar parts of these dipolar 
ions is consistent with the recent discussion by 
Frank® of the abnormally high entropies of 
vaporization of rare gases and hydrocarbon gases 
from water. Frank suggests that when a rare gas 
atom or a non-polar molecule or a polar molecule 
with non-polar parts dissolves in water at room 
temperature it modifies the water structure in 
the direction of greater ‘‘crystallinity.’’ This 
effect increases with the size of the solute mole- 
cule. Our own study leads to similar conclusions 
and simply suggests that this ‘‘freezing’’ of the 


water (Frank’s calculations are for infinite dilu- 
tion) increases with concentration in the range 
of this study. For glycine the structure breaking 
interaction, always present, is not compensated 
by the “‘iceberg”’ effect because of the small non- 
polar part. 

The high apparent heat capacities of amino 
acids and fatty acids and alcohols in aqueous 
solutions have attracted comment." In a homol- 
ogous series of liquids the increment in C, per 
CHe group is only 6 to 8 cal./mole while in 
aqueous solution the increment in ®(C,) is 20 to 
30 cal./mole. Similar behavior is shown by 
aqueous solutions of the rare gases, the effect 
also increasing with the size of the rare gas atom. 
The explanation for these large heat capacities 
given by Frank® should apply equally well to 
solutions of the amino acids. With increase of 
temperature the melting of the quasi-crystalline 
patches of water, which the entropy values 
indicate increase in size with increasing size of 
the hydrocarbon residue, will add an increasingly 
positive term to the heat capacities. 

We may also call attention to the viscosities 
of aqueous solutions of glycine and a-alanine." 
Previous discussions of the effect of solutes'® ® on 
the structure of water have pointed to the rela- 
tion between the viscosities of solutions and the 
structure weakening or strengthening influence 
of the solute. At 18° the relative viscosities of 1 m 
solutions of glycine and a-alanine are 1.153 and 
1.310, respectively. Qualitatively such a differ- 
ence would be expected in terms of the present 
discussion. More extensive data, both for en- 
tropies of dilution and viscosities, would be of 
interest for a further discussion of the effects of 
amino acids on the structure of water. 


“J. T. Edsall, J. Am. Chem. Soc. 57, 1506 (1935). 
See also reference 4, p. 168. 

15 G, Hedestrand, Zeits. f. anorg. Chemie 124, 153 (1922). 

16 J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 
(1933). 
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This paper deals with the separation of the contributing 
terms of the energy and entropy of activation for viscous 
flow and the identification of some of the molecular con- 
stants entering into these terms, as a step toward the 
prediction of the viscosity of liquids from first principles. 
The energy of activation for hole formation AH¢* is found 
to be determined mainly by the magnitude of the dis- 
persion energy and by the extent of displacement of the 
hole forming molecules from their equilibrium position. The 
latter, in the form of the viscosity-pressure coefficient A Vt, 
is, at present, not predictable from any molecular con- 
stant. The energy of activation for motion into the hole 
AH? was found for many substances to be numerically 
equal to the excess energy of vaporization AE yap’ and is 
thus an additional measure of restricted external rotation. 
The cases in which AH} >AE,,,/ offer strong evidence for 
the requirement of deformation of the flowing molecule 
against internal potential energy barriers. The rotational, 
translational, and cooperative terms which contribute to 
the entropy of activation have been separated, but can so 


far be determined only from viscosity data. Numerical 
examples show that the stereometric arrangement of 
molecular structure determines the magnitude of AVt and 
AH?! and thereby the viscosity of liquids to a far greater 
extent than chemical composition (except for OH groups). 
The existence of aliphatic fatty acids as double molecules 
in the liquid state over a wide range of temperatures is 
confirmed viscosimetrically, while the viscosity data of 
aliphatic alcohols suggest the presence of distinguishable 
multiple molecules only at very low temperatures but the 
existence of a continuous OH-bond network structure at 
ordinary temperatures. One of the important consequences 
of the rate process theory of flow is the recognition that the 
viscosity of a liquid is determined by the (very small) 
concentration of molecules in relatively shallow potential 
energy walls from which the activated molecules are 
preferentially recruited. Viscosity is therefore not a bulk 
property in the commonly accepted sense and depends only 
to a minor extent on the structure of a liquid, i.e., on its 
state of order. 





HE object of this investigation was to 

separate the terms contributing to the 
energy and entropy of activation for viscous flow 
occurring in Eyring’s' equation 

ih 


= exp (—AS+/R) exp (AH:/RT) 
AoAsA? 


Nh 
7? exp (—AS}/R) exp (AHt/RT) 


(where \j,23=intermolecular distances in the 
three space coordinates, \=length of path over 
potential energy barrier, #=Planck constant, 
V=molar volume (used where A=),), ASt 
=entropy of activation, and AH+=heat of 
activation), in order to facilitate the prediction 
of the absolute value of the viscosity of liquids 
from first principles. 


1. THE HEAT OF ACTIVATION AH; 
AH is also defined by 


a In (Vn) 
an=R(———) : (1) 
_ 01/T Jp 


* Presented in part at the annual meeting of the Society 
of Rheology, New York, October 26, 1945. 

'Glasstone, Laidler, and Eyring, Theory of Rate Proc- 
esses (McGraw-Hill Book Company, New York, 1941). 
R. E. Powell, Ph.D. Thesis, Princeton (1943). 


It has been suggested by Eyring! and his 
co-workers that AH}; consists of two elements: 
the energy required to form a hole in the liquid 
structure and the energy required to move into 
this hole. The latter is 


dln 7 
ai=R( ) : (2) 
01/T/y 


so that the energy of hole formation is 


AH?" = AH;— AH; 


0 In (Vn) dln 7 
“ie a 
01/T P 01/T/y 


So far no inquiry has been made into the nature 
of either AH" or AHy}/. It is here proposed to 
derive the two parameters as follows: From the 
model underlying the rate process theory of flow 
we can deduce that AH¢" is the change in poten- 
tial energy between neighboring molecules ac- 
companying the displacement from their equi- 
librium position when forming a new hole. The 
displacement can be computed from the pressure 
coefficient of viscosity, as 


OAFt 
Op T 
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b=6+2.0(3 cos? @B—1)+ {[(sin 6)? cos 
—2(cos 6)? ?—2 cos? 6} 


for various values of @ and ¢. 
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Fic. 1. London’s orientation function. 


where AF}=free energy of activation for viscous 
flow and AV+=local volume change (increase) 
accompanying the process of activation. The 
linear displacement of the molecules forming a 
hole is, therefore, to a good approximation 
(1/N')[(V+AV;)!— V#]. Using Lennard-Jones” 
equation for the potential energy between a 
central molecule and its nearest neighbors we 
obtain for AH" the expression 


mets) -G)! 
~L(5)-G)]} © 


where A* = Ze*, the minimum value of the inter- 
action energy ¢ at the equilibrium distance r*, 
Z=number of nearest neighbors, V*=r**/v2, 
Vit=V+AV3, A*~1.2AEyap in many cases. Inas- 
much as (5) is based on London’s inverse sixth 
power law for molecular interaction* one can 
readily show that the magnitude of « is deter- 
mined by the factor 2Joa? in the case of 
non-polar compounds and contains the addi- 
tional term (2/3)(u4/kT) in the case of di- 
polemolecules, where Jo=first ionization po- 
tential (also derivable from ultraviolet spectra‘) 

* Lennard-Jones and Devonshire, Proc. Roy. Soc. A163, 
53 (1937); A165, 1 (1930). 

3 F, London, Trans. Faraday Soc. 33, 8 (1937). 


*W. C. Price, et al. Nature 148, 372 (1942); Proc. Roy. 
Soc. A174, 207 (1940). 


and a=polarizability. One could argue that 
Eq. (5) is valid only for spherical molecules 
or at least only for those of spherical fields 
of force and that a more complicated ex- 
pression is required to describe the potential 
energy field of linear polyatomic molecules. 
London’s recent equation for the dispersion 
energy in systems composed of linear molecules® 
takes into consideration the difference in polariz- 
ability vertical to and in line with the bond axes. 
The potential energy between neighboring mole- 
cules becomes thus a function of their angle of 
approach. Figure 1 is a plot of the value of the 
numerical coefficient before [oa? as a function of 
the angles of approach. It is seen that the poten- 
tial energy is at a minimum for parallel orien- 
tation of the bond axes, in agreement with the 
known x-ray diffraction patterns of m-paraffins 
and their derivatives. In calculating the poten- 
tial energy along these curves, it becomes 
apparent that the difference in potential energy 
between the top and the bottom is of such mag- 
nitude that there should still be a finite concen- 
tration of molecules which are randomly oriented 
with respect to each other. This difference is 
often of the same order of magnitude as AF}. As 
AF‘ is defined by 


c 


Vn Ct 
AF:=RT In (*)=--z7in —, (6) 
Nh 


5 F. London, Surface Chemistry II (AAAS, 1943), p. 141. 
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where ¢c; and ¢ are the concentration of the 
activated and normal molecules, respectively, it 
appears reasonable to deduce that the activated 
molecules are (preferentially) recruited from the 
randomly oriented molecules whose numerical 
factor in front of Ioa? is so nearly equal to ? that 
the use of the simple Eq. (5) is fully justified for 
molecules of any shape, provided 


| ay a] (Baw?) ~AF. 


The close agreement of the values of AH;* cal- 
culated by Eqs. (5) and (3) indicates the validity 
of the assumptions made in the derivation of (5). 

If we remember that the energy of vaporiza- 
tion is the energy required to make a hole of 
molecular size, we should expect the ratio 
AH?"/AEyap to be of the order of (or equal to) 
AV:/V. The data on Table I show this indeed 
to be the case, what cannot be said of the ratio 
AH;+/AE,., as has been proposed by Eirich and 
Simha® and accepted by Eyring* and many sub- 
sequent workers. A number of cyclic-hydro- 
carbons and similarly constituted liquids show 
AH+/AE,,p values larger than 1,’ but their AV+/V 
and AH*/AE,., values are of the usual order of 
magnitude, i.e., around 0.2. 

As Io, a, and V* (and AHy) can be determined 
independently of viscosity, AH+" of any liquid 
could be predicted if AV were also derived from 
first principles; this has, however, not yet been 
accomplished. But it is conceivable that one 
could achieve this aim on the basis of Alfrey’s 
model of the flow process.’ The peculiar tem- 
perature dependence of A//;" (it goes through a 
maximum as does AV+t) is not immediately 
derivable from either Eqs. (5) or (7) and may 
well be connected with the anisotropic thermal 
expansion of the unit cells of a liquid, as observed 
by x-ray diffraction.* Some specific cases of rela- 
tionships between structure changes of a liquid 
and AHt will be discussed later on.** 

° Eirich and Simha, J. Chem. Phys. 7, 116 (1939). 

* Eyring’s other relation AEyap/AF:=2.45 was purely 
fortuitious as it was observed in the course of this work 
that AF; passes often through a minimum and rises then 
constantly until Toit; AEvap/AFt taking on values ranging 
from 4 to 0. 

7A. Bondi, ACS-Meeting in print, 1945, Petroleum 
Division, p. 315. 

’ T. Alfrey, Rheol. Bull. 16, No. 1 (1945). 

° E. W. Skinner, Phys. Rev. 36, 1625 (1930). 

** The anomalous increase of AH: with temperature 


which has so far mostly been observed for low boiling 
liquids (see h—cs on Fig. 5a, and the lower alcohols on 


The other component of AH}, the energy of 
motion AH}! has a number of interesting proper- 
ties which point toward an understanding of its 
nature, and which have so far been entirely 
overlooked: (a) the strong temperature de- 
pendence and (b) the close relation between AH}? 
and the ‘“‘non-ideal”’ portion of AE yap. 

The first named of these properties is peculiar 
to heats of activation connected with the ac- 
tivation of internal degrees of freedom.'!® The 
number S/2 of degrees of freedom activated is 
to (a not very good)f approximation 


0AH? S 
(2) 4). 
OT Jp 2 


and the average characteristic frequency vt of 
these degrees of freedom could be given by 


S 
— (9) 


The values of S/2 and yvt¢ thus calculated are 
assembled in Table II. One sees that S/2 
increases quite regularly with molecular weight, 
as would be expected, and that v¢ is of the order 
of Raman shifts and infra-red bands ascribed to 
external and internal rotation.'' The motion of 
larger molecules past each other seems to require 
the activation of normally (at that temperature) 


TABLE I. The heat of activation for hole formation AH?" 
(in kcal./mole) calculated by Eqs. (3) and (5), and 
its relation to the energy of vaporization. 











AHt* from equation AVt AHt* AHt* 
(5) V AEy®* AEy 





CCl 
n-hexane 
n-octane 
n-decane 
Toluene 
n-propanol 
n-pentanol 
Oleic acid 


0.244 0.237 
0.105 0.100 
0.101 0.090 
0.097 0.090 
0.165 0.155 
0.130 0.090 
0.168 0.115 
_— 0.125 


Be Nas 
Neem ooo= 
SSSSSSzE 


NrFOorooo- 
eEsksesi |e 
SSyssasi\~ 








* AEy® =AEy —AEyi. 


Fig. 7) may also be connected with the fact that the con- 
centration ¢c; of activated molecules starts in all of these 
cases to rise above about 1 percent in the same tem- 
perature range. Mutual interference of the activated zones 
may, therefore, be responsible for these peculiarities. 

10 TP. D. Eley, Trans. Faraday Soc. 39, 168 (1943) and 
private communications. 

+ See discussion between R. M. Barrer and D. D. Eley 
in Trans. Faraday Soc. 39 (1943). 

1G. Herzberg, Molecular Spectra I1 (D. Van Nostrand 
Company, New York, 1945), p. 496 seq. 





TABLE II. The internal energy contribution AH?’ to the energy of activation for viscous flow. 
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Neo-pentane 
n-hexane 
n-octane 
n-decane 
Benzene 
O-xylene 
Methyl alcohol 
Ethyl alcohol 
n-propyl alcohol 
n-amyl alcohol 
Acetone 
Carbon tetrachloride 


Ethyl iodide 


AH} 


keal./m 


0.690 
0.805 
1.060 
1.610 
0.462 
0.695 
1.810 
2.430 
3.200 
3.910 
0.570 
0.660 
0.780 


0/30 


AHti 
keal./m 


at 
< 





AEyi 
keal./m 





0.0 

0.570 
0.670 
1.250 
0.0 

0.115 
1.650 


2.850 
3.550 
0.400 
0.433 
0.635 


50/100 
40/75 
40/75 
40/75 
50/75 
60/75 
40/75 


40/75 
40/75 
40/75 
40/75 
40/75 


0.000 
0.450 
1.040 
0.740 
0.355 
2.030 
2.340 
3.320 
3.670 
0.730 
0.195 
0.870 
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TABLE IIA. Values of AH+, AH?’, AHz* and AV:/V of several petroleum-fraction (lubricating oils). 








Oil 





t 
°C 


AHt 
kcal. 


Ati 
kcal. 


v 





A-285 


A-409 


N-364 


N-393 


N-435 


N-666 


38/54 
54/99 


38/54 
54/99 


38/54 
54/99 


38/54 
54/99 


38/54 
54/99 


38/54 
54/99 


8.0 
6.30 


14.80 
11.30 


7.90 
6.77 


10.50 
7.80 


11.63 
9.44 


12.80 
10.30 


4.0 
2.43 


ust 
Ne 
tO — 


mw 


wor 
— OD + 00 ne lo oo) 
Who = OO “IW 


et ee 
00 Ge 


43 


68 


64 


49 


AHt* 
AEy 








0.180 
0.180 


0.235 
0.190 


0.145 
0.146 


0.167 
0.185 


0.177 
0.160 


0.125 
0.130 


0.100 


0.167 
0.109 


0.169 
0.126 


0.087 
0.079 





Waterman analyses of oils used: 


Rings per molecule _ 
C-atoms in paraffin side chains 


N-364 


N-393 


N-435 


N-666 





2.48 


15 


2.75 
17 


3.13 
19 


3.40 


34 








Notes: n =average number of C-atoms per molecule (note the fairly good proportionality between ” and S/2). AEy was estimated from the 
heat of evaporation data of “‘similar’’ pure hydrocarbons and corrected to the temperature of viscosity measurement using ACp(liq-gas) =0.2 Cp liq. 
which has been found a rather good approximation, for hydrocarbons of this molecular weight range. 

The constants of this table have been computed from viscosity data by R. B. Dow, J. S. McCartney, and C. E. Fink, J. Inst. Petr. 27, No. 213, 


301 (1941). 


still unexcited external rotation and sometimes 
even internal rotation, i.e., actual deformation 
of the molecule. There seems to be a way to dif- 
ferentiate between the two, and this involves the 
second of the properties of AH} mentioned 
above, the relationship to the non-ideal portion 
of the energy of vaporization. According to 
Hildebrand” the ideal entropy of vaporization 
at a vapor concentration of 0.0202 mole/1 is 
20.2 e.u. The excess entropy of vaporization ASy’ 
leads to the non-ideal portion of the energy of 


2]. H. Hildebrand, J. Chem. Phys. 7, 233 (1939). 


vaporization AEy’ by 


AEy? = T(ASyay — 20.2) = TASy’. (10) 
Table II contains a number of AEy’ values, com- 
pared—where the Hildebrand-temperature falls 
within the range of the viscosity measurements 
—with AH;’. The close similarity between A//:’ 
and AEFy’ is quite striking. Some of the incon- 
sistencies in the homologous series are probably 
caused by inaccurate viscosity and/or vapor 
pressure data. 








0.120 


0.163 
0.145 


0.160 
0.100 


0.167 
0.109 


0.169 
0.126 


).087 
).079 


|-666 
3.40 
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If we accept the suggestion by Halford* that 
ASy? is closely related to the length of free swing 
required for rotational motion in a liquid, we may 
say that whenever AH}/<AEy’ the activated 
rotation is a normally unexcited or restricted 
external rotation.* When AH}’>AEy’ as in the 
case of neo-pentane and other more complex 
hydrocarbons, we must assume AH}? to be 
caused by the activation of internal rotations, 
i.e., of the individual motion of molecule seg- 
ments past (internal) potential energy barriers. 
As restricted internal rotation would hardly be 
affected by change from the liquid to the gaseous 
state, it would not contribute to ASyap. 

In either case AH} arises from insufficiency of 
the available free volume for the molecule to 
execute free rotational motions. Considering that 
AH; decreases to zero for non-polar substances 
at elevated temperatures, it appears reasonable 
to introduce the postulate that the molecules of 


[ 


a liquid will encounter no obstacles to free 
motion whever the effective free volume ex- 
panded to a certain fraction, V,°/V, of the 
molecular volume. It is here proposed to write, 


therefore 
| b3V;\3 
ais'=B{1-(=—) , (11) 
VP 


where E,=energy of the restricted rotator, 
6;=free angle ratio, V;=free volume. Discussion 
of this equation is beyond the scope of this paper 
and will appear at a later date, suffice it to say 
that no unreasonable assumptions are required 
to give a good fit of the experimental data. 

To sum up: AH} is more complicated to derive 
from first principles than AH", and can in many 
(but not all) cases be predicted—for one tem- 
perature—from the excess entropy of vaporiza- 
tion. [Its temperature coefficient can in the case 
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Fic. 2. Relation between AS; and AH;/T at T= 303°K for aliphatic hydrocarbons @), 
alcohols (3), and acids (2). Slope of broken line = 1.0. 


* R. S. Halford, J. Chem. Phys. 8, 496 (1940). 


*In this case Eq. (8) would either have to be thought of as applying to more than one molecule, or to be aban- 


doned entirely. 
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of hydrocarbons be very approximately com- 
puted from the number of C-atoms per molecule. ] 
Where strong dipole forces are active AH} 
changes but little with temperature. 


2. THE ENTROPY OF ACTIVATION 


In most previous investigations little attention 
has been paid to the entropy of activation for 
viscous flow ASt, since its numerical value 
deviates but little from zero for most of the 
“common”’ liquids. The more viscous materials 
which have been examined in the course of this 
study show, however, AS} values of between — 10 
and +100 e.u. which cannot possibly be ignored. 
One can observe fluids of high AH; to show also 
large ASt values. There is, indeed, within any 
one group of compounds a degree of propor- 
tionality between AS: and AH; (Fig. 2) which 
reminds one of the Barclay-Butler rule con- 
necting the heat and entropy of vaporization. 
Frank suggested’® that the Barclay-Butler rule 
has the following meaning: ‘‘When a vapor con- 
denses, a molecule, so to speak, falls into a poten- 
tial well, i.e., the greater the loss of energy on 
condensation, the narrower it is and the smaller 
the free volume, therefore the greater the loss of 
entropy. In addition, the deeper and narrower 
the well, the greater the interference tends to be, 
not only with moving back and forth, but also 
with squirming into different positions. This also 
means loss of entropy.’ Applying this to the 
process of activation for flow, we should expect 
the same partition function changes to con- 
tribute to ASt which contributed to AH;. The 


Fyc. 3. Alfrey’s model; rotational flow. 


47, M. Barclay and J. A. V. Butler, Trans. Faraday 
Soc. 34, 1445 (1938). 
% H. S. Frank, J. Chem. Phys. 13, 793 (1945). 
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greater group specificity of the ASt vs. AH; rela- 
tionship points to the existence of a factor, which 
is specific for the flow process and does not affect 
vaporization, the cooperation of molecule pairs 
as depicted by Alfrey’s* model of the flow (Fig. 
3). The probability for two of the few molecules 
of extra energy content AH; to collide and align 
“properly” for this minuet is certainly very small 
and bound to contribute a sizeable negative term 
to AS}. This term and the term for intra- 
molecular cooperation lumped together as AS;¢ 
would be the only sources of entropy change once 
AH;:/T is zero. The intercept of the AS; vs. AHt/T 
curves with the AS} axis would, therefore, be AS’. 

As would be expected, ASt¢ is constant at low 
temperatures but increases at high temperatures 
reaching relatively quickly another constant 
value, which is then identical with AS}, sug- 
gesting that at elevated temperatures, the entire 
entropy change is caused by cooperative factors. 
Relationship between AS+* and chemical con- 
stitution will be discussed later. Prediction of its 
magnitude by other than semi-empirical rules 
seems at present not possible. 

Other terms contributing to AS; are: AS;" 
which can be derived from the change in effective 
free volume and would be 


’ 


AS;+*=R In ( 


63 V;+A ~) 
53Vy 


AS+" contributes about 1.5 to 10 e.u. The higher 
values are associated especially with liquids 
composed of spherically shaped molecules, whose 
free volume is usually small and AV} quite large. 
Another term is AS}'°t, which is caused by the 
activation of external rotation. The origin of this 
rotational contribution, the lack of free swing, 
has been mentioned above in connection with 
AHy;'. The same considerations which led to the 
conclusion (in Eq. (11)) that only a fraction of the 
total energy of a rotator has to be supplied for 
the process of activation, depending on the 
available free volume, apply to the magnitude of 
AS+"°t. In addition we have to consider the pos- 
sibility that the flow process might require 
rotation only around a single preferred axis. 
There is some experimental evidence for this 
hypothesis: If we compare the entropy of activa- 
tion for dipole orientation AS+(€) and for viscous 
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TaBLE III. Comparison of the processes of activation for flow and for dipole orientation. 











t 


Substance Process <= 


AHt ASt 
kcal. e.u. 


ASt(e) —ASt(n) S(rot —1)** 
e.u. e.u. 





n-propanol Flow 25 


Orientation 25 


Flow 25 
Orientation 25 


n-butanol 


Flow 30 
Orientation 30 


Cocoanut oil 


Flow* 10 
Orientation 10 


Ethyl abietate 


Flow 20 
Orientation 20 


* Extrapolated viscosity data used. 


4.53 4 
6.10 5 
7 
4 


5.11 6.7 


4.61 75 


40 8.2 


22.7 


23.84 
26.30 


19.15 
18.60 











** Rotational entropy contribution calculated for rotation around longest axis of molecule. 


flow ASt(n), we shall find ASt(€)—ASt() often 
to equal just the rotational entropy contribution 
calculated for rotation of the dipole around the 
longer axis of the molecule (see Table III), a 
motion which is in accord with our present views 
of the rotation of dipoles in an alternating electric 
field. Had ASt(7) contained the complete AS**, 
ASt(e) would have been smaller than ASt+(n) 
instead of larger, as is always observed. Which 
of the three possible external rotations con- 
tributes to AS+(n) will for the time being be a 
matter of trial calculations.** As our knowledge 
of the small (difficulty observed) temperature 
dependent rotational Raman-shifts of liquids 
increases we shall approach the possibility of 
predicting AS;"°t from such spectrographic data. 
Similar reasoning applies to the computing of the 
contribution AS} of activated internal rotation 
(or torsional oscillations). 
From the relation 


F: AFyzinternal 
AS:=R In >; . 


follows that the slope of the curves in Fig. 4 
represents the ratio of the change in internal 
energy (during activation) to AH;. It is in 
accordance with the views expressed above that 
—by these curves—AEi"*/AH; should tend to 
zero at high temperatures (high V;/V ratios). 

To sum up the results of this analysis of ASt: 


**From the work of Kauzmann and Eyring [J. Am. 
Chem. Soc. 62, 3113 (1940) ] it could be concluded that 
AS3°t of the moving segment of long linear molecules is of 
the order of 7 e.u. 


It is found that of the many terms contributing 
to AS only a few can be separated, let alone be 
predicted. By analysis of viscosimetrically deter- 
mined data one can separate the contribution of 
cooperation factors ASt*, and the translational 
term AS;**. In the future AS" and AS}? may 
become predictable from Raman or infra-red 
spectra in the 10 to 100 cm™ range. 

We must admit that the analysis of AH; and 
AS: which—it is believed—has here been tried 
for the first time has failed of its objective, the 
prediction of the absolute value of the viscosity 
of pure liquids from available physical property 
data, but the results obtained appear to provide 
sufficient new insight into the mechanism of the 
flow process to evaluate profitably the influence 
of molecular structure upon viscosity. 

In the following paragraphs the influence of 
the two most prominent factors, molecular 
geometry and hydrogen bonding, upon AHt, 
AS+, and AV: will be discussed in some detail. 


3. MOLECULAR GEOMETRY AND VISCOSITY 


A number of investigators have tried to cal- 
culate the viscosity or the viscosity-temperature 
coefficient from atomic or bond increments in a 
manner similar to the molar volume, parachor, 
refraction, etc.’* If it was not already obvious to 
a large number of workers in the field, the vast 

16 E. C. Bingham and collaborators, J. Rheol. 1, 372 
(1930); 3, 113, 479 (1932). J. N. Friend, Phil. Mag. 34, 
643, 810 (1943); 35, 57, 136 (1944). H. Umstatter, Kolloid 
Zeits. 102, 232 (1943). Others reviewed in “ist and 2nd 


Report on Viscosity and Plasticity to the Amsterdam 
Acad. Sci.,’’ Amsterdam, 1935 and 1938. 
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Fic. 4. Plot of AS; vs. AH:/T over a wide temperature range for a number of typical compounds. Slope 
of broken line= 1.0. 


@) Polyiso-butylene 

(2) n-triatetracontane 

(3) 10-nony! nonadecane 

(4) di-isobutyl-naphthalene 

(5) tri-(cyclohexyl-ethyl)-methane 


amount of viscosity data of pure hydrocarbons 
which have accumulated during the past decade” 
has made it quite clear that these attempts have 
been futile. The various isomers of the hydro- 
carbon differ often in viscosity to such an extent 
that bond additivity is ruled out entirely. 

On examining the equations of the previous 
sections for the basic causes of this behavior of 
liquids, we find that there are, indeed, a number 
of bond or atom additive properties contained 
in a number of parameters: the polarizability a 
in AH;", the mass per molecule in the rotational 
energy and 63, of AH+’ and AS}'°t and lastly the 
molar volume in the quotient Nh/V. There are, 
however, a number of decisive factors which are 


17 See in M. P. Doss, Physical Constants of the Principal 
Hydrocarbons (Texaco, New York, 1943). 


©) tri cresyl phosphate 

(7) 2,2- di-pseudo cumente-propane 
castor oil 

(9) cutyric acid 

(0) n-amyl alcohol 


TABLE IV. Relationship between the volume requirement 
AV: of the activated molecule and molecular structure 
(all from viscosity/pressure data below 
1000 atmospheres). 








AV} 


Substance > V 





(a) Flexible structures 
n-decane 

di-n-amy] ether 
n-octadecy! alcohol 
Oleic acid 

Oley! oleate 

Linseed oil 

Castor oil 


Petroleum fractions* 


(b) Rigid structures 
Neo-pentane 
Cyclohexane 
Chlorobenzene 
Tetrachlorodipheny]! 
Pentachlorodipheny1! 
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TABLE V. Effect of spatial configuration on viscosity of hydrocarbons and organo silicone derivatives. 








7, 
poise 


AHt 
kcal. 


AVt 
cm? 





0.0022 
0.00148 


0.0023 
0.0011 


0.0029 
0.0019 


0.00375 
0.00334 


n-pentane 
(AS; = 14.0 e.u.) 


Neo-pentane 
(AS; = 2.4 e.u.) 


n-hexane 
(AS; = 17.5 e.u.) 


Neo-hexane 
Ethyl-trimethyl-Silane* 


0.0082 


Cyclo-hexane 
0.0043 


(AS; = 2.2 e.u.) 


0.00783 
0.0054 


0.0448 
0.0176 


0.0131 
0.00541 


n-decane 
(AS; = 28.3 e.u.) 


Cyclodecane 


Decamethy] tetrasiloxane** 


Octamethyl cyclo tetrasiloxane** 0.0220 
0.0070 


n-Co¢ 


C:—C—Cz, 0.20(S) 
| 0.066(S) 
C; 
0.160(S) 


60 0.050(S) 


20 
60 


7.30(S) 
0.45(S) 


0.099(S) 


0.291(S) 
0.055(S) 


0.484(S) 
0.063(S) 


Cs 
LC >-@},¢ 
[ C; : C 
n—Co2 
(Cio)3°C 


80 


20 
80 


20 
80 
245 


G.-C-C-C-C-C-C—C, 
G GG G&G 


0.0656(S) 


0.0079(S) 


15.3 
17.5 


1.30 
34.4 
43.4 
7.3 
0.3 


19.01 
9.35 


4.37 


6.56 
5.06 


7.87 
6.00 
3.89 








* Viscosity data from: F. C. Whitmore et al., J. Am. Chem. Soc. 68, 475 (1946). 
** Viscosity data from: C. B. Hurd, J. Am. Chem. Soc. 68, 364 (1946), and D. F. Wilcock, ibid. 68, 691 (1946). 


not at all bond additive: AV;, which determines 
the magnitude of AH", depends upon the flexi- 
bility of the molecule, and is the smaller a 
fraction of the total molecular volume the more 
flexible a molecule, that is, the smaller a mo- 
lecular segment can move freely and act as unit 
of flow. Table IV gives a number of typical 
examples for the relationship between molecular 
structure and the value of AV+. Attention is 
called to the extreme values for rigid spherical 


structures like the chlorinated diphenyls on the 
one hand and the uniformly small AV¢ of the 
linear molecules such as the fixed oils on the 
other. A remarkable feature of the important 
ratio AV+/V is its independence from the con- 
stituent atoms of the molecule, being a func- 
tion only of the relative distribution in space 
(Raumerfuellung) and—as mentioned above—of 
the flexibility of a molecule’s structure. The non- 
additive elements of AH; are E, and 63, which 
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Fic. 5c. AH; vs. temperature curves of 
polycyclic hydrocarbons. 


(1) 2, 2 dipseudo cumene propane 

(2) 8-dihydro diamyl anthracene 

(8) 6-dihydro dioctyl anthracene Fic. 5. The temperature curves of the heat of acti- 
(4) tri(ethyl cyclohezane) methane vation AH:. These curves show quite well the de- 


é : pendence of the temperature sensitive component 
© 8-dihydro diethyl anthracene ; of AH+; namely, AH?/, upon molecule size and struc- 
tri-octadecyl-benzene (for comparison) tural complexity. 
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TABLE V(A). Effect of position of ring on viscosity. 











t Viscosity AHt ASt 
~ stokes kcal. e.u. 
n-octadecane (for 
comparison) 20 0.0579 3.84 — 3.0 
80 0.0184 3.64 —3.8 
1-phenyl dodecane 20 0.0629 4.36 — 1.6 
80 0.0197 3.52 —4.1 
2-phenyl dodecane 20 0.0688 4.74 —0.4 
80 0.0197 3.82 —3.2 
4-phenyl dodecane 20 0.0782 5.40 +1.0 
80 0.0196 3.85 —3.2 
6-phenyl dodecane 20 0.0880 5.94 +3.5 
80 0.0197 4.30 —1.9 


| 
! 
| 
| 





contain the moments of inertia—and their 
product—being clearly determined by the rela- 
tive distribution of the molecular structure in 
space, and the free volume V; ‘‘which is due to 
the thermal displacements of the center of 
gravity of the molecule from its equilibrium 
position.’’'§ The free volume is thus directly 
related to the expansion at the melting point; 
the curious relationship between viscosity func- 
tions and the entropy of fusion AS; becomes 
thereby somewhat clarified. The data on Table V 
show that whenever AS; is large, AHt and AS} 
are small and vice versa. The behavior of the 
spherical molecules neo-pentane and cyclohexane 
—whose AS; has nearly the theoretical value R 





TABLE VI. Demonstrating the relative independence of viscosity function from chemical composition and the dominating 
effect of molecular structure. 
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—is quite characteristic of this general trend. In 
other words V; is—by AS; and AV;—also essen- 
tially a function of molecular geometry and not 
of chemical character.t| The temperature coef- 
ficient of AH}’—and thereby of AH}—is according 
to Eq. (11) determined by the expansion coef- 
ficient and V;. The former is within any one 
group of compounds proportional to the re- 
ciprocal molecular weight, and varies quite 
generally inversely as the cohesive energy E”’. 

These relationships are well recognizable in 
Figs. 5a-c. 

The effect of most dipoles, excepting OH and 
NH, on AH; and on viscosity in general, is quite 
small. This is not at all surprising for the size 
of the dipole interaction terms occurring in both 
AH;" and AH} depends very strongly on the 
distance between the dipoles. With the exception 
of those dipoles which contain hydrogen this 
distance of closest approach is so large as to 
reduce the magnitude of the dipole interaction 
contribution to but a few hundred cal./mole, as 
is readily seen on Table VI, where polar mole- 
cules and structurally similar hydrocarbons are 
compared with each other.* It appears from 
these data that one need not consider chemical 
characteristics at all if one wants to synthesize 








Substance M ‘ 


poise kcal. kcal. 





ut 


U AFt 
d.u. kcal. 


sg 
e 





Di-n-amyl ether 158 20 
n-undecane 156 20 
Methyl stearate 298 50 
n-heneicosane 296 50 
di-ethylhexyl sebacate 428 38 
9-octyldocasane 422 38 
Tricresyl phosphate 368 60 
Tri(ethylbenzene) methane 329 60 
Ethylabietate* 330 40 
Di-isoamyl dihydro-anthracene* 322 40 


0.0108 3.66 2.87 —2.7 1.2 0.3 
0.0117 3.74 2.94 —2.8 0 
0.039 5.23 4.11 —3.5 1.8 0.6 
0.037 5.25 4.00 —3.8 0 
115 5.87 5.10 —2.4 3.6 2.4 
15 5.95 5.15 — 2.6 0 
135 6.15 9.03 8.7 (4.3) 2.3 
106 5.97 8.73 8.3 0 
41 7.26 15.30 25.6 1.8 0.6 
.20 7.82 18.7 36 0.5 








* These two compounds differ somewhat in structure. A better hydrocarbon analog of ethylabietate could, however, not be found in the literature. 
The relatively close numerical vlues of their viscosity constant demonstrate quite well the material effect of rigid condensed tricyclic structures. 


t The dipole moments are estimated values. 


's Kincaid and Eyring, J. Chem. Phys. 6, 620 (1938). 


t The relation between AH; and Raumerfuellung, i.e., the three-dimensionality of molecular structure, has first been 


discovered by J. Pirsch, Ber. 70, 12 (1937). 


*It should also be noted that y?/r* is the negative value of the orientation energy contribution to E” at full parallel 
orientation,of the dipoles. This interaction term will tend to zero with increasing angle of mutual orientation. As there 
's obviously always a finite (and quite sizeable) concentration of non-oriented dipole molecules from which c; would be 


recruited, the influence of dipoles on viscosity is, as a rule, quite small; as is brought out by the figures above. 
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Fic. 6. Relationship between AS+° (from AS; vs. AH;/T intercept) and molecule 
size for long chain hydrocarbons (1) and cyclic hydrocarbons (2). 


a substance of given viscosity characteristics but 
only the stereometric arrangement of molecular 
structure (provided there are no OH groups and 
not too many cumulative dipoles in the molecule). 

The cooperative entropy factor ASt* is 
seemingly directly proportional to the size of 
the unit of flow (OH compounds excepted). Fig. 6. 
The balance of the entropy terms depend on the 
same factors as the corresponding energy con- 
tribution. 


4. HYDROGEN BONDED LIQUIDS 


The well-known anomalies caused by hydroxyl 
groups, which have been quite conclusively 
ascribed to hydrogen bonds!* are naturally also 
reflected in the viscosity function. The simplest 
case is represented by the fatty acids. Electron 
diffraction and vapor pressure data®®*! have now 
well established that the hydrogen bond, coupling 
two carboxylic acids to a dimer unit, is so well 
stabilized (by resonance) that the double mole- 
cules exist not only in the liquid but even in the 
vapor state, at least for the lower members of the 

19S. Gladstone, Textbook of Physical Chemistry (D. Van 
Nostand Company, New York, 1943), p. 114. 

* J. Karle and L. O. Brockway, J. Am. Chem. Soc. 66, 
574 (1944). 


27H. L. R. Ritter and J. H. Simons, J. Am. Chem. Soc. 
67, 757 (1945). 


fatty acid series. It is, therefore, but an additional 
confirmation of the reality of the double molecule 
to find that the viscosity parameters of fatty 
acids correspond remarkably well to those of 
hydrocarbons of the structure and molecular 
weight of the fatty acid double molecule, as seen 
from Table VII. Because of the relatively large 
accumulation of dipoles in the double molecule 
of the lower members of the series, the viscosity, 
AH:, and ASt of the latter are somewhat higher 
than would otherwise be expected. The con- 
sequently more pronounced contribution of the 
dipole interaction terms to AH" and AH; cause 
the quite small temperature coefficient of AH: 


TABLE VII. Evidence for the existence of fatty acids as 
double molecules by comparison of viscosity data with 
those of hydrocarbons and esters. 





M — poise kcal. 





0.011 (3.60) 
0.011 3.51 
0.0091 3.54 


(4.28) 
4.22 


Propionic acid ( 20 
Propionic anhydride 130 20 
n-decane 1 20 


Valeric acid 20 
n-tetradecane 198 20 


nD ip 
NAN 


0.022 
0.022 


0.043 (4.88) 
0.045 4.78 


0.318 (6.60) 
0.330 6.30 


nw NNN 


n-heptylic acid 20 
n-octadecane 20 


Oleic acid 20 
Oley! oleate §32 20 
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in ee 
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Data in parentheses are calculated for the double molecule. 
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TABLE VIII. Energies and entropies of activation of aliphatic alcohols. 
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2.79 
2.41 


1-octadecanol 85 
140 


ga 
a) 
oot 


2.0: 


0.192 
0.199 


14 : 0.087 at 70°C 
5 0.092 at 150°C 











Notes: SEyi refers to temperature at which vapor conc. =0.0202 mole/1 AVt/V is always taken at 30 and 75°C, respectively, except where 


noted otherwise AEy® =AEy —AEy/. 


* The extremely high values of AHt and ASt for the symmetrical trialkyl carbinols seem readiiy ascribed to the spherical shape of these com- 
pounds. But only viscosity pressure data—which are still sadly missing—will provide through AVt, AHt/ and AHt* values a full understanding 


of the mechanism of flow of such fluids. 


t The data of these two isomers were inserted to show the interference of hydrocarbon radical structure with OH-network formation. 


and thereby the bent AS; vs. AH;/T curves 
shown in Fig. 9. Progressive dissociation of the 
dimer with increasing temperature would affect 
those parameters in the same direction. 

The state of molecular aggregation of alcohols 
is less well defined’ than that of the fatty acids. 
At ordinary temperatures the bond between 
hydroxyl groups is not strong enough to form 
distinguishable double molecules. They rather 
form a continuous network structure throughout 
the liquid phase. Contrary to ordinary dipole 
systems, for which the potential energy dif- 
ference between the fully oriented and random 
configuration is usually less than 2 kcal., the 
hydroxyl group network is more probable than 


randomness by a potential energy difference of 
more than 10 kcal.” 

The motion of any alcohol molecule involves, 
therefore, passage of OH-groups past each other. 
Since the displacement per jump is not larger 
than with other molecules, AV and AH¢" are of 
normal magnitude, Table VIII, but the motion 
in the strong fields requires a AH;/ considerably 
greater than that of molecules of comparable 
size. As “the jumping pairs’’ now have a pre- 
determined configuration, ASt* of (the lower) 
alcohols is not very important, and AS; becomes 
relatively high, compensating the large AH} 
value somewhat. The viscosity of alcohols is 

#2 H. Harms, Zeits. f. physik. Chemie B43, 257 (1939). 
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Fic. 7. AH: vs. temperature curves of aliphatic alcohols. 


nevertheless quite a bit higher than of other 
compounds of similar structure and molecular 


weight. Low temperature (at —75°C) x-ray 
diffraction data?* indicate the presence of true 
hydrogen bridges (as in fatty acids) of 2.7 to 
2.94 O—H-—O distance in methyl and ethyl 
alcohols. At room temperature these are no 
longer observable. If this should mean that 
alcohols form distinguishable double or otherwise 
multiple molecules at low temperatures, whose 
motion would not involve mutual passage of 
OH-groups, one would expect AH} to decrease at 
the same low temperature instead of increase. 
Calculation of AH}; from the careful viscosity 


% G. G. Harvey, Phys. Rev. 55, 1115 (1939). 


measurements of Tonomura** down to —98°C 
shows this actually to be the case, as depicted 
on Fig. 7. Much more revealing, indeed sur- 
prising, are the AS vs. AHt/T curves of these 
alcohols, shown in Fig. 8. They show the regular 
linear behavior at lowest temperatures until they 
reach a sudden inflection point, where within a 
fractional unit of AH+/T, ASt increases by 1.6 
to 3.0 e.u. Remembering that AS; is in the form 
e4S/k a measure of the population of a given 
energy state, we compute from these entropy 
changes that at the inflection point the number 
of units of flow (but not ct/c!) increases (by 


* T. Tonomura, Sci. Rep. Tohoku Imp. Univ. [1] 22, 
104 (1932). 
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Fic. 8. Alcohols AS vs. AHt/T curves. 


increasing the temperature) by a factor of 2.2 
for methyl-alcohol, 2.4 for ethyl alcohol, and 4.5 
for n-propyl] alcohol. The inflection temperatures 
are at about 10°C, —50°C and —30°C, re- 
spectively. The inflection of the butyl and amy] 
alcohol curves are more gradual, and less well 
defined. If the interpretation of the behavior of 
alcohols is correct, we should expect to observe 
similar effects of the dissociation of fatty acids 
upon ASt and AH; at elevated temperatures. 
There are but few high temperature viscosity 
data of fatty acids available, and dissociation 
of these in the liquid state takes place at tem- 
peratures at which the increase of AS}* starts 
already to flatten out the ASt curves, but never- 
the less the curves on Fig. 9 indicate the general 
trend to be as expected. Viscosity data of com- 
pounds containing more than one OH-group per 
molecule are very scarce. The most peculiar of 
these is castor oil. The AH}? values of this liquid, 
calculated from the viscosity data of Hersey 
and Shore®® are presented on Table IX. The 
stepwise decrease of AH}! with increasing tem- 
peratures suggests successive disengagement of 
the hydroxyl groups with increasing disorder in 
the liquid. 


** Hersey and Shore, Mech. Eng. 50, 221 (1928). 


5. VISCOSITY AND LIQUID STRUCTURE 


The structure of liquids, which by x-ray 
diffraction and their theoretical interpretation 
has acquired the status of a reality today, has 
hardly been referred to in the quantitative 
discussion. It is, however, not accidental that 
the rate process theory of flow should be able 
to give a good account of the relationship be- 
tween molecular structure and viscosity without 
explicit consideration of the structure of the 
liquid. The reason for this lies in the basic 
premise of the theory: the proportionality 
between fluidity and the concentration of 
activated molecules ct. The smallness of ct has 
never before been sufficiently emphasized. We 
have mentioned in Section 1 that ct is small 
enough to be entirely recruited out of the small 
random oriented population even in such 
liquids, which are otherwise known as well- 
ordered systems, such as n-paraffins, fatty acids, 
etc. at ordinary temperatures. The activated 
molecules have since the inception of the theory 
been thought to occur only in the immediate 
neighborhood of lattice disturbances, the “‘holes.”’ 
Just as the frequency of lattice imperfections in 
crystals is not @ priori connected with the 
geometry of the lattice, such condition should 
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Fic. 9. Acids ASt vs. AHt/T curves. 


not be expected to exist in liquids. An indirect 
connection does exist through the proportionality 
between the number of holes in the lattice and 
exp [(—¢"+¢')/kT ], where ¢* is the energy of 
hole formation and ¢‘ the potential energy of a 
molecule on an interlattice position, both of 
which are affected by the number of nearest 
neighbors and thereby the lattice arrangement. 
We shall similarly expect AV: (and by it AH" 
and AS}‘t) to be affected by the number and 
arrangement of the nearest neighbors. Further 
refinement of the theory, designed to predict 
AVt and its temperature and pressure function 
will, therefore, have to take the details of the 
liquid structure into account. 


SUMMARY 


1. The energy of activation for viscous flow 
consists of two independent terms: the energy of 
hole formation AH; and the energy of motion 
into the hole AH;’. AH;* has here been derived 
as the change in potential energy between one 


molecule and its nearest neighbors, which accom- 
panies the displacement from equilibrium posi- 
tion incident to the formation of a hole. With the 
displacement given by the pressure coefficient 
of the free energy of activation (@AF:/dp),=AVt 
London’s interaction energy law is shown to 
predict AH;" values of the right order of mag- 
nitude. AH}, the heat of activation of the vis- 
cosity isochore is shown to be often of the same 
magnitude as T X Hildebrand’s excess entyopy 
of vaporization, a parameter which might also be 
called the non-ideal portion of AEyap. Both may 
be considered as a measure of restricted external 
rotation in the liquid state. Where AH’ is 
materially larger than AEyap’, the molecular 
structure as well as the temperature coefficient 
of AH?’ suggest that the process of motion 
requires deformation of the molecule against 
internal potential energy barriers. 

2. The entropy of activation ASt is separated 
into a cooperative term AS}* which is derived 
graphically from the AS: vs. AH+/T curve, a 
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translational term AS;*" which is calculated from 
the volume increase A Vt and the free volume, and 
a rotational term AS;"°t, which represents the 
entropy counterpart of AH;’. 

3. Numerical examples are provided to show 
that the viscosity of liquids depends largely on 
the stereometric arrangement of molecular struc- 
ture and only secondarily on the nature of the 
constituent atoms and resulting dipole moments 
(OH groups excepted). This result is shown to 
be caused by the determining influence of the 
purely geometrical factor AVt upon AH}* and 
AS:**, and the presence of the product of the 
moments of inertia in AH;/ and AS}"°t. These 
findings are believed to be significant for the 
planning of the syntheses of substances with 
given viscosity characteristics. 

4. Viscosity data confirm the existence of 
fatty acids as double molecules in the liquid 
state at room temperature. The viscosity be- 
havior of alcohols suggests the presence of dis- 
tinguishable multiple molecules at low tem- 
peratures (below —50°C), and quite well-defined 
transition from this form to the continuous 
OH-bond network structure which prevails at 
ordinary temperatures. 

5. It is pointed out that in all ordinary liquids 
the concentration ct/c of activated molecules is 


TABLE IX. AHt?/ values of castor oil. 








eC: 20/40 40/60 60/80 80/100 100/133 
AH; kcal./mole: 8.75 8.28 5.33 5.34 2.60 








so small (10~* to 10-*) that the random oriented 
portion of otherwise well-ordered liquids can 
usually provide their entire number. This causes— 
contrary to expectation—the structure of a 
liquid, i.e., its state of order, to play only a minor 
role in the calculation of viscosity. 
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Absorption Measurements at 15 Megacycles* 


J. R. PELLAM AND J. K. GALT 
Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received July 8, 1946) 


Equipment developed by the Massachusetts Institute of Technology Radiation Laboratory 
has been applied to the measurement of sound velocity and absorption in liquids at 15 mc/sec. 
Pulses of one microsecond duration are generated by a transducer, which also picks up the 
resultant echoes from a plane reflector. Velocity measurements are made by determining the 
distance the transducer must be moved to delay the received echoes by a specified increment. 
Absorption measurements are made by determining the attenuation necessary to keep the 
receiver signal constant as the transducer is moved. The attenuation factor can be measured 
to an accuracy of about 5 percent and sound velocity to about 0.05 percent. Measurements 
in homologous series of organic liquids are reported. 





I, INTRODUCTION 


HE electronic pulsed-circuit techniques and 

equipments developed during the war are 
adaptable to the measurement of velocity and 
absorption of sound in liquids.! Preliminary work 
of this nature took place at the M.I.T. Radiation 
Laboratory,? and the pulse principle has also 
been applied’ to the examination of mechanical 
flaws in metals. For velocity determination in 
liquids the pulse method is at least as accurate 
as optical diffraction and acoustic interferometer 
methods previously employed,‘ and considerably 
more accurate for absorption measurements. It 
is direct and convenient and therefore particu- 
larly adapted to investigations involving large 
numbers of liquids. The present paper reports 
measurements of sound velocity and absorption 
in selected organic liquids at 15 mc/sec. and at 
several temperatures. 

The scheme is essentially to use the liquid 
sample as a ‘‘storage medium”’ for short sound 
pulses and to measure the time delay and attenu- 
ation undergone by the sound in traversing a 
known path within the liquid. The acoustical 
pulses are generated from electrical pulses by 


* This paper is based on work done for the Office of 
Scientific Research and Development under Contract 
OEMsr-262. 

1 This has been reported at the April 1946 meeting of 
the American Physical Society: J. R. Pellam and J. K. 
Galt, abstract A12. 

2M. Cefola, M. E. Droz, S. Frankel, E. M. Jones, G. 
Maslach, and C. E. Teeter, Jr. Radiation Laboratory, 
M.I.T. Report 963, March, 1946. 

3F. A. Firestone, Metal Progress 48, 505 (1945). E. N. 
Simon, Metal Progress 48, 513 (1945). 

4L. Bergmann, Der Ultrashall (Edwards Brothers Re- 
print, Ann Arbor, Michigan, 1944). 


means of a crystal transducer and are converted 
back to electrical form upon completing a transit 
through liquid; the effects of changing the acous- 
tical path length can be compensated for electri- 
cally. Hence the increased delay produced by an 
increase in path length gives a measure of sound 
velocity; the attenuation which must be removed 
from the electrical circuit to balance acoustical 
losses in the additional distance provide a meas- 
ure of absorption. 

The variation in acoustic path length is ac- 
complished by mounting the crystal transducer 
on a movable support riding within a tank con- 
taining a sample of the liquid. The 15 mc/sec. 
sound pulses, after being generated from electri- 
cal pulses by the transducer, travel through the 
liquid to a plane reflector and back again to 
re-excite the crystal at a later time. The resulting 
delayed electrical signal passes through a cali- 
brated attenuator to a receiver, the output of 
which is fed to an oscilloscope equipped with a 
special sweep indicating time delays; the start 
of this sweep is synchronized with the generation 
of the original pulse, and calibrated accurately 
with respect to time. The change of the signal 
position on the screen as the transducer mount 
is moved with respect to the reflector is a direct 
measure of the corresponding delay; simultane- 
ously, attenuation is inserted or removed by 
means of the calibrated attenuator to keep the 
receiver input signal constant. 

Velocity is obtained directly from the slope of 
distance traveled plotted vs. pulse delay, and ab- 
sorption is obtained from the slope of compen- 
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sating electrical attenuation vs. distance. Since the 
method depends only upon observing the effects 
of differences in acoustical path length, knowl- 
edge of exact distance traveled by the sound is 
unnecessary. 

Results obtained in this manner may be re- 
garded as physically equivalent to those obtained 
by continuous wave methods, since it may be 
shown that differences in behavior between 
pulsed sound and continuous sound in liquids 
should produce negligible effects. For example, 
no significant dispersion in liquids has been 
detected at these frequencies‘ (nor any pulse 
distortion during these measurements) so that 
the group velocities measured by pulsing should 
not differ measurably from the phase velocities 
ordinarily dealt with. Similarly, although atten- 
uation of sound in liquids varies as the square of 
the frequency, the band width associated with 
the pulses used is small enough to prevent 
appreciable effect on absorption measurements.* 


Il. EQUIPMENT 
(a) Electrical Equipment and Circuits 


Figure 1 is a block diagram of the equipment 
and Fig. 2 is a photograph of the assembled 
apparatus. The timing sequence is as follows: 


* The exponential pressure attenuation coefficient (a), 
defined by 
b= poe ** 


will have a fractional error given by 


Aa _ (~)" Pa _ (2) 

a v] adv 2v 
if measured over a frequency range of width (Av). The 
last term results from placing (@) proportional to frequency 
squared, as is frequently the case, It may be shown that, 
for pulses of much longer duration than the period of a 
single oscillation, Av/2» is just the reciprocal of the number 
of cycles per pulse. The pulses employed in the measure- 
ments reported here contained about 15 oscillations, corre- 
sponding to an error of roughly one part in 250, 
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A crystal controlled circuit in the oscilloscope 
starts the scope sweep and simultaneously sends 
a trigger to the pulse generator. The pulse formed 
by the generator passes through a simple im- 
pedance matching network to the transducer; 
when the resulting sound pulse returns to the 
transducer after its round-trip within the liquid, 
the delayed electrical pulse is formed which 
passes back through the matching network. 
Since the transmitter is off when the echo signal 
arrives, very little leakage loss occurs and the 
main portion travels through the attenuator to 
the receiver. 

All triggering and time measurements are 
accomplished by the crystal-controlled circuit in 
the oscilloscope (DuMont Type 256B A/R 
Range Scope), so that synchronization occurs 
automatically. In addition to setting off the 
pulse generator and starting the sweep, this 
circuit also provides marker pips spaced accu- 
rately every 12.192 microseconds thereafter 
(corresponding to 2000 yards radar range). 
Minute examination of the front edge of the echo 
signal is possible by means of an expanded sweep 
with variable delay and a speed of 1 inch per 
microsecond. Crystal control is necessary pri- 
marily for timing the marker pips accurately; 
one pip is developed for each oscillation, and 
since the crystal frequency is accurate to better 
than 1 part in 10,000 the same is true for the 
range marker positions. 





Fic. 2. Apparatus. Left: A/R range scope; center: 
liquid delay tank; and right: pulse forming network, 
oscillator, attenuator, and receiver. 
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Fic. 3. Liquid delay tank. 


The detailed discussion of the design of the 
pulse generator circuit has been given elsewhere.5 
The trigger from the oscilloscope passes through 
a stage of amplification to excite the second stage, 
a blocking oscillator. The output of the blocking 
oscillator, after being clipped to form a video 
pulse, is applied through a switch tube to the 
15 mc/sec. r-f oscillator. The latter is driven into 
oscillation for the duration of the pulse, and the 
resultant r-f pulsed signal feeds through two 
stages of amplification. The final signal level 
put out to the transducer may be controlled by 
means of a variable power supply to the amplifier. 

Since the electrical impedance presented by 
the crystal is very different from the 70-ohm 
characteristic impedance of the transmission 
line leading from the pulse generator, an imped- 
ance matching network® is provided to improve 
the response of the transducer. A small coil is 
used both as an autotransformer to match the 
line and as a means for tuning out the dead 
capacity of the crystal to give a purely resistive 
load (see Fig. 3). A coil of about 3 microhenries 
total inductance with a 6 to 1 turns ratio fulfills 
both requirements. 

The receiver includes: (1) four stages of ampli- 
fication at 15 mc/sec.; (2) a 15 mc/sec. local 
oscillator and mixer; (3) two stages of amplifica- 
tion at 30 mc/sec. (each stage consisting of three 
tubes with circuits tuned to slightly different 

5S. Frankel, 
645-8. 

6 The authors are grateful to Messrs. Jones and Frankel 


of the Radiation Laboratory for the design constants used 
in this network. 


M.I.T. Radiation Laboratory Report 


frequencies for broad banding); (4) a detector; 
(5) two stages of video amplification; (6) a video 
amplifier in the oscilloscope. The 30 mc/sec. 
stages were used only because of equipment 
availability; a receiver operating at 15 mc/sec. 
throughout is now under construction. 

The calibrated attenuator has steps of 1 db 
and up. Shielding is adequate to prevent leakage 
from introducing appreciable errors, even though 
attenuations as high as 110 db are inserted. 


(b) Sound Delay Tank for Liquids 


A sample of the liquid being investigated is 
placed in a sturdy brass tank (A) (see Fig. 3) 
within a constant temperature bath (B). The 
bottom of the tank is sufficiently thick (3’’) to 
accommodate a milled track (C), with adjustable 
side plate for tightening, and upon which the 
transducer mount (D) travels. Motion is pro- 
duced by turning the phosphor bronze microm- 
eter screw (£) which threads through a stainless 
steel nut (F) set in the transducer mount. A 
hand-wheel (G) is provided and a dial (/7) which 
indicates travel in thousandths of an inch. Cali- 
bration of the screw against size blocks makes it 
possible to read the transducer position signifi- 
cantly to 0.0001”. 

The sound pulse (/) is reflected at the end of 
the tank by a disk (J) of stainless steel, ground 
and lapped optically flat. This reflector is sup- 
ported by three adjustment screws (K), and 
may thereby be set perpendicular to the sound 
beam by aligning for maximum returned signal 
strength. 
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The transducer element (L) is an X-cut 15 
mc/sec. quartz crystal plate, }’’ in diameter and 
gold plated on both sides. This crystal is held in 
place within a cartridge (M) by a Lucite bead 
(NV) so that a diameter of about 1 cm is exposed 
to the liquid. A metal foil (O) between the bead 
and the quartz establishes electrical contact with 


the back surface. 


(c) Temperature Control 


Temperature equilibrium is maintained by 
balancing the effects of dry ice in the bath 
against the warming action of a thermostatically 
controlled electric heater. A water pumping sys- 
tem improves temperature uniformity outside 
the measuring tank, and a motor-driven stirrer 
(P) maintains constant temperature within the 
liquid sample itself. Temperature is read by 
thermometer (Q) and must be maintained to 
within less than 0.03°C for velocity measure- 
ments accurate to 0.01 percent. 


Ill. MEASUREMENT TECHNIQUES 


The velocity of sound within most organic 
liquids is of the order of 1.5 10° cm/sec. which 
corresponds at this frequency to about 100 
wave-lengths per cm. Hence the crystal radiating 
surface is about 100 wave-lengths across, so that 
the pulses travel out within a very sharp beam 
(0.7° to the first minimum). It may be shown 
that this spread produces less than 4 db of 
geometrical attenuation over the range em- 
ployed; therefore readings of attenuation are 
suitable, without corrections, for computing 
liquid absorption coefficients directly. 

Velocity readings are taken at 2000-yard 
(radar) range intervals by placing the returned 
signal pip at definite positions with respect to 
successive scope range markers, and recording 
the setting of the micrometer screw. By properly 
adjusting the scope these markers appear as 
very short (dark) breaks in the sweep, and if 
consistent positioning criteria are adopted, accu- 
tate range difference readings result. It turns out 
that careful adjustment of attenuation is neces- 
sary for consistent velocity readings. Conse- 
quently simultaneous readings for velocity and 
attenuation are preferable. 

This interdependence between measurements 


exists because the signal pulses never build up 
instantaneously, even though they are sharp; 
as a result the front edge of the pulse may 
actually be shifted with respect to the range 
marker by varying the attenuation level. There- 
fore by adjusting attenuation to bring the signal 
peak always to the same arbitrary height on the 
screen (chosen within the linear range of scope 
amplification), the identical portion of the signal 
may be obscured by successive range markers 
Once this condition is satisfied, the crystal posi- 
tion may be set more precisely if a predetermined 
attenuation is then arbitrarily removed; this 
allows setting on a lower, and therefore steeper, 
portion of the pulse front. Consistent results are 
obtained in this manner providing the same 
quantity is removed for each range reading, and 
a value of 10 db was found convenient. 

It is evident that if this adjustment of the 
returned signal to a pre-determined level is 
attained by means of the attenuator leading to 
the receiver (see Fig. 1) without disturbing the 
gain of the remainder of the system, absorption 
data follow as a by-product of the velocity 
measurements. 

For accurate absorption readings the sharpness 
of the beam requires that (1) the reflector disk 
be adjusted carefully, and (2) temperature gradi- 
ents within the liquid sample be avoided. The 
proper operation of the stirrer (P) is absolutely 
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Fic. 4. Plot of sound attenuation vs. distance. 
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TABLE I, Experimental values of the pressure attenuation and velocity in liquids. 
Frequency f=15 mc/sec. p= po exp [— ax]. 


































(1) (2) (3) (4) (S) yar) = ™ 
Temp. Velocity Attenuation (a/f*) X10" sec?/cm =" pg “ 
Substance i * cm/sec. cm! experimental classical oes ont 
ure bromide 2.0 0.905 X 105 0.684 (41%) 304 --- _- 
3Br 
Ethyl bromide 2.0 0.959 0.136 (45%) 61 10) \ 
CH;CH2Br 10.0 0.932 0.14 62 = +0.010 { — 0.0036 
20.0 — 0.14 62 10 
30.0 — 0.18 80 —} 
n-propyl bromide 2.0 1.0326 0.088 (+5%) 39 11 — 
G 3CH2CH-Br 
n-butyl bromide 2.0 1.083 0.11 (410%) 49 13 -- 
CH;(CH2)2CH2Br 
Methy] iodine 2.0 0.888 0.555 (41%) 247 10) ) 
CH;I 6.5 — oT 262 —| — 0.0031 
20.0 — 0.711 316 — 
30.0 0.815 ont " 334 o +0.0087 
35.0 —- 765 340 — 
40.0 — 0.776 345 —| 
— ae 2.0 0.9235 0.09 (410%) 40 12 — 
sCHel 
n-propyl iodide 2.0 0.982 0.122 (+5%) 54 14 -- 
3CH.CH2I 
n-butyl iodide 2.0 1.022 0.109 (4+5%) 48 17 
CH;(CH2)2CHel 
n-propy! chloride 2.0 1.158 0.095 (410%) 42 8 
CH;CH.,CH.Cl 
n-butyl chloride 2.0 1.213 0.242 (42%) 108 10 






CH;(CH:2)2CH:2Cl 
Methyl acetate (Pract) 19.0 1.195 0.16 (43%) 71 6 















CH;COCH; 
Ethyl acetate (Pract) 
H;COC;H; 24.3 1.145 0.163 (43%) 72 8 = 
Amy] acetate (Pract) 























Methyl alcohol 2.0 1.186 0.10 (415%) 45 15\ —0.011 — 0.0032 
CH;0H 19.3 1.123 0.084 (410%) 37 14} 
Ethyl alcohol 2.0 1.246 0.052 (410%) 23 29 
CH;CH,0OH 
n-propyl! alcohol 2.0 1.288 0.195 (42%) 87 55\ — 0.0028 
CH;CH:CH.OH 27.5 1.1976 0.157 (43%) 70 36) _0,0086} 
n-butyl alcohol 2.0 1.324 0.240 (45% 54 65 
CH;(CH:2)2CH,OH 
n-amyl alcohol 2.8 1.320 0.348 (42%) 155 113) \ 
CH:(CH;);CH:;OH 286 —*1.224 0.231 106 58/ —0.014 f  —0,0029 
Benzene 7.5 1. ? 9\ 
C.He 21.8 1.3116 1.79 (41% 800 8 +0.011 — 0.0036 
Nitrobenzene a0 1.518 0.165 (45%) 73 ia 
NO:C.H;5 23.8 1.462 0.178 (+4%) 79 14 +0.0048 — 0.0023 
Chlorobenzene is 1.333 0.276 (41%) 123 10 






Ce 
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n-hexane 0.174 (45%) 

CH;(CH2)«CH; 
n-heptane 22.4 1.150 0.18 (410%) 80 10 
CH;(CH2)sCHs 













Carbon tetrachloride 
Ch, 23.5 0.924 1.20 (4+1% 533 20 
























essential to the temperature uniformity required for level and non-linear absorption curves result. 
accurate absorption measurements. If the stirrer is The heat developed by attenuation within the 
stopped there is an immediate change in signal beam contributes to these refraction effects. 
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IV. RESULTS 





(a) Linearity of Absorption Curves 


Curves of attenuation vs. range are givenin 
Fig. 4 for three sample liquids, indicating the 
degree of linearity found (over a range exceeding 
100 db for CCl,). No effects attributable to 
geometrical divergence of the beam are evident; 
the limit of accuracy appears to depend upon 
the smallest attenuation increment used (1 db) 
as suggested by the small systematic variation 
about a straight line. 






(b) Summary of Measurements 








Part I of this paper is concerned only with the 
presentation of measurements in the form of raw 
data, interpretation being reserved for Part II. 
However, these results are given in terms of 
quantities or parameters possessing physical 
significance in view of Part II. 

Table I is a summary of measurements ob- 
tained for the various series of related organic 
compounds investigated. The velocity is given 
in column (2) in units of cm/sec. and the ab- 
sorption coefficients (@) in column (3) in units 
of cm~!. (a) is the exponential pressure attenua- 
tion coefficient defined by 










p=pe™. (1) 






Since for a number of liquids (a) varies nearly as 
the frequency squared, a convenient form for 
comparing measured results with those computed 
from classical (viscosity) theory is given in 
column (4). The classical absorption values were 
calculated considering viscosity effects alone; 
heat conduction will increase these values by 
about 5 percent. Columns (6) and (7) give the 
experimental temperature coefficients for (a) and 
velocity, respectively. 

These particular series of homologous com- 
pounds were chosen with a view to observing the 
gross effects of substitutions within molecules 
(such as the revelation for example that replacing 
an —H in benzene by a —Cl cuts absorption by 
a factor of 5, and that an —NOz cuts it by 9). 
The most complete array examined was the 
alkyl halide group, of which 10 members were 
included (the 2°C was chosen just below the 
boiling point of methyl bromide). 
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Fic. 5. Sound pressure absorption coefficient a 
vs, temperature. 






(c) Accuracy 


Velocity results are given here to an accuracy 
of about 1 part in 1000. It is believed, however, 
that the pulse method is inherently capable of 
yielding as high accuracy for velocity measure- 
ments as the best diffraction techniques. On the 
other hand, whereas previous techniques were 
often incapable of exceeding 25 percent accuracy , 
in absorption measurements, the absorption is 
given here to about 5 percent on the average. 
Since the limit of this accuracy is set by the 
smallest attenuation increment, best results 
occur for highly absorbing liquids since they 
provide a greater attenuation range. Thus for 
liquids such as CCl, (a) is given to about 1 
percent, but is accurate only to about 10 percent 
for liquids exhibiting low absorption. 


(d) Temperature Dependance of Absorption 


The pulse method for measuring (a@) is suffi- 
ciently accurate for investigating the effects of 
temperature. For example Fig. 5 shows the 
absorption coefficient for methyl iodide plotted 
against temperature (for values up to the boiling 
point). Temperature coefficients may be obtained 
from such curves and, as shown in Part II of this 
paper, have considerable significance with regard 
to the mechanism of sound propagation. 


V. REMARKS 


(a) Anticipated Measurements 


A new apparatus for measuring velocity and 
attenuation by the pulse technique has just been 
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completed. It is designed for operation at liquid 
helium temperatures. Preliminary measurements 
have already been made for liquid nitrogen and 
argon, and measurements in liquid helium will 
shortly be under way. 
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The free volume theory of liquids is applied to problems connected with sound propagation. 
The equation of state derived by Tonks for assemblies of hard elastic spheres may be written 
PVa=3.NkT, where the available volume Va is the difference between the actual volume and 
the volume corresponding to the closest possible packing. It is shown that the molar available 
volume is approximately given by the increase in molar volume of a substance on heating from 
0°K to the temperature considered, so that the pictorial concept of the available volume has a 
real physical basis. Using the Tonks equation, expressions are derived for the sound velocity, 
and the temperature and pressure coefficients of sound velocity. These expressions are com- 
pared with experimental results. Expressions are also given for the cubic expansion coefficient 
and the difference of the specific heats on the free volume model. The Herzfeld-Rice derivation 
of the sound absorption caused by phase shifts in the vibrational specific heat is carried through 
for the Tonks equation of state. The temperature and pressure coefficients of sound absorption 
are discussed. Remarks are made on the relation of molecular structure to sound absorption. 
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I, INTRODUCTION 


HE purpose of this paper is to describe the 
velocity and absorption of ultrasonic waves 

in liquids in terms of the free volume theory of 
liquids and the relaxation theory of absorption 
processes. It is emphasized at the beginning that 
simple models of liquids may not be expected to 
give close numerical agreement with observation. 
We should, however, require that the model give 
the correct order of magnitude of physical quan- 
tities and their variations. This test we propose 
to apply to the free volume theory of liquids. 
The implications of the theory with respect to 
the mechanical properties of liquids have not 
*This paper is based on work done for the Office of 


Scientific Research and Development under Contract 
OEMsr-262. 


previously been considered, to the author's 
knowledge. 

The important experimental facts regarding 
ultrasonic propagation in liquids are summarized 
below: 


(a) The velocity of sound in a liquid is commonly of 
the order of five to ten times the velocity in the vapor at 
the same temperature. For example, the velocity of sound 
in liquid oxygen! at 90°K is 910 m/sec., while the velocity 
in oxygen vapor? at the same temperature is 180 m/sec. 

(b) The temperature coefficient of sound velocity 
c1(dc/8T)» is of the order of —3 10-3 °C-; in gases the 
temperature coefficient is of the same order but is positive. 

(c) The pressure coefficient of sound velocity c“(dc/dp)r 
is of the order of +5107 cm?/dyne. 

(d) The absorption of sound waves in liquids usually 
exceeds the classical value calculated on the basis of 


1H. W. Liepmann, Helv. Phys. Acta 9, 507 (1936). 
2 A. van Itterbeek and P.Mariéns, Physica 4, 207 (1937). 
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Ea 
A 
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Fic. 1. Intermolecular separation L and free length Ly on Eyring model. 


ordinary viscosity losses; in the case of CSe, for example, 
the observed value is several thousand times greater than 
the classical value. 

(e) The temperature coefficient of sound absorption 
a"(da/dT)» is of the order of +0.01 °C~. Both positive 
and negative signs occur, whereas on classical theory the 
sign is always negative, as the viscosity has a negative 
temperature coefficient. 

(f{) The pressure coefficient of sound absorption 
a'(da/dp)r is —6X10~ cm?/dyne for toluene, the only 
liquid for which a determination has been reported. 


The two traditional concepts of the liquid 
state regard a liquid as either a very dense gas 
or else as a solid in which the long range order 
between the molecules has broken down. Both 
of these concepts account qualitatively for the 
two distinguishing properties of liquids; namely, 
that a liquid combines the small shear modulus 
characteristic of gases with the small compres- 
sibility characteristic of solids. To explain in a 
convenient and natural manner the origin of 
other physical properties of liquids, such as 
fusion, viscosity, dielectric constant, density, 
and solution effects, sometimes the gas and 
sometimes the solid point of view is used. 

In understanding the propagation of sound 
in liquids the gas model is particularly suitable. 
lt turns out that sound velocities can be pre- 
dicted qualitatively from this model; and, 
further, the absorption of sound in liquids can 
be related to the absorption processes in gases, 
which are fairly well understood. On the other 
hand, the absorption processes in solids are 
poorly explained at the present time. 


II. SOUND VELOCITY 


The relationship between the velocity of sound 
in the liquid and gaseous states has been made 
plausible by the work of Eyring and his col- 


laborators,’ who have given a simple and useful 
explanation of the relation in terms of the gas 
model. The relation is derived by Eyring in a 
qualitative way and is not quite correct; later 
on we shall consider a rigorous equation of state 
for this model from which the sound velocity 
can be calculated unambiguously. 

Eyring pictures the sound wave as traveling 
with infinite velocity within a molecule and with 
gas kinetic velocities through the space between 
the molecules of the liquid (Fig. 1). This space 
is termed the ‘‘available volume.’’ The molecules 
effectively short-circuit a part of the path of the 
sound wave. If the intermolecular separation is 
L and the “free length’”’ Ly, the molecules short- 
circuit all but the fraction L;/L of the path. The 
sound velocity in the liquid is then given by 


Ciiq = (L/L) Cons; (1) 
where kinetic theory gives 
Cons” _ Yeas(RT/M), (2) 


with Ygas=C>p/C,, the ratio of specific heats in 
the gas. 

The reasoning leading to Eq. (1) from the 
specified properties of the model is not rigorous. 
It is possible to derive a rigorous equation of 
state for the free volume model, and this has 
been done by Tonks‘ in a paper on the statistical 
mechanics of assemblies of hard elastic spheres. 
Recently, Rice’ and Frank* have drawn attention 
to the to the importance of Tonks’ work in relation to 

~ 9H, Eyri Eyring and J. O. Hirschfelder, J. Phys. Chem. 41, 
249 (i937): |. F. Kincaid and H. Eyring, J. Chem. Phys. 
5, 587 (1937), 6, a (1938); J. O. Hirschfelder, D. P. 
Stevenson, and H. Eyring, J. Chem. Phys. 5, 896 (1937). 

5 Tonks, Phys. Rev. 50, 955 (1936). 

5Q. K. Rice, J. Chem. Phys. 10, 653 (1942); 12, 1 


(1944); 14, 335 (1946). 
6H. S. Frank, J. Chem. Phys. 13, 478 (1945). 
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theories of the liquid state and have considered 
some of the thermodynamic consequences of 
the model. 


The Tonks Equation of State 


The Tonks equation of state for a classical 
three-dimensional gas of hard elastic spheres in 
the limit of close packing is 


pV(1—63) = NT, (3) 


where @ is the packing fraction with respect to 
the closest possible packing. We have @= Vo/V 
=1—(V.,/V), where Vo is the minimum possible 
volume and V4, the available volume, is defined as 
V— Vp. For close packing as in a liquid 


1-9 V./3V. (4) 


In this approximation the equation of state 
becomes 


pbVa=3NkT. (5) 
The form 
p=(3RT/M) pi, (6) 
where 
a Or | 
—=--— (7) 
P1 Pp Po 


is also useful. Here po is the density corresponding 
to the closest possible packing. 

The remarkable and simple form Eg. (5) 
which the equation of state assumes makes it 
easy to deduce various properties of the model. 
This form does not appear to have been stated 
explicitly before. 

Before using Eq. (5) we shall consider its 
physical significance as applied to actual liquids. 
The pressure p is understood to be made up of 
two terms, one representing the effect of the 
attractive intermolecular forces and the other 
representing the external applied pressure. Let 
us imagine the intermolecular potential energy 
to be altered by removing the attractive part of 
the potential while leaving the repulsive part; 
the “internal pressure’? may then be regarded 
as that pressure which would compress the 
system to the volume actually occupied. The 
implied relationship between internal pressure 
and temperature is plausible on the free volume 
model, since the internal pressure does balance 
the kinetic pressure of the molecules regarded as 
moving freely within the volume. It is not ob- 


vious, however, that the thermal expansion at 
constant external pressure should be given cor- 
rectly by the Tonks equation. This amounts to 
the assumption that the balancing pressure due 
to potential energy is independent of the volume. 
The experimental checks of Eqs. (10) and (15) 
below may be taken as justification of this 
assumption within the generous limits of error of 
this paper. 

It is important to note that what is called the 
available volume V, in this paper is not identical 
with the free volume V; as customarily defined, 
although the same model is used. The available 
volume is the difference between the actual 
volume and the minimum possible volume and 
is given by V(1—8@); the free volume is defined*-* 
as the volume of possible motion of the center 
of a single molecule, and is given, per mole, by 
V(1—6')*. For close packing V;=V.(V.?/27V”). 
The available volume is a natural quantity to 
use in the equation of state; the free volume is a 
natural quantity to use in the partition function. 
The numerical values of the free volume given 
by Kincaid and Eyring are defined in terms of an 
approximate equation for sound velocity in 
liquids; the use of the more exact equation 
below will necessitate small changes in the 
numerical values. 

The velocity of sound is given in the usual 
way by c?=(0p/dp)s=(p1/p)?(9/dp1)s, so that 


V C,+3R\'i/3RT\3 
~ (EVE) 
Va C. M 
V 3V1ia\ ? 
(E)Caym 
Va Yeas 


where Y1iq has been set equal to (C,+3R)/C,, in 
accordance with Eq. (16). It may be noted that 


Cliq and Cyas are to be taken at the same tem- 
perature. 








Calculation of Sound Velocity 


It is not practical to give a direct comparison 
of Eq. (8) with measurements of sound velocity 
in liquids since there is no really accurate inde- 
pendent method for determining the available 
or the free volume (see, however, Kincaid and 
Eyring*). We may, however, make an estimate 
of the available volume by considering the 
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TaBLeE I. Comparison of available volume as calculated 
a expansion and sound velocity data for 
300°K. 








Calculated avail- 


Available volume increments able volume in 





cc/mole cc/mole derived 
1 Ve V3 from 
0°K to fusion melting thermal sound 
melting volume point to expansion velocity 
Substance point change 300°K [Eq. (10)] [Eq. (8)] 
CHBr; 8 10 (2) 20 19 
CCl 11 4 5 20 28 
Benzene 10 10.3 20 23 25 














Notes: Thermal expansion and fusion volume change data are from 
the Landolt-Bérnstein tables. Vi is based on a linear extrapolation of 
measured values obtaining between 80°K and 200°K; this overesti- 
mates the change at low temperatures (since by the Nernst theorem 
the expansion coefficient must vanish at 0°K) and underestimates the 
change at high temperatures. No allowance has been made for any 
effective ‘‘zero-point” available volume. 


increase in volume between absolute zero and the 
temperature concerned. 

If we are to take the free volume picture 
seriously there should exist some correlation 
between the molar available volume as calculated 
using Eq. (8) from sound velocity measurements 
at room temperature and the actual increase in 
molar volume of the material between 0°K and 
room temperature. Such a relationship does in 
fact exist. 

We write 

V=VotVa, (9) 


where Vo is the molar volume corresponding to 
the closest possible packing in a geometrical 
sense, and -V, is the molar available volume. Let 


Va=VitVet Vs, (10) 


where 


V,=molar volume change on heating solid from 0°K to 
the melting point; 

V2=molar volume change on fusion; and 

’3=molar volume change on heating liquid from melting 
point to the temperature concerned. 


Unfortunately the expansion data from which 
the V’s can be calculated are inadequate, and a 
certain amount of guesswork enters into any 
estimate of the available volume based on ther- 
mal expansion data. A comparison of available 
volumes as calculated from sound velocity and 
thermal expansion data is given in Table I. The 
agreement is as good as can be expected. 

This evidence shows that the pictorial concept 
of the available volume has a real physical basis. 
It is in principle possible to give an approximate 


calculation of sound velocities from thermal 
expansion data alone. 


Temperature Coefficient of Sound Velocity 


We may determine the available volume from 
the sound velocity at one temperature and using 
this value compare the calculated with the ex- 
perimental value of the temperature coefficient 
of sound velocity c~'(dc/dT),. The velocity of 
sound in liquids has usually a negative tempera- 
ture coefficient; that is, the velocity decreases 
with increasing temperature. This is contrary to 
the position with gases, which show a positive 
temperature coefficient. In terms of the free 
volume picture of a liquid the effect of the 
increased distance the molecules have to travel, 
owing to thermal expansion, dominates the 
increase in molecular velocity. 

We may derive an expression for the tempera- 
ture coefficient of sound velocity on the free 
volume model. From Eqs. (2) and (8), 


c'(8c/8T) p= V-"(8V/8T) »— Va-"(8 Va/8T) » 
+(2T)-'+ (2y1iq) "(OV1ig/9T), (11) 


where the variation in 7gas is neglected. Now AV 
is essentially equal to AV,, since the expansion 
of the molecules themselves may be neglected 
in comparison with the expansion of the available 
volume. Writing A for the temperature coeffi- 
cient of cubic expansion,’ we have 


c-'(8c/8T)» = (2T)-!— A[(V/Ve) —1] 
+(2y1iq)"(Onig/9T)p; (12) 


this equation may be written in terms of ob- 
servables by using Eq. (8): 


Yeas i Cliq } 
-'(ac/@T)p= (27) A} ( ) (=)-1] 
3Y iq Cgas 


+ (2¥1iq)"(OrV1iq/9T)p. (13) 





A detailed comparison of experimental values 
of the temperature coefficient with calculated 
values from Eq. (13) is given in Table II. The 
purpose in giving this and other tables is not in 

7 The conventional symbol for both the cubic expansion 
coefficient and the pressure attenuation coefficient is a; 


in this paper, A is used for the cubic expansion coefficient 
to avoid confusion. 
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TABLE II, Temperature coefficient of sound velocity c~(dc/8T),—comparison of Eq. (13) deduced from the Tonks 
equation with experimental values (data-tor T=300°K). 








Temp. coef. of 
Yliq™ sound vel. °C71 
(Oyliq/OT )p exper. 
? —4.8X 107 
—0.6X 1073 
+1.2 
—0.6 
—0.8 


Cubic exp. 
coef, °C-1 


3 1.16 107% 
4 1.24 
6 1.20 


Cp/Co 
s liq calc. 


—1.9X 10-3 
7 


Substance ga 


(1.2) 1 
1.08 1 
1.24 1 
1.08 1 
(1.2) 1 
111 1 
1 
1 
1 
1 
1 


~ 
= 


Cliq/Cgas 


7.47 
7.00 
5.54 
6.93 
7.42 
6.12 
6.35 
5.40 
8.96 
1 6.92 
117 3.50 





Amy] acetate 

Benzene 

Carbon disulfide 

Carbon tetrachloride 

Chlorobenzene 

Ethyl alcohol 

Heptane 

Methyl iodide “a. 

Nitrobenzene (1. 
(1. 


6 1.23 
5 0.98 
8 1.10 
(1.2) 19 1.24 
(1.3) 1.24 
(1.3) 0.83 
Toluene 4 

Water 0 


3 
4 
5 
A 
3 
| 

1 

3 

3 
3 
0 


NW UN G2 Ga HE On SG 
woOodrK Ar Or OU 


* 
0. 








Notes: Values of Cp/C» in parentheses are assigned arbitrarily; other values from E. Hiedemann Ultraschallforschung (Walter de Gruyter and 
Company, Berlin, 1939). Cubic expansion coefficients are from the Landolt-Bérnstein tables. Temperature variation of Yliq from E. B. Freyer, 
J. C. Hubbard, and D. H. Andrews, J. Am. Chem. Soc. 51, 759 (1929). Experimental values of the temperature coefficient of sound velocity are 
taken from Hiedemann’s book and from the preceding paper by J. R. Pellam and J. K. Galt. The values used for cgas are calculated from Ea. (2). 


order to emphasize specific cases of agreement or 
disagreement, but rather to exhibit the con- 
sequences of the free volume equation of state. 
It seems likely that Eq. (5) will find wide use in 
problems of the liquid state, so that it is well to 
understand some of the limitations and capa- 
bilities of the equation. 

It is seen that the calculated values are in 
order of magnitude agreement with the measure- 
ments, but that the calculated values under- 
estimate the value somewhat. It is significant 
that the order of magnitude of the results derived 
from the simple model of hard elastic spheres is 
correct, and that the existence of an ‘‘anomalous”’ 
positive temperature coefficient for water is cor- 
rectly indicated by the theory. This agreement 
is attributed to the use of experimental values 
of the thermal expansion coefficient, rather than 
the theoretical values given by Eq. (15) below. 


Pressure Coefficient of Sound Velocity 


It is of interest to consider the pressure coef- 
ficient of sound velocity c~1(dc/dp)r. Using the 
Tonks equation we find 


c'(0¢/0p)r =Brl (V/Va) -1] 


Y zas , Cliq 
eel (5) (=) 4} 
3Y liq Cgas 


where Br is the isothermal compressibility and 
is equal to y1iq times the adiabatic compressi- 
bility. The variation of yiiq with pressure should 
probably be considered, but relevant data could 


(14) 


not be found. A comparison of experimental with 
theoretical values is given in Table III, from 
which it appears that the predictions of the 
theory are slightly low. The positive sign of the 
pressure coefficient arises because the sound path 
through the available volume is shortened by 
compressing the volume. 

The fact that both the temperature and pres- 
sure coefficients of velocity as calculated are 
somewhat low suggests that the available volume 
is overestimated by the sound velocity equation. 


III. DISCUSSION OF FREE VOLUME THEORY 


It is appropriate at this point to discuss the 
basis for the application of the Tonks equation 
of state to actual liquids. Five main assumptions 
underlie this application: 


(a) It is assumed that the internal structure of the 
molecule is “‘hard’’—that is, the internal compressibility 
of a molecule is neglected in comparison with the com- 
pressibility of the system as a whole. This assumption is 
justifiable.* 

(b) It is assumed that the molecules are spherical. It 
should be possible to generalize the derivation of the 
equation of state to describe other shapes. 

(c) It is assumed that the compressibility of the mole- 
cules with respect to each other—that is, the inter- 


* The internal compressibility of a molecule Bmoie is of 
the order of L/k, where L is the length of the molecule 
and k is the force constant. The force constant is given 
by k=(2xf)?m, where the frequency f is a fundamental 
vibration frequency of the molecule and is of the order of 
1000 cm or 3X10" c.p.s.; m is the reduced mass of the 
molecule and is of the order of 3X10-* g for a light 
molecule. With these values k~10* dynes/cm, so that 
Bmole~10-" cm?/dyne, which is to be compared with 
Bgas™~10- cm?/dyne under the equivalent pressures ob- 
taining in a liquid. 
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molecular compressibility—is purely kinetic in origin. The 
repulsive forces are imagined to rise sharply at the 
boundary of the hard sphere. This is quite a drastic as- 
sumption. 

(d) It is assumed that the internal pressure arising from 
the attractive intermolecular forces is approximately 
independent of volume. This is also a drastic assumption. 

(e) The Tonks derivation is based on classical me- 
chanics, and no account is given of quantum effects. It is 
almost certain that quantum effects are of considerable 
importance in liquids, since the spacing between the energy 
levels of a molecule confined within its own free volume 
is of the order of kT for room temperature. 


In the light of the assumptions which have 
been made it is of interest to consider further 
consequences of the Tonks equation. We have 
already seen that the equation leads to values of 
the temperature and pressure coefficients of 
sound velocity which are correct to within an 
order of magnitude at least. 


Thermal Expansion Coefficient 


The measure of agreement with experiment 
attained by the calculated values of the tem- 
perature coefficient of velocity may be attributed 
in part to the use of experimental values of the 
cubic expansion coefficient when substituting 
in Eq. (13) to obtain the results displayed in 
Table II; this is particularly true with respect 
to water. We may, however, calculate the cubic 
expansion coefficient from the Tonks equation of 
state, obtaining 


A=V-\(AV/8T)»=(V2/V)(1/T). (15) 


TABLE III. Pressure coefficient of sound velocity 
¢(dc/8p)r—comparison of calculated values from Eq. 
(14) with experimental values. 








Adiabatic 
compressibility 
cm?2/dyne 


59.0 10-” 


Pressure coefficient of 
sound velocity cm?/dyne 
experimental calculated 


3.7X10- =1.7 107° 


Substance 


Carbon disulfide 
Carbon tetra- 
chloride 74.5 
Ether 138.7 
Ethyl bromide 88.1 
Methyl acetate _ 80.9 
enzene 64.9 
Toluene 68.2 





NO DY dD KY W bo 
Snore 








Notes: Values of the adiabatic compressibility are from Hiedemann’s 
k. Experimental values of the pressure coefficient of sound velocity 
are from J. C. Swanson, J. Chem. Phys. 2, 689 (1934), except for the 
value for methyl acetate, which is from P. Biquard, Comptes rendus 
(Paris) 206, 897 (1938). 


®E. Bartholomé and A. Eucken, Trans. Faraday Soc. 33, 
45 (1937). 
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TABLE IV. Cubic expansion coefficient—comparison of 
calculated values from Eq. (15) with experimental values 
—V-(8V/8T), in °C. 








Cubic expansion coefficient 
experimental calculated 


1.24X10-* 0.95x 107 
1.20 1.19 
1.23 0.98 
0.98 0.84 
1.10 0.98 
0.18 0.52 
1.24 1.07 
0.12 1.45 
0.97 (1.1) 


= 
» 
= 


Substance 





Benzene 

Carbon disulfide 
Carbon tetrachloride 
Chlorobenzene 
Ethyl alcohol 
Mercury 

Methyl iodide 
Water 

Water (400°K) 


NN WAwW FE ww 
QW Of OO 


— 
~~ 











Notes: The data are for T =300°K. The values of the ratio of the 
volume to the free volume (V/Va) are from Table II. 


A comparison of experimental with theoretical 
values is given in Table IV. The agreement is 
excellent for the organic liquids, but is poor for 
water and mercury, as might be expected. We 
have neglected entirely the role of the metallic 
electrons, so that the present discussion is often 
inapplicable to liquid metals. Water is anomalous 
in general behavior at room temperature and 
below, but near the boiling point the calculated 
value of the cubic expansion coefficient is in 
good agreement with experiment. 


Difference of Specific Heats (C,—C,) 


Using Eq. (5) we may compute the difference 
(C,—Cy,) for the Tonks equation, finding 


C,—C,=p(dV2/8T)p=3R (16) 


instead of simply R which obtains for an ideal 
gas. A comparison is given in Table V of experi- 
mental with theoretical values of the difference 
(C,—C,). The agreement in the case of the 
liquefied gases is good except for the case of 
hydrogen, where quantum effects are of im- 
portance. For the organic liquids the calculated 
(Cp— Cy») of 6 cal./mole/°C is generally low. With 
regard to C,, one might expect on the basis of 
the free volume model that C, would be equal in 
the liquid and gaseous states at the same tem- 
perature. It is well known,’ however, that C, is 
greater in the liquid state. 


IV. SOUND ABSORPTION 


The two principal mechanisms responsible for 
the absorption of sound in liquids are the ordinary 
shearing viscosity and the energy loss associated 
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with incomplete excitation of vibrational energy 
states of the molecules. In many liquids the 
latter mechanism is the dominant cause of 
attenuation. A detailed bibliography of the 
literature concerned with incomplete vibrational 
excitation in gases has been given by Richards.° 
The theory developed for gases was first applied 
to liquids by Mandelstam and Leontovic,'® and 
later extended by Kneser." Herzfeld has given 
an interesting discussion of absorption in liquids, 
while Tisza has considered the important im- 
plications of the subject with respect to the 
hydrodynamical equations. 

The program of this section is to give a 
physical order-of-magnitude estimate of the 
maximum absorption expected from relaxation 
processes; to go through a detailed calculation 
of absorption for the Tonks equation of state; 
to consider the temperature and pressure coef- 
ficients of absorption ; and, finally, to discuss the 
relation between sound absorption and molecular 
structure. 

It is possible to make a simple estimate of 
the maximum absorption caused by relaxation 
of a vibrational mode. The average value of 
the energy in a sound wave is of the order 
of ApPAV. Now Ap/p=—AV/V+AT/T for an 
ideal gas, where the first term on the right 
represents “‘mechanical”’ effects and the second 
term “thermal” effects of the pressure. For an 
adiabatic process the ratio of the thermal term 
to the mechanical term is (y—1), where y is the 
ratio of the specific heats. Let C; be the specific 
heat associated with the: relaxing degree of 
freedom; then C;/C, represents the fraction of 
thermal energy in the mode concerned, and 
(C;/C.)(y—1)/y is an approximate upper limit 
to the fraction of the total energy of the sound 
wave dissipated in one cycle. We should of 
course find absorptions of this order of magnitude 
only in the neighborhood of the relaxation fre- 
quency characterizing the exchange of energy 
between translational and vibrational degrees of 
freedom. 

Let us compare this estimate of the peak ab- 


®W. T. Richards, Rev. Mod. Phys. 11, 36 (1939). 

10, Mandelstam and M. Leontovic, Comptes rendus 
(Doklady) Acad. Sci. U.S.S.R. 3, 111 (1936). 

11H. O. Kneser, Ann. d. Physik 32, 277 (1938). 

2K, F. Herzfeld, J. Acous. Soc. Am. 13, 33 (1941). 

%L. Tisza, Phys. Rev. 61, 531 (1942). 


TABLE V. Comparison of experimental values of (C,—C,) 
with calculated values from Eq. (16). 











(Cp —Cr) 

Cp Cp/Crv cal./mole/°C 

Substance expt. expt. expt. calc. 

Argon (87°K) 10.0 pa a 6 
Oxygen (90°K) 13.0 7 5.3 «66 
Nitrogen (76°K) 13.6 2.0 66 6 
Hydrogen (20°K) 6.9 Be as 
Benzene 31.6 1.44 9.7 6 
Carbon tetrachloride 30.7 1.46 9.7 6 
Ether 40.6 1.33 10.0 6 
Ethyl alcohol 27. 1.18 42 6 
n-heptane 49, 1.19 7.8 6 
n-propyl! alcohol 35.6 1.16 5.0 6 
Toluene 37.8 1.34 10.0 6 











Notes: The data are for 300°K unless otherwise specified. Experi- 
mental values of Cp/C» are from L. Bergmann, Der Ultraschall (Edwards 
Brothers reprint, Ann Arbor, 1944), 


sorption with some experimental values. For CO, 
gas we have (y—1)/y~0.23, and C;/C,~0.3, so 
that about 0.07 of the energy is dissipated per 
cycle. This means that the intensity is down 1/e 
in about 14 wave-lengths; using data for high 
frequencies we should estimate 8 wave-lengths. 
Measurements by Fricke“ show that at the re- 
laxation frequency of 20 kc/sec. the intensity is 
down by 1/e in about 4 wave-lengths, indicating 
that our order of magnitude estimate is generally 
correct. The famous example of a highly absorb- 
ing low viscosity liquid is CSe, for which the 
measured absorption" at 1 mc/sec. is about 2000 
times the classical value as calculated from vis- 
cosity and heat conduction losses alone; yet at 
this frequency the intensity is only down by 1/e 
in 50 wave-lengths. At 7 mc/sec. the intensity’® 
is down by 1/e in about 12 wave-lengths; this 
value appears quite reasonable in the light of the 
above estimates. 


Absorption in Liquids on Free Volume Model 


It is instructive to calculate the absorption in 
detail for the free volume model of the liquid 
state, applying the method of Herzfeld and Rice"’ 
to the Tonks equation of state, Eq. (5). In a 
sense the calculation is but a formal extension 
of the treatment of absorption in gases developed 
by Herzfeld and Rice, and others. Let C; be the 
specific heat of the relaxing mode of vibration, 

“4 E. F. Fricke, J. Acous. Soc. Am. 12, 245 (1940). 

15 J. Claeys, J. Errera, and H. Sack, Trans. Faraday Soc. 
33, 136 (1937). 

16S, Parthasarathy, Current Sci. 6, 501 (1938). 


17K, F, Herzfeld and F. O. Rice, Phys. Rev. 31, 69! 
(1928). 
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and C, the specific heat of all the other degrees of 
freedom of the molecule. Let 7’ denote the 
effective temperature of the energy associated 
with C;, and T the temperature associated with 
C,. Losses owing to ordinary viscosity and heat 
conduction are not considered; it can be shown 
that these are additive to the first approximation. 
From the Tonks equation of state, 


p=3RTpi/M; (1/p1)=(1/p)—(1/p0). (17) 
The force equation gives the particle acceleration 
du/dt= —(1/p)(dp/dx) 


= —(3RT/M)(pi/p)L(0 In T/dx) 
+(0 In pi/dx)], (18) 


while the conservation equation is 

0 In p/dt= —du/dx =(p/pi)(0 In p;/dt). (19) 
From the first law of thermodynamics 
C(dT /dt) + CT" /dt) =(pM/p:*)(dp1/dt). (20) 


The relaxation time constant 7’ is introduced* 
through the equation 


C,(0T’/dt) =C(T—T’)/r’. (21) 


We consider small changes in p;, T and 7”, 
and suppose that the small changes are of the 
form exp [j(wt—kx) ], where k =k: — jke. On sub- 
stituting in Eqs. (17) through (21) we find an 
expression for k?. For the technique of the cal- 
culation Herzfeld and Rice many be consulted. 
The imaginary part of k? is —2k,ks, where 


2kiks = w(M 3RT) (p pi)” 
x 0 


WT oo 
; -}. (22) 
1+(wr)? y°y* 





Here 
t=7'(C:/C,") (y?/y"); v°=C>p°/C,°; 
ans = C"/C,*. 


The low frequency specific heat C,° is given by 
C:+C,., while the high frequency specific heat 
C,” is simply C,. 

Now ki&w/ciiqg, where the sound velocity C1iq 
is given by Eq. (8). Also, ke=a, where a is the 
ordinary pressure attenuation coefficient in the 
expression p= po exp [—ax_] for the damping of 

* For simplicity all of the modes of internal motion are 


supposed to relax at the same frequency, in which case C; 
includes the whole of the vibrational specific heat. 


a plane wave. We have ae 
fo 
Cy* C.° fP+fo? 
where f=w/2m and fy=1/2mr7r; here C,—C,=3R 
from Eq. (16). 
The maximum absorption per wave-length 
occurs when f= f» and is given by 
T 4 Cp— Ce Ce 
C,* Cc," 
which is of the same form as the previous esti- 


mate reached by a qualitative argument. 
At low frequencies (f<fo) we have 








a= (rf? feu) — 


(Chiq /f) mex = 





(24) 


C,° 


so that at low frequencies a/f? is constant for the 
attenuation caused by the relaxation mechanism, 
just as for viscosity and heat conduction. These 
results are in agreement with those of Kneser 
and Herzfeld, when slight differences in de- 
finitions are taken into account. 





Co- 
a= f?(r/fociia : (25) 


Temperature Coefficient of Sound Attenuation 


Experimental determinations of the tempera- 
ture coefficient of attenuation in different liquids 
show both positive and negative values. Where 
relaxation is a minor factor, a negative coef- 
ficient is found owing to the decrease in vis- 
cosity with increasing temperature. Where 
relaxation is the predominant factor, a positive 
coefficient is expected when the frequency is 
above the relaxation frequency, and a negative 
coefficient when the frequency is below the 
relaxation frequency. This is because a rise in 
temperature increases the probability of energy 
transfer and thus increases the relaxation fre- 
quency. If f>fo, an increase in fo brings the 
region of maximum absorption nearer to the 
measurement frequency; the opposite holds if 
f<fo. The actual situation, however, may be 
more complicated. 

Considering only relaxation losses we find 
from Eq. (23): 


a~'(da/dT) ,=+ fo-(9fo/dT) > 
+C-(8C;/8T) »— 2C.-"(8C»/8T) » 
—Ciiq '(OCiiqg/9T) p, (26) 
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TABLE VI. Temperature coefficient of sound absorption—interpretation of results in terms of Eq. (26)— 


frequency 15 mc/sec. 














Temp. coef, 
a1(da/dT )p exper + 
Substance “—s classical Tentative interpretation 

Ethyl bromide +0.010 7 rere factor is change in fo, the relaxation frequency. fy<15 
mc/sec. 

Methyl iodide -+ 0.0087 30 Fundamental vibrational frequency of CH;I molecule is fairly low 
(534 cm), so that specific heat change is not important; fo is 
dominant factor. fo<15 mc/sec. 

Methy! alcohol —0.011 3 Change in viscosity contributes only —0.004; velocity contributes 
-+0.003; specific heat term is positive, so that at least —0.010 must 
be accounted for by fo change. fo>15 mc/sec., in agreement with 
Herzfeld,* who concludes on other grounds that fyo~3 X 10* mc/sec. 

n-propyl alcohol — 0.0086 2 As for methyl alcohol, fo>>15 mc/sec. 

n-amy] alcohol —0.014 1S Changes in both viscosity and fo are probably effective here. The 
absorption is nearly classical. fo>>15 mc/sec. 

Benzene +0.011 90. The evidence of the temperature coefficient of absorption is that 


fo<15 mc/sec., yet Herzfeld indicates that fo~2000 mc/sec. It is 
possible that these figures represent two different relaxation modes, 
the mode with a weak absorption peak near 10 mc/sec. having a 
predominant effect on the temperature coefficient. Such a peak is 
suggested by a plot of @ vs. f based on collected results. Kneser® 
gives fo>55 mc/sec. 











* See reference 12 in text. 
b See reference 11 in text. 


where the plus sign before the term fo'(0fo/dT) p 
applies to the case f>fo, and the minus sign 
applies when f<fo. The intermediate case is 
more complicated. Consider now the magnitude 
of the term arising from fo: Landau and Teller'® 
show that under certain assumptions 


fo <exp [—3(wav)?/8(M/2kT)"8], (27) 


where a@ is the atomic radius, M the mass of the 
molecule, and v the frequency of the vibrational 
mode. The exponential factor gives essentially 
the probability of exchange of vibrational and 
translational energy in a collision, and comes out 
to be of the order of 10~* to 10~®, in agreement 
with the experimental results for gases, where a 
collision frequency of the order of 10!° collisions 
per second is associated with relaxation fre- 
quencies of the order of 10° cycles per second. 
From Eq. (27) we have 


fo (0fo/dT) p=q/3T, (28) 


where g is the argument of the exponential factor 
in the above equation. Since g is of the order of 
10, we have fo (0fo/dT) p>~ 10. 

Let us consider the order of magnitude of the 
other terms in Eq. (26). The term arising from 
Ciiq is of the order of 3X10-* and may usually 
be neglected. The two terms in the specific heat 
may be combined in One, since AC,AC;, giving 


18. Landau and E. Teller, Physik. Zeits. Sowjetunion 
10, 34 (1936); cf. C. Zener, Phys. Rev. 38, 277 (1931). 


(0C;/dT) »(C;"—2C,“). Since the specific heat 
C; usually corresponds to a vibration, it is given 
by an Einstein function. The relative slope of an 
Einstein function varies from very steep at low 
temperatures to flat at high temperatures. For 
molecules containing at least two non-hydrogen 
atoms we are up near the flat part of the curve 
at room temperature, and the contribution of 
the change in specific heat to a(da/dT), will 
be of the order of 5X10-* or less. It therefore 
appears that the term arising from fp often may 
be the most important single contribution; in 
some cases, however, the term in the specific 
heat may be significant. 

In liquids for which the viscosity 9 is a major 
factor contributing to the attenuation a term 
arising from 7 of the order of 10-? must be con- 
sidered. Including viscosity, the attenuation a 
as given by Eq. (23) is increased by the additive 
Stokes term 82? f?n/3c*p. 

Measurements of the temperature coefficient 
of absorption have been reported by Pellam and 
Galt,!® Bazulin,?® Sérenson, Baumgardt,” and 

19 J. R. Pellam and J. K. Galt, “Ultrasonic propagation 
in liquids. I. Application of pulse technique to velocity and 
absorption measurements at 15 megacycles,” J. Chem. 
Phys. 14, 608 (1946). 

20P, Bazulin, Comptes rendus (Doklady) Acad. Sci. 
U.S.S.R. 14, 273 (1937). 

21 C, Sérenson, Ann. d. Physik 27, 70 (1936). 


%E. Baumgardt, Comptes rendus (Paris) 202, 203 
(1936). 
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Hunter.% An interpretation of the results of 
Pellam and Galt is given in Table VI. 

Study of the table suggests that the tem- 
perature coefficient of sound absorption is a 
useful quantity in the study of molecular re- 
laxation times. The indications are not always 
unambiguous, but in the absence of more com- 
plete experimental knowledge of the acoustic 
absorption spectrum of liquids the temperature 
coefficient provides pertinent information. 


Pressure Coefficient of Sound Attenuation 


The pressure coefficient of sound attenuation 
in toluene has been determined by Biquard™ who 
measured the absorption at 7 mc/sec. between 
pressures of 1 to 500 atmos., giving a—'(da/dp)r 
= —6X10-'° cm?/dyne. At one atmosphere a/f? 
=80X10-"7 sec.2 cm, as compared with the 
classical value 8X10-'’, suggesting that the 
principal part of the absorption may be ascribed 
to relaxation processes. On the assumption that 
f<fo and that changes in specific heat may be 
neglected, we have from Eq. (25) 


a '(da/dp) r= —fr(Ofo/dp)r—c(dc/dp)r. (29) 


Now fo at constant temperature is proportional 
to the collision frequency, so that 


— fo (0fo/dp)r =L;“"(OL;s/dp) r 
=—(V/3V.)Br, (30) 


where L; is the free length as in Fig. 1. The mag- 
nitude of the term in fo is 1X10~'°, while the 
term in the sound velocity c is 3X10~°. The cal- 
culated value of a '(da/dp)r comes out to be 
—4X10-'° cm?/dyne, in fair agreement with the 
experimental value —6X10-!° cm?/dyne. 

The physical interpretation of the result is that 
the molecules are brought closer together by the 
application of pressure, thus increasing the col- 
lision and relaxation frequencies. If f<jfo an 
increase in fy moves the region of maximum ab- 
sorption away from the working frequency. Also, 
the wave-length is increased by pressure, so that 
the absorption per unit length will decrease even 
if the absorption per wave-length is constant. 


* J. L. Hunter, J. Acous. Soc. Am. 13, 36 (1941). 
* P. Biquard, Rev. d’Acoustique 8, 130 (1939). 


Sound Absorption and Molecular Structure 


According to present ideas the magnitude of 
the non-classical part of the attenuation is deter- 
mined by the position and magnitude of the 
relaxation frequencies of the molecule. It is dif- 
ficult to calculate these frequencies since they 
are dependent upon the detailed coupling 
between different modes of motion. In gases it is 
well-confirmed experimentally that the relaxa- 
tion mechanism accounts for most of the non- 
classical absorption. Relaxation frequencies in 
pure gases are of the order of 10'—10° c.p.s.; 
since ‘‘collisions’’ in liquids are at least 10* more 
frequent than in gases we would expect relaxation 
frequencies in liquids to be of the order of 10’ 
c.p.s. The indication of the present evidence is 
that the actual frequencies are of this order or 
higher. The existence of well-defined vibrational 
states in liquid molecules is known from infra-red 
and Raman spectra evidence. 

A review of the absorption of sound in various 
liquids at 15 mc/sec. shows several regularities 
with regard to molecular structure. 

(a) The three most absorbing liquids of low 
viscosity are non-polar—CSe, CCl4, benzene. 


Substance Gexper/ classical 
CS: (7 mc/sec.) ~600 
CCl, 25 
Benzene 90 


Other non-polar liquids are highly absorbing— 
heptane, hexane, toluene. This tendency suggests 
that the coupling of translational and vibrational 
motion is small in molecules with small electro- 
static interactions, giving a long lifetime to the 
vibrational state and a low relaxation frequency. 

(b) Highly associated liquids tend to show 
low, near-classical absorption at 15 mc/sec. The 
alcohols and water are good examples of asso- 
ciated liquids, and they show the following 
values: 


Substance 


Methyl alcohol 
Ethyl alcohol 
Propyl alcohol 
Butyl alcohol 
Amy] alcohol 
Water 


Gexper /Aclassical 


This tendency suggests that the formation of 
molecular groups is conducive to the transfer of 
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vibrational energy and therefore leads to high 
relaxation frequencies. This is expected. In the 
limiting case of a solid there is little distinction 
between translational energy of molecules with 
respect to one another and vibrational energy 
internal to a molecule. In solids we have relaxa- 
tion times of the order of 10-™” sec., as estimated 
from heat conductivity data. 


CONCLUSIONS 


It is found that the Tonks equation of state as 
applied to liquids is a suitable framework for 
making close order-of-magnitude estimates con- 
cerning quantities involved in sound propagation. 
The detailed numerical agreement is usually un- 
satisfactory, thus indicating the limitations of 
the model. In particular the temperature and 
pressure dependence of the compressibility is 


CHARLES KITTEL 


given only roughly by this model. On the other 
hand, the thermal expansion coefficient is given 
quite accurately; also, it is possible to attach 
physical significance to the concept of the 
available volume. 
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A theoretical evaluation of some of the normal frequencies of vibration of the isomeric 
octanes has been made. This report gives the skeletal frequencies for those of the eighteen 
isomeric octanes which lead to representations having 7 or fewer frequencies. The forces are 
assumed to be valence forces. The force constants used are 4.12105 dynes per cm for the 
C—C bond stretching, and 0.37 X 10° dynes per cm for the bending between two C—C bonds, 
the values being based on computations by the author for n-propane, using the experimental 
data of E. F. Barker and Violet L. Wu. The angle between two C—C bonds was taken as 
107° whenever possible. The results of this investigation are summarized by means of charts, 
showing the location of bands calculated. For the molecules reported all normal skeletal 
frequencies for a valence force model were evaluated, and attempts at correlation between 
experimental and theoretical bands have been made. 


























imation, and the results obtained with it cannot 
be expected to be identical with experimental 
values. It is found, for example, that in many cases 
zero frequencies arise, caused by the fact that, 
for the force field chosen, the displacements of 
the atoms from the equilibrium positions give 
rise only to terms of third or higher order in the 
potential function. If additional forces, such as 
those restricting free rotation about bonds, were 
introduced these zero frequencies might be re- 
moved. The actual frequencies, though not ex- 
actly zero, should, if the valence force model is 
any good at all, be very small. The frequencies 
obtained are to be regarded as approximations, 
the force model being too crude to obtain exact 


INTRODUCTION 





T has been shown? that for small oscillations 
which are assumed to be simple harmonic the 
frequencies obtained on the basis of classical 
mechanics are identical with those obtained by 
quantum mechanics. For this reason, it is still 
customary to use classical theory of small oscil- 
lations to obtain the normal frequencies of vibra- 
tion of molecules. 

In obtaining the frequencies belonging to 
normal modes of vibration, it is assumed that the 
potential energy is a function of the displacement 
coordinates, and, in order to express this function 
explicitly, it is necessary to make use of a definite 
force model. The earliest model used was the 
central force model, in which the potential energy 
depends only on the changes in the distances 
which separate the atoms of the molecule. It was 
found, however, that such a model is not ade- 
quate to represent the force field, as far as the 
study of vibration spectra is concerned. In the 
valence force model the atoms are connected by 
valence bonds, and it is assumed that a change 
in angle between two bonds, as weil as a change 
in the length of a bond, results in an increase in 
potential energy. The valence force model has 
been more successful in the treatment of vibra- 
tion spectra, though it is still, at best, an approx- 












results. 

In the following discussion, a molecule is 
assumed to consist of m atoms bound to equi- 
librium positions by valence forces. The atoms 
are bound together in such a way that a change 
in length of a valence bond by an amount AL 
gives rise to an increase in potential energy of 
3Kcec|AL|?, and a change A@ in the angle between 
two valence bonds meeting at a single atom 
results in an increase in potential energy of 
4K,(A0)?. Hence, the total increase in potential 
energy, because of an arbitrary change, is 


V=4 > Kec|AL|?+43 > Ko(A8)?, 













' This research was completed at the Harrison M. 
Randall Laboratory of Physics, University of Michigan, 
Ann Arbor, Michigan, 1945. 

2D. M. Dennison, “Infrared spectra of polyatomic 
molecules,” Rev. Mod. Phys. 3, No. 2 (April, 1931). 


where the summation extends over all valence 
bonds in the first term and over all valence 
angles in the second term. 
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The isomeric octanes constitute a group of 
saturated hydrocarbons, each possessing 8 carbon 
and 18 hydrogen atoms. The direct application 
of the above theory to a molecule of that type 


would lead to a secular determinant of order 72, 


and, even with the aid of group theory, a com- 
plete determination of the frequencies would 
present a task too formidable to be attempted. 
It has been found, however, that the skeletal 
frequencies obtained by considering the hydrogen 
atoms frozen to the carbon atoms are not much 
changed. In the present discussion, all hydrogen 
atoms are considered frozen to the carbon atoms, 
thus reducing the original system in each case to 
a system of 8 mass points, made up of the com- 
bined masses of 1 carbon and 3 or fewer hydrogen 
atoms. The present paper gives the skeletal fre- 
quencies of those isomeric octanes which belong 
to point groups of higher symmetry than the C,. 
For the C, group, the frequencies belonging to 
the lower ordered representations have also been 
worked out, but will not be reported here. 

In general, the atoms constituting a hydro- 
carbon tend to stretch out and form a rather 
symmetric pattern, as has been indicated by 
measurement on oil films. The molecules are 
assumed to consist of zigzag chains arranged as 
symmetrically as possible, and since the effect 
of the C—H bonds is eliminated, the distinction 
between optical isomers is ignored. No systematic 
experimental spectra are available which distin- 
guish the optical isomers, and it is felt that 
significant information is obtained about the 
skeletal frequencies with the forms chosen. The 
angles in the main chains of the molecules are 
taken as 107° in the equilibrium position, though 
in some cases this had to be altered because of 
symmetry considerations. Values for the angles 
between valence bonds ranging from 99° to 116° 
have been used by other investigators for various 
molecules involving CH; groups.’ Although the 


3K. W. F. Kohlrausch, Der Smekal-Raman Effekt 
(Verlagsbuchhandlung, Julius Springer, Berlin, 1938). 
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Fic. 1. Skeletal form of n-octane. 81°= 82° = 83° =84°=8s° =86° = 107°. 





Gi gi 
A4 Re 
B& 
5 7 
CHy CH, 





angle would vary with position, 107° was chosen, 
somewhat arbitrarily, as representative, this 
value having been used with some success as the 
tetrahedral angle for certain hydrocarbons. Four 
types of masses occur, representing the masses 
of C, CH, CHe, and CHs; groups, respectively. 
The values for these masses were determined 
from the masses of the hydrogen and carbon 
atoms given by R. T. Birge.* The values used are 


C=19.93X10- g, 
CH =(19.93+1.673) X10 g = 21.6 X10- g, 
CH2= (19.93 +2 1.673) X10- g 
= 23.310 g, 
CH3=(19.93+3 X 1.673) X10 g 
= 24.9 10-* g. 


The force constants used are based on computa- 
tions by the author on m-propane, using the 
experimental data of E. F. Barker and Violet L. 
Wu.® The values chosen as best suited for the 
octanes are 


For C—C stretching, Kcc=4.12 X 10° dynes/cm. 
ia 











For C bending, Ke= .37X105 dynes/cm. 
C 
TABLE I. Table of characters. 
Infra- 

Car E C2 oh i Raman _ red No. 
TO ~AG 1 1 1 ae IA 7 
TO ~Be 1 —-—1 1 i DP fA 2 
TO~A, 1 1-1 -1 F A 3 
To ~B, 1 —1 1 —1 F A 6 
OR 0 T 0 7 
Ur 8 0 8 0 
+1+2 cos dr 3-1 1 —3 
x(R) 24 0 8 0 
= 18 2 8 0 








‘Raymond T. Birge, “A new table of values of the 
CTT physical constants,” Rev. Mod. Phys. 13, 233 
1941). 
os. F. Barker and Violet L. Wu, ‘The infra-red absorp- 
tion spectrum of propane,” J. Chem. Phys. 9, No. 7, 490 
(July, 1941). 
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NORMAL FREQUENCIES OF OCTANES 


MOLECULES BELONGING TO THE 
POINT GROUP C2, 


N-Octane® (Table I, Fig. 1) 


Representation [™: The values obtained for 
the 7 bands belonging to this representation are 


1 1 
—=200.0cem—, —=955.7 cm™, 


Ay Ns 


1 1 
—=224.3 cm, —=1053 cm", 
Ao As 


1 1 
—=443.4cm, —=1054cm~), 
X3 Xe 


1 
—=1127 cm—. 
Az 


Representation I'®: Because of the choice of 
potential function, the values obtained for the 
2 bands belonging to this representation are 

1 


—=0cm™', —=0cm". 
Ai Ae 


Representation I: Because of the choice of 
potential function, the values obtained for the 
3 bands belonging to this representation are 


Representation I: The values obtained for 
the 6 bands belonging to this representation 


TABLE II. Table of characters. 











+1+2 cos or 
x(R) 2 
= 1 








5 Details of the calculation for this, as well as for the 
other molecules are given in a thesis submitted by the 
author to the Horace H. Rackham School of Graduate 
Studies, University of Michigan, 1945. 


Fic. 2. Skeletal form of 2-5-dimethylhexane. 6:°=82° 
= B3°=f,°=Bs° = Be° = Br°=Bs° = 107°. 


1 
— = 86.16 cm—, 


At 


1 
— = 344.0 cm, 


Ae 


1 
— = 396.6 cm™', 


X3 


1 
—= 986.3 cm—!, 


Aa 


1 
— =996.7 cm—', 


As 


1 
—=1104cm—. 


Ae 


2-5-Dimethylhexane (Table II, Fig. 2) 


Representation T®; 


The values obtairied for 


the 6 bands belonging to this representation are 


1 
—=(0cm—, 
Ai 


1 


—=217.9cm", 


Ae 


1 
— = 367.1 cm™, 
3 


Representation [: 


1 
—=851.3 cm™, 


4 


1 
—= 1032 cm, 
As 


1 
—=1146cm™. 
Xe 


The values obtained for 


the 3 bands belonging to this representation are 


1 1 
—=0cm"', — 
1 r 


Representation T®:; 


= 376.0 cm-!, 


1 
—=1113 cm. 


As 


The values obtained for 


the 4 bands belonging to this representation are 


1 
—=0cm", 
Ai 


1 
-—=(0cm"', 
No 


- 


Representation [; 


1 
—=257.9 cm, 
Na 


1 
—=1061 cm“. 
As 


The values obtained for 


the 5 bands belonging to this representation 
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3 bands belonging to this representation are 
1 1 
—=0cm™, —=0cm", 
Ai Ae 


Representation I: The values obtained for 


Fic. 3. Skeletal form of 3-4-dimethylhexane. 61° = 2° 
= B3°=B,°= 107°. 85°=B6° =B7° =Bs° = 126.5° 


1 
=152.6cm—!, —=883.7 cm—!, 


4 
1 
=372.6cm~!, —=1124cm—!. 
2 As 
1 


—= 426.8 cm, 
3 


3-4-Dimethylhexane (Table III, Fig. 3) 


Representation [!®: The values obtained for 
the 7 bands belonging to this representation are 


1 1 
—=222.2cm—!, —=670.9cm~, 


Ai 4 


1 
—=1005 cm—, 


1 
—=276.4cm—!, 
Ae As 


1 1 
—=429.7 cm—!, —=1274cm"—', 
A3 Ae 


1 
—= 1369 cm. 


Az 


Representation [': Because of the choice of 
potential function, the values obtained for the 
2 bands belonging to this representation are 


1 
—=0cm", —=0cm"!. 
1 Ao 


Representation [': Because of the choice of 
potential function, the values obtained for the 


the 6 bands belonging to this representation are 


1 1 
—=146.lcm™, —=856.5 cm, 
Ay Ng 


1 1 
—=255.9cm—, —=1031 cm", 
r Xs 

1 1 
—=554.9cm—!, —=1208 cm. 
A3 Ae . 


MOLECULES BELONGING TO THE 
POINT GROUP C:, 


4-Methylheptane (Table IV, Fig. 4) 


Representation ['): The values obtained for 
the 7 bands belonging to this representation are 


1 1 
—=109.6cm—!, —=743.6cm~, 
A1 N4 


1 1 
—=239.2cm™!, —=984.2 cm, 
Ae 5 

1 1 
—=410.5cm—, —=1089 cm, 
A3 Xe 


1 
—=1288 cm“. 
Az 


Representation I: The values obtained for 
the 6 bands belonging to this representation 


TABLE III. Table of characters. 








Infra- 


Cor a 2 , i Raman red 


TY LAG it #£ IA 
T2~Bg _ 1 DP IA 
TOYA, 1 1 F A 
ro ~B,, 1 F A 





OR 

Ur 

+1+2 cos OR 
x(R) 


— 
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od fc or 


m are 
1 3 5 
CH; cH, 


Fic. 4. Skeletal form of 4-methylheptane. Bi° = B2° = B;° = 64° =B;° = 107°. 
Be° =Br° = 126.5°. 


2-3-4-Trimethylpentane (Table V, Fig. 5) 


1 1 Representation ['™: The values obtained for 


——— Scm™, —— em™, the 6 bands belonging to this representation are 
4 


1 1 
1 aa oi ee oan 
=288.5cm—!, —=1034cm—', “Taian :. taeda 1 


As 


1 1 
—=251.3cm—, —=1027 cm™, 


Ae As 


1 
—=1280 cm. 


Xe 


Representation [®: Because of the choice of —=427.3cm—, —=1288 cm—. 
potential function, the values obtained for the As As 


2 bands belonging to this representation are , , 
_ . Representation ['®: The values obtained for 


1 1 the 5 bands belonging to this representation are 
—=0cm", —=0cm". 
Ai Ae 


* 29313 cm, * 2889.3 cm, 
Representation I: Because of the choice of Ai A 

potential function, the values obtained for the 

3 bands belonging to this representation are 


1 1 
—=374.3cm—, —=1238 cm—. 
Ao As 
1 1 1 
—=0cm", —=0cm", —=0cm. 
1 for Ai Ae 3 — = 587.9 cm—!, 


ition As 


TABLE IV. Table of characters. ° : 
Representation [: The values obtained for 

ite, the 3 bands belonging to this representation are 
gy Raman _ red 


rO~A ) #2 1 1 

rOwB, on 1 DP A —=0cm"', —=1111 cm. 
PO~A» ; —-1 DP IA Ai As 

rO~B, -1 DP A 











1 
Ae 


° 0 = 286.9 cm, 
R 8 
+1+2 cos dr 1 
) s 
8 


Representation I: The values obtained for 
the 4 bands belonging to this representation 
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Fic. 5. Skeletal form of 2-3-4- trimethylpentane. “ge “tal 


= B3° = By° = Bs° = Be° = 87° = 107°. Bs°= 89° = 126 
are 

1 

—=0cm"), =457.2 cm—!, 
Ai As 

1 1 

—=0cm™", —=1115 cm". 
Ae vi 


2-Methyl-3-Ethylpentane (Table VI, Fig. 6) 


Representation T!™: The values obtained for 
the 7 bands belonging to this representation are 


1 
—=175.8cm—, 
A Ag 


= 584.1 cm—!, 


1 1 
—=311.2cm—!, —=1017 cm, 


Ae As 
1 1 
—=410.2cm™, —=1043 cm™, 
As Xe 
1 
= 1364 cm™!. 
hz 


Representation ['@: The values obtained for 
the 5 bands belonging to this representation are 


1 

—=0cm—, —=957.1 cm™, 
Ai Na 

1 

—=252.4cm—!, —=1206cm—. 
Ao As 

1 

—= 549.5 cm“, 

As 


Representation I: Because of the choice of 
potential function, the values obtained for the 


2 bands belonging to this representation are 


1 
—=0cm™", —=0cm". 
Ai Ae 


Representation I: The values obtained for 
the 4 bands belonging to this representation are 


1 1 
—=0cm™, —=450.5 cm, 
A1 Xs 
1 1 
—=0cm™, —=1103 cm". 
Ne N4 


MOLECULES BELONGING TO THE 
POINT GROUP Dz, 


2-2-3-3-Tetramethylbutane (Table VII, Fig. 7) 


Representation !™: The values obtained for 
the 3 bands belonging to this representation are 


1 1 
—=247.4cm—!, —=1284cm—. 
Ai Xs 

1 
—=639.4cm—!, 

Ae 


Representation I: Because of the choice of 
potential function, the value obtained for the 1 
band belonging to this representation is 


1 
—=0cm—. 
Ai 


Representation I: The values obtained for 
the 2 bands belonging to this representation are 


1 1 


—=484.1cm~!, —=847.0 cm“ 
Ai Ae 


TABLE V. Table of characters. 











Cov E C2 Oz oy Raman _ red No 
TOCA, 1 1 1 Lb £ A 6 
TO~B, 1 —1 1 ' Br aA .) 
Te~A>» 1 1 i =] DF if 3 
To~B, 1 —-1 1 -1 DP A 4 
OR 0 Tv 0 0 
UR 8 2 2 4 
+1+2 cos or 3-1 1 1 
x(R) 24 —2 2 4 
a 18 0 2 4 











val 
rep) 


val 
rep’ 


wel 
oct 
bar 


the 


To), 
T®). 


rT 


for 
are 


for 
ire 


or 
re 


I 


mW | Ss 
Is 
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Representation I (doubly degenerate): The 
values obtained for the 3 bands belonging to this 
representation are 


1 1 
—=353.6cm-!, —=1270cm—. 
Ai As 

1 

—=570.7 cm“, 

Ae 


Representation I“ (doubly degenerate): The 
values obtained for the 3 bands belonging to this 
representation are 


1 1 
—=209.5cem—!, —=1123 cm. 
Ai As 

1 

—=657.4cm™, 

Ae 


SUPPLEMENTARY EVALUATIONS 


The frequencies for m-butane and n-hexane 
were evaluated by the method described for the 
octanes for the purpose of determining how these 
bands merge with those of n-octane. 


N-Butane (Table VIII, Fig. 8) 


Representation [™: The values obtained for 
the 3 bands belonging to this representation are 











1 1 
—=359.6cm—!, —=1092 cm. 
Ai As 
1 
—=979.7 cm~', 
de 

TABLE VI. Table of characters. 

Infra- 

Cov E Ce Oz gy Raman _ red No. 
TOYA, 1 1 1 , £ A 7 
TO~B, 1 —-—1 1 i BP a 5 
TO ~As 1 1 1 -—-1 DP IA 2 
TO~Bs 1 —-—1 1 —] DP A 4 
oR 0 T 0 0 
Ur 8 2 4 6 
+1+2 cos or 3 —1 1 1 
x(R) 24 —2 4 6 
= 18 0 4 6 








Representation [': Because of the choice of 
potential function, the value obtained for the 1 





Fic. 6. Skeletal form of 2-methyl-3-ethylpentane. 61° = 82° 
=63°=63°= 107°. Ba° =Bs°=Be° =B7° = 36:35". 


band belonging to this representation is 


1 
—=0cm". 
Ai 


Representation [!: The values obtained for 
the 2 bands belonging to this representation are 


1 1 
—=294.3cm-!, —=982.1 cm—. 
Ai Ae 


N-Hexane (Table IX, Fig. 9) 


Representation [™: The values obtained for 
the 5 bands belonging to this representation are 


1 1 
—=261.6cm™, —=1034cm~}, 
Ay N4 

1 1 
—=348.0cm—!, —=1119 cm". 
Ae As 

1 

— = 987.4 cm—!, 

Az 


Representation ['@: Because of the choice of 
potential function, the value obtained for the 1 
band belonging to this representation is 


1 
—=0cm~"'. 
At 


Representation [': Because of the choice of 
potential function, the values obtained for the 2 
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bands belonging to this representation are 


1 
—=0cm™, —=0cm". 
Ai Ae 
Representation [': The values obtained for 
the 4 bands belonging to this representation are 


1 1 
—=149.1em—', —=972.4cm—', 
At Az 


1 1 
—=427.00cm—, —=1092 cm—. 


Ao Na 
SUMMARY AND ANALYSIS OF RESULTS 
Internal Consistency of the Theoretical Spectra 


The calculated bands for all the molecules fall, 
roughly, into two main regions, one below 500 
cm and the other between 800 cm~ and 1300 
cm, although there is some intermixture as 
well. The bands in the higher region may be 
assigned to modes in which the motions of the 
masses are mainly along the valence bonds, 
resulting in changes in the lengths of these bonds, 
while the bands in the lower region may be 
assigned to modes in which the motions of the 
masses are mainly perpendicular to the valence 
bonds, resulting in changes in the valence angles. 
This is borne out by the evaluation by the 
discussed method of the skeletal frequencies of 
the ethane and n-propane molecules. In the case 
of ethane, the motion for the frequency is a pure 
stretching of the valence bond, and gives rise to 
a band at 965 cm™, and, for n-propane the 
motion*in"which the angle between the valence 
bonds is changed gives rise to a band at 375 cm™. 
On the whole, the calculated spectra show no 


very distinct systematic displacement of the 
bands. (See Chart I.) 


CHARLES O. 


AHONEN 


CH 
33 





7. Skeletal form of 2-2-3-3-tetramethylbutane. 
Br° = Bs° = Bo° = 810° 


FIG. 
Bi° = B2° = B3° = By° = Bs° = Be° = 107°. 
= Bi1°=Bi2° = 111°50’. 


From the spectra of the plane zigzag chain 
molecules, n-butane, n-hexane, and n-octane, 
may be seen the types of shift that result when 
new members are added to the chain of carbons. 
The single band at 965 cm7 in the spectrum of 
ethane, corresponding to a stretching motion of 
the valence bond, becomes multiple in the others 
in a somewhat irregular fashion, and, starting 
with n-butane, bending vibrations appear at 
lower frequencies, which again spread out into 
several more bands in the spectra of the longer 
chains. These results are to be expected from 
general considerations, the additional bands re- 
sulting from the addition of more elements to 
the chain. 

The main conclusion to be drawn is that small 
changes in form, and possibly also small changes 
in the potential function constants, result in 
large variations in the positions of the bands. 
Theory would lead us to expect no very system- 
atic pattern for the bands in the spectra of the 
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TABLE VIII. Table of characters. 








Inira- 
i Raman _ red 


a IA 
1 DP IA 
1 
1 





F A 
F A 





octanes. It becomes evident that comparison of 
experimental with theoretical spectra will be 
difficult, and that for the most part only quali- 
tative correspondence may be expected. An exact 
agreement with experiment would, of course, be 
extremely unlikely, since the exact forms of the 
molecules are not known, the forms used having 
been somewhat arbitrarily chosen to give the 
most symmetry possible consistent with chemical 
evidence. Also, the valence force model is itself 
an approximation. 


Correlation with Experimental Spectra’ 


Any attempt to establish correlation of the 
theoretical spectra with spectra found experi- 
mentally meets with certain difficulties. The 
observed Raman and infra-red spectra are far 
from complete. The infra-red spectra do not 
include frequencies much below 700 cm, and 
the Raman spectra, although somewhat more 
complete as to range covered, are available only 
for a few molecules. The purity of the experi- 
mental material is often questionable, especially 
in the case of the Raman spectra, so that ex- 
traneous bands due to impurities may be present. 
Both the Raman and the infra-red spectra show, 
besides the fundamental vibrations, many har- 
monics and overtones, as well as combination 
bands. Usually the fundamental bands should 
be stronger than the overtones or harmonics, but 
exceptions have been known. Acetylene is an 
example of a case in which the harmonic is very 

’ The infra-red spectra for all molecules are given by 
R. A. Oetjen, “The infrared spectra of the isomeric 
octanes,” Thesis (University of Michigan, 1941). The 
Raman spectra for 2-5-dimethylhexane and 3-4-dimethyl- 
hexane are given by J. W. Murray, J. Chem. Phys. 2, 618 
(1934). The Raman spectrum for n-octane is given by K. 


W. F. Kohlrausch, Der Smekal-Raman-Effekt (Verlags- 
buchhandlung, Julius Springer, Berlin, 1938), p. 154. 


pronounced, being stronger than the funda- 
mental. In addition, the observed spectra show 
not only the skeletal frequencies (i.e., frequencies 
arising from carbon vibrations), but also those 
resulting from the hydrogen vibrations. In most 
cases the hydrogen bands are readily identified, 
because of their position, and because of the fact 
that they repeat from molecule to molecule in 
the same pattern. For example, in the region 
covered by the charts, there are two infra-red 
bands at about 1380 cm™ and 1470 cm™, and 
three Raman bands at about 1280 cm~, 1340 
cm-!, and 1460 cm=, which clearly repeat in all 
the spectra, and so may be classed as very likely 
caused by hydrogen vibrations, and excluded 
from further discussion. (See Chart IT.) 

On the whole, the observed skeletal frequencies 
show the same type of irregularity as is found 
for the calculated frequencies, and the general 
groupings are qualitatively alike. The lack of 
pattern in the experimental bands may be 
regarded as a definite correlation between experi- 
ment and theory, even though of a negative 
character. The absence of trends in the experi- 
mental spectra of the octanes might have seemed 
somewhat surprising, since in certain classes of 
spectra a particular atom and bond are known 
to give rise to a characteristic frequency, which 
persists even when terms are added to the mole- 
cules. This is observable in the case of the C—H 
vibrations, as well as in cases where a double 
bond occurs, such as C=C, the so-called ethyl- 
enic linkage, for example, the frequency for which 
persists throughout a series of compounds. 

Possible correlation of bands in the theoretical 
and experimental spectra of the isomeric octanes 
for which all normal skeletal frequencies were 
evaluated is indicated in Chart II. The identifi- 
cation of lines is merely tentative. The following 
observations may be of interest. 

N-octane.—The experimentally observed skele- 
tal vibrations appear to fall into two main 


2 
Ay 
1 3 
CHs CH, 


Fic. 8. Skeletal form of n-butane. 81° =82° = 107°. 
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TABLE IX. Table of characters. 
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groups, corresponding to stretching and bending 
motions, as predicted. There seems to be some 
correlation of fundamental lines in the two spec- 
tra, as well as a possibility of establishing a 
correlation for the overtones in the infra-red. 
The two lines at 344.0 cm and 396.6 cm in 
the theoretical spectrum are repeated at 723 
cm and 768 cm in the experimental spectrum, 
with decreased intensity. Whether this is an 
actual case of correlation or not is uncertain. 


Infra-red Calculated 
152.6 


257.9 
372.6 
426.8 
883.7 
1061 
1124 
Raman Calculated 


217.9 


367.1 
376.0 


851.3 
1032 


1113 
1146 


Observed 


Not measured 


750 
922 
955 
1046 
1091 
1173 
1229 


Observed 


265 
313 
444 


778 
836 
962 
1048 
1150 
1176 
1302 


2-2-3-3-tetramethylbutane.—Possible correlation 





Infra-red Calculated Observed : j z 
86.16 of lines in the theoretical and experimental 
344.0 Not measured spectra is given in the table below, including 
396.6 overtone bands, as indicated. 
2(344.0) 723 
2(396.6) 768 Infra-red Calculated Observed 
882 209.5 
915 484.1 Not measured 
986.3 970 657.4 
996.7 788 
1104 1081 847.0 796 
1140 803 
2(484.1) 931 
Raman Calculated - Observed 999 
200.0 250 1123 1184 
224.3 
443.4 420 4-methylheptane, 3-4-dimethylhexane, 2-3-4-tri- 
500 methylpentane, and 2-methyl-3-ethyl pentane.—The 
810 frequencies fall into two main groups, corre- 
ta sponding to stretching and bending motions, as 
955.7 950 before. Identification of lines is extremely uncer- 
1030 tain, but the correlation indicated in Chart I] 
1053 1060 is regarded as possible. 
1054 1070 
2-5-dimethylhexane.—The same type of group- fe Ay 
ing of bands is observable in the theoretical and 1 A . A, : 
experimental spectra, in the regions covered by %"3 CH, CHa 


both, and there is a comparatively good possi- a ee OE! Loe nay 
bility for correlation of lines. 1° = 2° = 83° = By° = 107°. 
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Ethane 
n-Butane 
n-Hexane 
n-Octane 

4-Methylheptane 
2-5-Dimethylhexane 
3-4-Dimethylhexane 
2-3-4-Trimethylpentane 
2-Methyl-3-Ethylpentane 


2-2-3-3-Tetramethylbutane 





500 


Cuart I. Shows theoretical values of 1/X in cm~. Infra-red active ?. Raman active |. Active in both }. 


The agreement between theory and experi- 
ment is not very good, but in view of the fact 
that small changes in form, and possibly also 
small changes in the potential function constants, 
result in considerable shifts in the location of the 
bands, it is, perhaps, what might be expected. 
The lack of pattern in both the theoretical and 
experimental bands is a definite correlation. 


n-Octane 

Experimental 
4-Methylheptane 
Experimental* 
2-5-Dimethylhexane 
Experimental 
3-4-Dimethylhexane 
Experimental 
2-3-4-Trimethylpentane 
Experimental* 
2-Methyl-3-Ethylpentane 
Experimental* 
2-2-3-3-Tetramethylbutane 
Experimental* 


* Raman Spectra not measured 


The tentative correlation established by this 
research is valuable, since it can give the type of 
motion belonging to the observed spectral bands, 
and this knowledge may be of considerable value 
in interpreting the spectra of the derivatives of 
the octanes. For example, if it is known that a 
certain band corresponds to a given motion of 
the atoms of an octane, the disappearance of 





Cuart II. Shows theoretical and experimental values of 1/A in cm™. Infra-red active }. Raman 
active |. Active in both j. 
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this band in a derivative of this octane would 
indicate that the atomic group whose motion is 
responsible for the band has been changed by 
the added terms, while the persistence of the 
band would indicate that the added terms are 
not connected to that group in such a way as to 
interfere with its motion. 
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sor J. C. Halford, of the Department of Chemis- 
try, University of Michigan, for his collaboration 
in deciding on the forms that might best repre- 
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1, INTRODUCTION 


ECENTLY some articles have been pub- 
lished dealing with the behavior of polymers 
of the coiling type suspended in a flowing liquid.! 
One of the issues is to establish a physical reason 
for Staudinger’s rule according to which the 
intrinsic viscosity of polymer solutions should be 
proportional to the molecular weight of the 
polymer. I would like to present here a detailed 
picture of the motion of the polymer molecule 
and its parts, which makes it easy to see why the 
rule should hold for the very simplified model 
currently substituted for the actual molecule. 


2. MOTION OF THE POLYMER MOLECULE 


The model is as usual a number of identical 
atomic groups connected to each other by rigid 
links, the consecutive ones making a definite 
angle with each other and able to turn freely 
around each other. It is also assunied that each 
of these groups is acted upon by a frictional 
force fv if in a solvent flowing along with the 
velocity v and that neither the links suffer friction 
from the medium nor any interaction between 
the different groups exists, due to their individual 
disturbance of the general flow. 

The velocity of the liquid in a rectangular 
x, y, system of coordinates may be represented 


1 John G. Kirkwood, J. Chem. Phys. 14, 180, 347 (1946); 


H. A. Kramers, J. Chem. Phys. 14, 415 (1946), a sequence 
to an article by J. J. Hermans, Physica 10, 777 (1943). 


by the vector V with the components 


aZ=V, 


V=) 0=1, (1) 
0=v, 


as it would result for instance from the motion of 
a plate in the x,y-plane moving with a constant 
velocity in the x-direction at a constant distance 
above another parallel plate at rest. The velocity 
gradient in the z-direction is a. 

A polymer molecule suspended in this liquid 
will be set in motion by the frictional forces 
acting upon each atomic group of the chain. If 
now we confine our attention to the center of 
gravity of the molecule, it can easily be seen 
that this center of gravity (apart from its 
irregular Brownian motion) will move with the 
velocity of the liquid as it is at the position 
where this center is momentarily situated. In 
this way we have taken care of the condition 
that the sum of all the forces acting on the whole 
molecule, averaged over all the configurations 
the chain can acquire is zero. 

Let us now take the center of gravity as the 
origin of our system of coordinates. Equation (1) 
will then represent the relative velocity of the 
liquid with respect to this center and will appear 
as represented in Fig. 1. It is evident from Fig. 1 
that the still remaining frictional forces will now 
tend to rotate the molecule as a whole around 
the y-axis. Call the angular velocity of such a 





roti 
due 
con 


anc 
wil 
hay 
mo 
gra 
has 


the 
the 
xy 


in 
cer 
ha’ 
wil 


In 
be 
rot 


Bu 


Fo 


ha 


as 


th 


tment 
egest- 
ssions 
rofes- 
iemis- 
ration 
repre- 
rated. 


1946 


(1) 


on of 
stant 
tance 
ocity 


iquid 
orces 
n. If 
er of 
seen 
1 its 
1 the 
ition 
1. In 
ition 
rhole 
tions 


; the 
n (1) 
' the 
pear 
ig. 1 
now 
yund 
ch a 


VISCOSITY OF POLYMER-SOLUTIONS 


rotation w. A group at the position x, y, z will, 
due to this rotation, have a velocity v with the 
component, 
ws 
yv=40 (2) 


— Wx 


and the frictional force acting on such a group 
will be F=f(V—v). In the steady state we will 
have to satisfy the condition that the average 
moment of this force around the center of 
gravity is zero. The instantaneous moment M 
has the components 


fw Do xy 
M=4fLl(a—w)  2—w Dd x"), (3) 
—f(a—w) Li yz 


the summation extending over all the groups of 
the molecule. Since in the average the products 
xy and yz are zero and 


r2 
Av (x?) = Av (2?) =Av (-). 


in which r is the distance of the group from the 
center of gravity, the average moment will only 
have a component in the y-direction and this 
will be equal to 


Av (M,) = f(a—2w) Av (z =). (4) 


In the steady state this average moment has to 
be zero and so we conclude that the molecule 
rotates as a whole with an angular velocity 


w=a/2. (5) 


But even so the liquid will still have relative 
velocities with respect to the different groups. 
For one of them this relative velocity V—v, will 
have the components 





as follows from (1) and (2) with w=a/2. 
The flow lines of this relative flow of the liquid 
through the skeleton of the molecule are repre- 





ei 








Fic. 1. Relative velocity of the medium with respect to 
the center of gravity of the molecule before rotation. 


sented in Fig. 2. The frictional energy-losses con- 
nected with this flow will appear in the form of 
an increased viscosity of the solution. 

The flowlines themselves are hyperbolas with 


the equation 
2° —x* =const. 


and the absolute value of the velocity 
a 
—(2?+x?)! 
2 


is constant on circles around the y axis, in- 
creasing proportionally with the distance from 
this axis. 


3. ENERGY LOSSES 


In the flow represented by Fig. 2 and Eq. (6) 
the frictional force on one of the groups is 


a 
f-(2?+x?*)! 
2 


and since this force is in the direction of the 
velocity, the work in unit time for one group is 


2 


a 
f—(? +x"). 
4. 


If m is the number of molecules per cc, the 
increase in energy losses per sec. and in 1 cc 
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due to the presence of polymer molecules will in 
the average be equal to 


w=nf—Av (rc (tte) =nf— Av(Xr). (7) 


This increase makes itself felt by an observable 
increase of the macroscopic viscosity from a 
value 7 to a higher value n’. Now in a medium 
of viscosity 7 the total energy loss per sec. and 
per cc for a velocity distribution as represented 
by Eq. (1) is 

w= na’, (8) 


This is easily checked by realizing that if one of 
the two parallel plates of surface S at a distance 
h is moving with the velocity v= ah with respect 
to the other, which is at rest, the frictional force 
on the upper plate is F=naS and the displace- 
ment of this plate per unit time is y=ah. Con- 
sequently the total amount of work dissipated 
in the volume Sh per sec. is Fu=yna*®Sh. If due 
to the presence of the polymer molecules the 
medium acts as if the viscosity is increased from 
n to 7’, the corresponding increase in dissipated 
energy will be 

w= (n’—7n)a’. (9) 


Equations (7) and (9) represent the same quan- 
tity in two different forms, so we conclude that 


n’ —=6nf 2. Av (7°). (10) 


It can now already be seen at once that the 
increase in viscosity will be proportional to the 
square of the number of links contained in the 
chain of the molecule. It is well known that the 
average of the square of the distance from 
beginning to end of the chain is proportional to 
N, denoting by N the total number of links. From 
this it follows immediately that the average of 
the square of the distance of one of the groups 
from the center of gravity of the whole molecule 
will also be proportional to NV. On the other hand 
the number of terms in the sum of Eq. (10) is 
N+1. Neglecting terms of the relative im- 
portance 1/N against 1, we thus come to the 
conclusion that the sum is proportional to N?. In 
a qualitative way this dependence is caused by 
the fact that in longer molecules the average 
distance of each group from the center of gravity 
is larger. Consequently each group comes into a 


P. DEBYE 


Fic. 2. Relative velocity of the medium with respect to 
the center of gravity of the rotating molecule. 


region in which the relative velocity of the 
medium with respect to the rotating molecule is 
also larger. Increasing the length of the molecule 
therefore means not only increasing the number 
of centers which contribute to the energy losses 
but also increasing at the same time the indi- 
vidual energy loss due to each center. 


4. INTRINSIC VISCOSITY 


In order to complete the calculation indicated 
in Eq. (10), we have to know what the average 
square of the distance is of one of the groups 
from the center of gravity. Obviously this aver- 
age will depend on the position of the group in 
the chain and we expect that for a group in the 
middle it will be smaller than for a group near 
the end of the chain. This is confirmed by a 
calculation not reproduced here, which although 
somewhat more involved proceeds along the 
same lines followed in calculating the average 
square of the distance from beginning to end, 
which is 

R? = +? Nat 
1—p 


In Eq. (11) a is the length of a link and # is the 
cosine of the angle between two consecutive 
links (3 for C—C bonds). 

If we number the N+1 groups of the chain 
from O to N and concentrate our attention on 
group number », it is found that the average 
square of the distance of this group from the 


(11) 
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center of gravity is 


tad v v 
oltine’ 12)1 aa 

“'< —p 3 N N 
The factor in front of the brackets being equal 
to R?, we see that for y»=0 or v=WN (the ends of 


the chain) r?=4R?, whereas a minimum value is 
reached for v/N= 3 (the middle of the chain) 


equal to r?=7;R*. Substituting (12) in (10) and 
performing the summation, which within the 
limits of accuracy can be represented by an 
integral, we find 


eae ae) 


1 1+ p 
sonal HP [3-e1—g) de (13) 
6 1-—p 0 


1 1+) 
=—nf——N’*a?. - 
36 1-—p 


From the increase in viscosity we obtain the 
relative increase in dividing by 7 and the intrinsic 
viscosity [7] in dividing by the concentration, 
which is » Nm if m is the mass of one of the groups 
of the chain. The final result for the intrinsic 
viscosity is 

1 f/nit+p 


[n]=—— ——Na’, (14) 
36 m 1—p 


which represents Staudinger’s rule that the 
intrinsic viscosity is proportional to N and 
therefore to the molecular weight of the whole 
molecule. Should we represent each group by a 
sphere of radius p and apply Stokes law, ac- 
cording to which in this case f=67np, we will 


find 
aw pR? 
[yjJ=—-—. (14’) 
6 m 

As is evident from the definition of [7] it has 
the dimension of a specific volume and this is 
according to the theory represented here by the 
quotient of a volume which is 7/6 times the 
product of the radius of the sphere, representing 
each group multiplied by the average square of 
the distance from beginning to end of the chain, 
divided by the mass of a single group. In order 
to show. by an example that the order of mag- 
nitude is correct we take a polymer of molecular 
weight 500,000, each single group having the 
molecular weight 50. We now have N=10,000 
and with a=1.54A (the C —C distance) it follows 
that R?=47,400.10—* cm?. Since m=82.10-* g 
we find that in order to explain a supposedly 
observed intrinsic viscosity ['y ]= 100 cc/g (which 
in the customary practical units, 1 g and 100 cc, 
corresponds to the intrinsic viscosity (1) it is 
necessary to represent each group by a sphere 
of radius p=0.33A. 

As to the streaming double refraction, which 
may be expected from the model, it is immedi- 
ately evident from Fig. 2 that in a second 
approximation the molecule will in the average 
lose its central symmetry, be elongated in the 
direction AA, and compressed in the direction 
BB. With bonds having different polarizabilities 
parallel and perpendicular to their direction this 
will lead to observable double refraction of the 
solution in motion with axes under 45° against 
the direction of flow. At the same time it is also 
evident with respect to the viscous behavior of 
the solution that in a second approximation and 
due to the distortion of the molecule propor- 
tionality between the shearing stress and the 
velocity gradient will cease to exist. 
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The Effect of Solvents on the Critical Concentra- 
tion for Micelle Formation of Cationic Soaps* 


M. L. CorrIn AND WILLIAM D. HARKINS 
The University of Chicago, Chicago, Illinois 
August 21, 1946 


HE effect of various solvents upon the critical con- 
centration for micelle formation has been discussed 
by Ward! and by Ralston and Hoerr.? The development of 
titrametric methods for the determination of critical con- 
centrations has made it practical to extend such studies.* 
With the titration technique it is possible to measure a 
number of critical concentrations under varying conditions 
far more easily than if a conductometric method is. used; 
by this latter procedure a large number of accurate meas- 
urements are required for a single determination. The 
work of Ward is restricted to sodium dodecyl sulfate, 
while Ralston and Hoerr consider only dodecylammonium 
chloride; the latter workers, however, do not report any 
values for the critical concentration. The work herein 
reported is restricted to cationic soaps since no indicator 
is available which gives a sharp end point with anionic 
materials in solutions containing organic solvents. 
Data for the critical concentration of dodecylammonium 
chloride in various solvent-water mixtures are plotted in 
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Fic. 1. The critical concentration for micelle formation in solutions of 
dodecylammonium chloride in various solvent-water mixtures. 
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Fig. 1; the indicator was dichlorofluorescein. With decy]! 
trimethyl ammonium bromide the measurements were ex- 
tended in the region of low solvent concentration. These 
data are plotted in Fig. 2; the indicator employed was 
Sky Blue FF. 

With decyl trimethyl ammonium bromide the decrease 
in critical concentration becomes more pronounced as the 
chain length of monohydric alcohol used as the solvent is 
increased. Dioxane also produces a slight initial decrease 
in critical concentration while no decrease was observed 
with ethylene glycol even at low solvent concentrations. 
A similar behavior was observed with dodecylammonium 
chloride. Isopropanol is but slightly different in its effect 
from n-propanol. At higher solvent concentrations the 
critical concentration of the detergent increases with 
dioxane, ethylene glycol, glycerol, and ethanol. At the 
highest concentrations of isopropanol, u-propanol, and 
n-butanol employed, the critical concentration was still 
decreasing. 

On the basis of the above results it seems probable that 
Ward's hypothesis, which attributes the initial drop in 
critical concentration to adsorption of solvent molecules 
on the “surface of the micelle’’ is untenable. 

It is, at present, difficult to correlate experimental 
observations with any dielectric constant data for the 
various solvent pairs. This is not unexpected since Ward 
has pointed out that a significant term in micelle formation 
is the dielectric environment about the polar head groups 
of the detergent molecules in the micelle; this quantity 
cannot, at the present time, be measured or even approxi- 
mately calculated. 

The term “interfacial energy”’ applied by many authors 
to the energy arising from an interaction of amphipathic 
molecules with surrounding solvent molecules is mislead- 
ing and tends to obscure the actual processes involved in 
such systems. Interfacial energy is both defined and meas- 
ured as the energy difference between the surface and bulk 
of a substance; there is no reason to consider that the 
interfacial energy represents the interaction energy of a 
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Fic. 2. The critical concentration for micelle formation in solutions of 
decyl trimethyl ammonium bromide in various solvent-water mixtures. 


640 





1946 


decyl! 
re ex- 
These 
1 was 


crease 
is the 
ent is 
crease 
erved 
tions. 
nium 
effect 
s the 
with 
t the 
and 
3 still 


» that 
op in 
cules 


ental 
r the 
Ward 
ation 
oups 
ntity 
rOXxi- 


thors 
athic 
lead- 
ed in 
neas- 
bulk 
> the 
of a 


LETTERS TO 


single molecule or portion of a molecule with the surround- 
ing solvent molecules. 

Since it is impossible at present to evaluate the energy 
terms involved in micelle formation, the experimental 
data herein reported can contribute only to the empirical 
knowledge of the properties of detergent solutions. It is 
hoped that the accumulation of further experimental 
information on a wide variety of properties of such solu- 
tions may eventually lead to a true understanding of 
their nature. 

* The work reported in this paper was done in connection with the 
Government Research Program on Synthetic Rubber under contract 
with the Office of Rubber Reserve, Reconstruction Finance Corpora- 
nA. F, H. Ward, Proc. Roy. Soc. A176, 412 (1940). 


2 A. W. Ralston and C. W. Hoerr, J. Am. Chem. Soc. 68, 851 (1946). 
3M. L. Corrin and W. D. Harkins, following letter. 





Determination of Critical Concentrations for 
Micelle Formation in Solutions of Cationic 
Soaps by Changes in the Color and 
Fluorescence of Dyes* 


M. L. Corrin AND WILLIAM D, HARKINS 
University of Chicago, Chicago, Illinois 
August 26, 1946 


N the first paper! on the use of dyes for the determina- 

tion of the critical concentration for micelle formation, 

a titration method developed by the authors was men- 

tioned. In this work with anionic soaps a cationic dye, 

pinacyanol chloride, was used. Attempts to use this same 
dye with cationic soaps were unsuccessful. 

The following dyes, whose spectra in alcohol are visually 
distinguishable from those in water, were found to be 
ineffective with either anionic or cationic soaps: Orange G 
(S189), Azophloxin 2G (S40), Bordeaux B (S123), Congo 
Brown G (S280), New Methylene Blue (51043), and 
Ponceau G (S37). 

It was thought desirable to attempt to find dyes, one 
of whose several molecular forms is preferentially solu- 
bilized in soap micelles. This involved an extensive investi- 
gation in which it was discovered that the following dyes 
may be used as indicators with cationic soaps: Sky Blue 
FF, Eosin, Fluorescein, Dichlorfluorescein, and acidified 
2,6-dichlorophenol indophenol. Another dye other than 
pinacyanol chloride applicable to anionic soaps, Rhoda- 
mine 6G, was also found. The values of critical concentra- 
tions determined by the use of these dyes are quite self 
consistent, as shown in Table I. 

These results are also in agreement with literature values 
obtained by other, much more laborious methods. The 
critical concentration for dodecyl ammonium chloride is 
0.0130 molar, as reported by Ralston and Hoerr,? and for 
decyl trimethyl ammonium bromide 0.065 molar by Scott 
and Tartar.? 

Critical Concentration by Change in Fluorescence.— 
Cationic soap solutions containing eosin or fluorescein are 
highly fluorescent when micelles are present. This fluores- 
cence is strongly quenched at concentrations below the 
critical; hence this property, as well as color, serves to 
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TABLE I. The titrametric determination of critical concentrations 
with various dyes. 








Anionic Soaps 
A. 


Cationic Soaps 
A. 


_ 
Decyl 
Dodecyl trimethyl 
ammonium ammonium 
chloride bromide 





~ 


Sodium Sodium 
Potassium dodecy! decyl 
Dye laurate sulfate sulfonate 
Pinacyanol 0.0235 0.00602 0.0400 
chloride 
Rhodamine 6G 0.0234 0.00612 0.0387 
Sky Blue FF 
Eosin 
Fluorescein 
2,6-dichloro- 
phenol- 
indophenol 











determine the critical concentration, With eosin the color 
is orange above and red below the critical concentration. 
Both of these dyes are anionic and are effective with only 
cationic soaps. 

Rhodamine 6G, a cationic dye and effective only with 
anionic soaps, also yields orange fluorescent solutions in 
the presence of micelles. When the soap concentration is 
reduced, the color changes to red and the fluorescence 
almost entirely disappears at the critical concentration. 

* The work reported in this paper was done in connéction with the 
Government Research Program on Synthetic Rubber under contract 
with the Office of Rubber Reserve, Reconstruction Finance Corpora- 
nM. L. Corrin, H. B. Klevens, and W. D. Harkins, J. Chem. Phys. 14 
480 (1946). 


2A. W. Ralston and C. W. Hoerr, J. Am. Chem. Soc. 64, 772 (1942). 
3A. B. Scott and H. V. Tartar, J. Am. Chem. Soc. 65, 692 (1943). 





On the Structure of Ozone 


W. H. EBERHARDT 
Georgia School of Technology, Atlanta, Georgia 
September 9, 1946 


SING central-force-field methods, Glockler and Mat- 
lack! have recently presented a discussion of the 
bond energies and molecular structure of ozone based on 
a correlation of bond energy and interatomic distance 
suggested by other oxygen-oxygen bonds. Since their dis- 
cussion is limited to acute angle models of the molceule, 
it is interesting to extend the technique to obtuse angle 
models. 

From electron diffraction results? the ozone molecule 
has been reported to be an obtuse isosceles triangle with 
O—O bond length 1.26+0.02A and bond angle 127+3°. 
The separation of the base oxygen atoms is thus 2.25A. 
If one assumes the dependence of bond energy on inter- 
atomic distance as suggested by Glockler and Matlack, the 
energy of the two valence bonds is 2X 3.06 =6.12 electron 
volts. The contribution of the base O—O bond to the 
total binding energy will be small, and indeed may well 
be negative owing to the repulsion of the electron clouds 
of the two base oxygen atoms. Hence, it may be con- 
sidered that adequate agreement with the total binding 
energy, 6.12 electron volts, may be obtained on the basis 
of this model as well as acute angle models. 


1 Glockler and Matlack, J. Chem. Phys. 14, 505 (1946). 
2? Shand and Spurr, J. Am. Chem. Soc. 65, 179 (1943). 
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Adsorption Isotherms for Hydrogen-Bonded 
Liquids on Charcoal 


D. H. VoLMAn* AnD I. M. Kiotz 
Northwestern University, Evanston, Illinois 
August 28, 1946 


N van der Waal’s adsorption on charcoal, a plot of the 

weight adsorbed as a function of the relative pressure 
(p/po) of the vapor gives a curve which is roughly the 
same for a large number of substances. A few exceptions 
are known, of which water is the outstanding example, 
in which marked deviations from the typical curve occur, 
particularly at low pressures. The underlying basis of 
these deviations becomes evident when it is realized that 
these exceptional substances all form strong hydrogen 
bonds in the liquid state. A stronger parallelism is then 
found to exist between adsorbability and condensation 
properties of a gas than has been realized previously. 

A plot of boiling point as a function of van der Waal’s a 
shows a region, bounded by the curves in Fig. 1, in which 
“normal” liquids are found. These are the substances 
which exhibit normal adsorption isotherms. In contrast, 





ae 


? 


“NORMAL” 
Liauios 











THE EDITOR 


water, ammonia, and methyl alcohol, which show abnormal 
isotherms, also stand outside the normal region in Fig. 1, 
the distance being a rough measure of the deviation from 
the typical isotherm. 

The deviation in the boiling-point curve is caused by the 
presence of strong hydrogen bonds in the liquid state, a 
condition which necessitates more energy to vaporize the 
molecules than would be required to overcome van der 
Waal’s forces alone. 

The deviations in the adsorption isotherm are also caused 
by hydrogen-bonding which may express itself in two 
different ways. In the region of low relative pressures 
(0.3-0.4), for water for example, the extent of adsorption ap- 
pears to be unusually small. The apparent anomaly, is ow- 
ing, however, to misleading values assigned to the relative 
pressure when one deals with a strongly hydrogen-bonded 
liquid. If there were no such bonding in liquid water, the 
boiling point and vapor pressure would be in line with 


- those of the other hydrides of the sixth-group elements 


and would be near — 80° to — 100°C, and about 20,000 mm 
(at room temperature), respectively. On the basis of this 
latter value for the vapor pressure if the liquid had no 
hydrogen bonds, the relative pressure of 0.3-0.4 is reduced 
to 0.0003-—0.0004 and the observed adsorption is normal 
even though small. Since the adsorption at low pressures 
on charcoal (not containing specially-prepared, surface 
complexes) does not reflect in any way the hydrogen- 
bonding properties of water, this vapor behaves as if it 
came from a normal liquid with a vapor pressure of 20,000 
mm. 

Above relative pressures of 0.3-0.4 there is a steep rise 
in the water isotherm. Clearly this must be owing to the 
formation of a network of hydrogen-bonded water mole- 
cules, anchored initially by the few adsorbed molecules. 
Under these conditions the adsorbed phase does reflect 
the properties of the liquid and hence the isotherm rapidly 
approaches the normal curve, 


* On leave, 1941-45, from University of California, Davis, California. 





